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The Triangulation of a Two-Dimensional Continuum for the Purpose 
of the Approximate Solution of Second-Order Partial 
Differential Equations 


L. Tasny-TSCHIASSNY 
Department of Electrical Engineering, University of Sydney, N.S.W., Australia 
(Received April 26, 1948) 


It is shown how a non-homogeneous, non-isotropic, electrically conductive, two-dimensional con- 
tinuum can be replaced by a network of resistors whereby the nodes of the network are situated at 
any arbitrarily selected points, and how Southwell’s relaxation method can be applied to this network 


for the solution of potential problems. 


I. THE GENERAL DIFFERENTIAL EQUATION 
CONSIDERED 


ECON D-ORDER partial differential equations 
for the unknown function V=V(x, y) of the 
general form 


df AV avy a av aV 
—(1.—+ r9—) + Yat w—) 
Ox Ox oy oy Ox oy 


—i(x,y)=0, (1) 


where yz, Yzy, Yy, and i(x, y) are given functions 
of x and y, can be interpreted by the analog of a 
plane non-homogeneous and non-isotropic elec- 
trically conducting sheet, the conductivity of which 
(given by the quantities yz, Yzy, Yy) 1S a two- 
dimensional symmetrical tensor varying from point 
to point, and the transverse loading of which by 
an external electric current density leaving the 
conducting sheet is represented by the function 
i(x, y) which varies from point to point. V(x, y) is the 
electric potential and —[y.(@V/dx)+y2,(dV/dy) | 
and —[y.,(0V/dx)+y,(9V/dy)] are the compo- 
nents in the x and y directions of the current density. 
Special cases of Eq. (1) are the quasi-plane- 
harmonic equation (yz;=y, and y,,=0), the two- 
dimensional Poisson equation (y.=y,=constant, 
Yx=0), and the two-dimensional Laplace Equa- 
tion (y,= y,=constant, yz,=0, i(x, y) =0). 
Methods of solving Eq. (1) approximately for 
given boundary conditions are very important for a 
large variety of mechanical and electrical problems. 


Il. THE REPLACEMENT OF A DIFFERENTIAL 
EQUATION BY A DIFFERENCE EQUATION 


When solving a differential equation for the 


function V(x, y) approximately as a difference 
equation, one concentrates on the values of the 
function V(x, y) at a sufficiently large number of 
specific points ” or s, called in the following the 
nodal points, and tries to find a set of values V, or 
V, for V(x, y) at these points that complies 

(a) with a difference equation, i.e., an equation that ex- 
presses a relation between the values of V(x, y) at a certain 
nodal point and some neighboring nodal points, 

(b) with the given boundary conditions. 
The difference equation is derived from the differ- 
ential equation by expressing the first derivative of a 
function as the quotient of the function values and 
the argument values at two neighboring points, 
and by applying this approximation to the higher 
derivatives in considering them as the first deriva- 
tive of the derivative of one degree less. 

Southwell’s relaxation method is a modern and 
very efficient method for finding a set of values of 
V(x, y) that complies with the difference equation.* 

Southwell pointed out that Prandtl’s Mem- 
brane** analogy when used as a physical representa- 
tion of the quasi-harmonic, Poisson, and Laplace 
equations, i.e., special cases of Eq. (1), may be con- 
sidered as being physically replaced by a network 

*R. V. Southwell, Relaxation Methods in Theoretical Physics 


(Clarendon Press, Oxford, 1946). 
** See reference 1, pp. 39, 86, 147. 
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Fic. 1. Net of triangles of different shapes and sizes. 


of strings connected at the nodal points, if the 
differential equation is replaced by the difference 
equation. In the conducting sheet analog used in 
this paper the conducting sheet may be considered 
as being physically replaced by a network of electric 
resistors joined at the nodal points. Expressed in 
terms of the conducting sheet analog the approxi- 
mate solution of a differential equation by a differ- 
ence equation may be described in the following 
way: 

(1) Replace the two-dimensional continuum by a 
network of electric resistors. The value of the con- 
ductance G,, of the resistor connecting the two 
nodes s and m depends on the type of the difference 
equation. 

(2) Find suitable values of the transverse node 
loading currents J, that are to be placed at the 
nodes in order to simulate the loading of the con- 
tinuum by the continuous current density i(x, y). 

(3) Consider the potentials V, and V, of the 
nodes as unknowns, work out for each node s the 


expression 
T,= > n(Vn— Vs) -Gas, (2) 
and put J, equal to J, for each node; 
I,=1.,. (3) 


The equations are simultaneous linear equations for 
the unknowns V,. Since many of the coefficients 
Gns are equal to nought, the equations are particu- 
larly suitable for solution by approximate methods. 

As far as the author is aware, previous methods 
have used only arrangements of the nodal points 
in regular patterns of triangles, squares, and 
hexagons.*** Certain artifices were necessary to 
cope with the irregularities arising at the bound- 
aries. Furthermore the general form of Eq. (1) in 
which y.# 7, and y.,#0 has not been given special 
attention. 


Ill. SOUTHWELL’S RELAXATION METHOD 


Southwell’s procedure for finding an approximate 
set of values V, complying with Eqs. (2) and (3) is 
this: 

(1) Approximate values of V, are assumed for 


all nodes. 


~_ for example, the Historical Survey in reference 1, 
p. 36 
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(2) I, is computed for all nodes according to 
Eq. (2) and the values of the so-called residuals R, 
are found for all nodes according to the formula 


R,=1,—I.. (4) 


(3) The potential of one node s with a large 
residual is changed by the value AV,. The change of 
the residual AR, of this node is found from the 
formula 


AR,= —AV,: > nGen, (5) 
and of all other nodes m from the formula 


AR,=AV,°-Gas. . (6) 


The values of the new résiduals are computed. 


(4) The procedure described in (3) is repeated 
for all nodes—always considering the new resid- 
uals—until all residuals have become negligible. 


IV. THE TRIANGULATION OF THE CONTINUUM 


The possibility of the replacement of the con- 
tinuum by a network of regular pattern raises the 
question whether it would not be possible to em- 
ploy for the same purpose a network of irregular 
triangles of different shapes and sizes with any 
specified position of the nodes (see Fig. 1). Once 
the values of the conductances of the resistors con- 
necting the nodal points are known Southwell’s 
relaxation method as described in Section III can 
be applied at once. 

The advantages of an irregular pattern of nodes 
are: nodal points can be arranged for on the 
boundaries and on interlayers where the conduc- 
tivity changes discontinuously (i.e., on the border 
lines of different materials), the nodes can be spaced 
differently at different parts of the field (grading of 
the net), and finally one can arrange for different 
directions of the lines connecting neighboring nodes 
at different parts of the field. The possibility of 
placing nodal points on the boundaries and inter- 
layers eliminates all special problems connected 
with them, no matter whether the values of the po- 
tential or the values of the flux density are specified 
on the boundaries or whether the boundary condi- 
tions are mixed. 

The author found that the question raised above 
can be answered affirmatively, and it is the purpose 
of the present paper to demonstrate this and to 
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Fic. 2. Replacement of the interior of a triangle 
by a delta of resistors. 
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give the necessary mathematical formula for find- 


ing the values of the conductances between the 
nodes. 


V. TRANSFORMATION FORMULAS FOR TENSOR 
CONDUCTIVITIES 


Let (x, y) be the axes of an arbitrary system of 
Cartesian coordinates and (x’, y’) the axes of that 
system of coordinates that coincides at the point 
considered with the principal axes of the conduc- 
tivity tensor at that point. Let y be the angle be- 
tween the x’ and x axes (see Fig. 3). Let yz and 
yy be the two different principal conductivities in 
the directions x’ and y’. If i,’, iy, iz, iy denote the 
components of the current density at the point 
considered and V the potential of this point, we 
obtain 


—iz = y2+(dV/dx’), (7) 
— ly’ =Vy’" (d V/dy’). (8) 


It follows from the superposition principle and a 
transformation of the coordinates that 


; OV OV 
—tg2 Ys TY (9) 

Ox oy 

OV OV 
—ty™ Yay? ty" (10) 

Ox oy 

where 
Ye =Y¥2Cosy+ yy’- sin*y, (11) 
Yy =yz2sin*y+ yy -cos*y, (12) 
Te — Ve’ 

Yzy =————:sin(2y). (13) 


By applying Kirchhoff’s first rule to an infinitesimal 
element at the point considered we obtain Eq. (1) 
as the differential equation for the potential func- 
tion V. 


VI. THE LEADING PRINCIPLE IN THE DERIVATION 
OF THE BASIC FORMULA 


In replacing Eq. (1) by a difference equation 
certain mathematical approximations are made. In 
the cases treated in the previous literature (i.e., 
for yz=‘Yy, Yx=0 and for regular patterns formed 
by the nodal points) the replacement of the differ- 
ential equation by a difference equation corre- 
sponds to the replacement of the continuum by a 
network of resistors (see Section II). The values 
of the resistors can be found either in the usual 
way from the difference equation, or, as can readily 
be demonstrated, from the concept of a physical 
network of resistors by the following procedure: 
It is assumed that the network is free from any 
transverse loading at the nodes and is placed in a 
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Fic. 3. Position of a triangle with respect to the 
axes (x, y) and (x’, y’). 


uniform field. It is stipulated that for any intensity 
and direction of the field and for a certain width 
perpendicular to the direction of the field the 
average flux equals the actual flux produced in the 
continuum, if the continuum is not loaded by trans- 
verse current densities. 

The second procedure just described can be 
applied to the general case (yz# Vy, Y2xy#0, and for 
the nodes arranged in an irregular pattern) and is 
the method used in this paper for deriving the 
basic formula. 

When carrying out the described procedure, one 
replaces the part of the continuum inside a tri- 
angle A, Aso, A; (Fig. 2) by a delta of resistors 
arranged between its vertices. If this is repeated 
for all triangles into which the continuum is di- 
vided, one always obtains two resistors in parallel 
that connect two nodal points, at least one of which 
is situated in the interior of the field, and one 
resistor between two nodal points situated on the 
boundary. By combining into one the two resistors 
connecting internal nodes, the final network that 
replaces the continuum is obtained. 


VII. THE DERIVATION OF THE BASIC FORMULA 


Let the triangle A,A2A; (Fig. 3) be subjected to 
a uniform field in an arbitrary direction. The com- 
ponents 2, and 7, of the current density in the plane 
are given by 


; OV OV 

2 >= -(+. —+10—), (14) 
Ox oy 

OV OV 

in -(1w—+1 wt (15) 
Ox oy 


The total current J,. leaving the interior of the 
triangle through the side A1Az is given by 


: OV OV 
I}2= —3 sings: (—+10—) 
Ox oy 


+43 cos¢s: (12+ w—), (16) 
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and similar equations hold for the corresponding 
currents J23 and J3;. 

If the continuum is to be equivalent to the delta 
of resistors, the current J, leaving the delta of 
resistors at the vertex A, must be given by 


T11,=3(1i2.+Ta1), (17) 


and similar expressions hold for J, and J3. 

On the other hand, the symbols g3, gi, and ge 
denote the conductances of the resistors connected 
between the nodes A1;A2, A2A3, and A3A,, respec- 
tively. The current J; is given by the expression 


OV OV 
I= (a: COSY3:—— +43 sings: —) “£3 
Ox oy 
OV OV 
— (a: COS¢2:——+d2 sings —) ‘go, (18) 
Ox oy 


and the currents J, and J; by similar expressions. 
By substituting in (17) for Ji. and J3, from the 
equations of the type of (16) and for J; from Eq. 
(17) we obtain an equation in which a simplifica- 
tion is possible by the help of 


a; sing; = — (dz sings.+az; sings), (19) 


and @; COSY; = — (dz COSY2+43 COS¢s). (20) 
Since the resulting equation must hold for any 
values of dV/dx and dV/dy, we obtain two equa- 
tions by equating the coefficients of dV/dx and 
aV/dy on the right- and left-hand sides of the equa- 
tion, respectively. 

By eliminating g2 and making consecutive use of 
the relations 








¢2— ~3=ait7, (21) 
a a3 
a, (22) 
sina, sinag 
$i-— ¢2=as+z, (23) 








BOUNDARY A 


Fic. 4. Division of the field into two parts. 
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we obtain 


Yet Vy Yz— Vu cos(¢i+ $2) 








223= -cotas+ 


sinas 


sin(¢i+ ¢2) 
= 


sinas 


By substituting for y:, y:,, and y, from Eqs. (11) 
to (13) we obtain as final result 








V2 t+ Vy’ Ys’ — Ye’ cos(¢1' + ¢2’) bs 
223 = - cota3+— ee , (25) 
2 sina; 
where 
gi =gi-y ) (26) 
4 
go’ =92—¥ | 


represent the inclinations of the sides towards the 
principal x’ axis of the conductivity tensor. Hence 
g; is invariant. Similar equations hold for g: and go. 

Mathematical investigations not described here 
show that equations of the type of (24) or (25) 
always result whatever starting point is taken for 
finding the values of gi, go, and g; and that for an 
assembly of deltas of different sizes and shapes the 
resulting network replaces the continuum in every 
respect, as long as the field is uniform. Hence 
Eq. (25) is the basic formula in the sense of Sec- 
tion VI. 


VIII. CONCLUSIONS FROM THE BASIC FORMULA 


As was explained in Section VI, there are always 
two conductances to be considered that connect 
two nodes unless both nodes are on the boundary. 
These two conductances are given by Eq. (25) and 
will be called partial conductances. Their sum will 
be called total conductance between the two nodes 
concerned. 

Formula (25) shows that the partial conductance 
consists of two additive terms, the first depending 
on the mean conductivity (y2:+yy’)/2=(vzt Yy)/2 
and the second on the degree to which the material 
is non-isotropic, i.e., the value (yz—y,’)/2. Par- 
tial, and consequently total, conductances can be 
positive, zero, and negative. For isotropic materials 
(yz’=Yy’) the partial conductance is positive, zero, 
and negative, if the angle opposite to the side con- 
sidered of the triangle is acute, right, and obtuse, 
respectively. 

Cyclic polygons in an isotropic medium exhibit 
interesting features. Let the polygon be divided into 
triangles drawn through one of its vertices. Two 
triangles adjacent to one diagonal form a cyclic 
quadrilateral. The two angles that are opposite to 
the common diagonal are supplementary, and there- 
fore the total conductance corresponding to the 
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Fic. 5. Square-section cable in square sheath. 
No nodes between boundaries. 


diagonal equals zero. Hence, a cyclic polygon in an 
isotropic medium is represented by a polygon of 
resistors along its sides without any resistor across 
the diagonals. This arrangement is known for 
squares and regular hexagons. 


IX. NOTES FOR THE PRACTICAL USE OF THE 
METHOD WITH SPECIAL REGARD TO 
SOUTHWELL’S RELAXATION METHOD 


The continuum between the boundaries is tri- 
angulated according to the viewpoints given in 
Section IV. If a boundary has sharp corners, nodes 
are placed at these corners. If the sharp corner is 
concave, the network is arranged in such a way 
that the total conductance of a resistor connecting 
the corner node with a node inside the medium is 
nearly zero. At convex sharp corners the conduc- 
tance of resistors leaving the corner between the 
two perpendiculars to the boundary should not 
differ very much from zero. 

As to the conductivity, a reference value equal 
to unity may be assumed and an appropriate multi- 
plying factor be applied in the final result. For non- 
homogeneous media the value of the conductivities 
at the center of gravity of each triangle considered 
or the arithmetic mean of the values of the con- 
ductivities at its vertices is substituted in Eq. (25). 
For computing partial conductances or for checking 
their values the relation 


1—cota;-cota2 
cota;= (27) 
cota,;+cotae 





is useful, where ai, a2, and a3 are the angles of a 
triangle. Problems in connection with three-dimen- 
sional fields of axial symmetry can be solved by 
considering the plane field of a section containing 
the axis and substituting in (25) for yz. and yy the 
product of the actual conductivities at the indi- 
vidual points and the distance of the points from 
the axis. If the conductivity at a point depends on 
the current density at that point, as is the case in 
the analog of magnetic flux problems in ferro- 
magnetic materials, successive approximations lead 
to the final solution. 


VOLUME 20, MAY, 1949 


Regarding transverse loading by current densi- 
ties two simple procedures are possible. One is to 
place at each vertex of a triangle a load equal to one 
third of the area of the triangle multiplied by the 
current density of the loading at its center of 
gravity and to sum up the loads at each node. The 
other is to sum up the areas of all triangles meeting 
at a certain node and to multiply one third of the 
total of the areas by the current density of the load- 
ing at the node. Cyclic or regular polygons are to 
be divided into triangles for this purpose. 

If Southwell’s relaxation method is used, one 
can work on a graph of the net, as is customary in 
this method. The values of the total conductances 


between the nodes and the sum of the total con- 


ductances meeting at a node are put down once and 
for all, generally in ink. The potentials of the nodes, 
the currents established in the conductances be- 
tween the nodes, and the residuals at the nodes are 
put down in pencil and amended during the process 
of the liquidation of residuals. 

As to the grading of the net, a special procedure 
which is a generalization of an artifice developed by 
Southwellf should be given and explained in terms 
of the conducting sheet analog. The continuum 
is divided into two parts by two sets of arbitrarily 
selected reference nodes which lie on two contours 
approximately parallel and reasonably close to each 
other. Figure 4 shows an example of the double 
set of nodes. They are marked by asterisks. The 
procedure is to split up the field problem into two 
separate field problems, one being the problem of 
field A between the boundary A and the set of 
reference nodes nearer to the boundary B; and the 
other being the problem of Field B between the 
boundary B and the set of reference nodes nearer 
to the boundary A. The nodes for field A are the 
nodes marked by crosses and asterisks and the 
nodes for field B are the nodes marked by circles and 
asterisks. The two problems are solved separately, 
the only link being that the potentials of the nodes 
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Fic. 6. Square-section cable in square sheath. 
Three nodes between boundaries. 


t See reference 1, p. 98 and Fig. 40. , 
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marked by asterisks must be the same in both 
partial problems. By this arrangement the require- 
ment is met that the values of the potentials and of 
the gradients in the two partial problems are ap- 
proximately equal near the reference nodes. The 
solutions of the two partial problems taken together 
are the solution of the main problem. 


X. EXAMPLE 


As an example, the capacity of a square-section 
cable in a square sheath with an homogeneous 
isotropic dielectric of dielectric constant unity will 
be computed. This example was treated by Gandy 
and Southwell.tf The ratio of the two sides of the 
square is 2-v2. Figure 5 shows the triangulation of 
one-eighth of the cross section without any nodes 
between the boundaries. In the graph double the 
values of the partial conductances and the arith- 
metic means of these values, i.e., the total conduc- 
tances are entered. If C is the required capacity, 
one obtains at once without any liquidation of resid- 
uals 4rC=4X0.547+8 X0.773 =8.372. The value 
given by Gandy and Southwell is 4rC=6.715, or 
the error is +24.7 percent. 

Figure 6 shows the treatment of the same prob- 
lem with three nodes between the boundaries, and 
the final solution found by the relaxation process is 
entered in the figure. The numbers under arcs are 
total conductances, the numbers in rectangles their 
sums, the numbers at the left-hand sides of nodes 
the potentials of the nodes, the numbers near ar- 
rows the currents established in the total conduc- 


tf See reference 1, p. 133, and R. W. C. Gandy and R. V. 
Southwell, Phil. Trans. Roy. Soc. A238, 453-75 (1940). 
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tances, and the numbers at the right-hand sides of 
nodes the residual currents. The capacity of the 
cable can be found either from the total current 
leaving the conductor or the total current entering 
the sheath. The first alternative leads to 


8 xX601+4X597 
4nrC =— =7.196, 
1000 





and the second to 


8X (228+417)+4x 506 
4rC= = 7.184. 
1000 





The deviation of the mean of 4%™C=7.19) from 
Gandy’s and Southwell’s value is +7.1 percent. 

The author computed the capacity for different 
arrangements and different numbers of nodes. For 
6 nodes and a triangulation that gave closer meshes 
near the vertex of the inner square than at other 
points of the field the error was +3.5 percent. For 
a square net without grading and 15 nodes the 
error was +2 percent. 


XI. CONCLUSION 


A procedure has been developed according to 
which any differential equation of the type of 
Eq. (1) can be solved approximately by a relaxa- 
tion method whereby the use of only five different 
formulas is necessary, viz., the use of formulas (2), 
(4)-(6), and (24) or (25), as the case may be. The 
shape of the boundary and the type of the boundary 
conditions are irrelevant as far as the procedure is 
concerned. 
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On the Velocity-Dependent Characteristics of High Frequency Tubes 
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Department of Physics and Institute for Atomic Research, Iowa State College, Ames, Iowa 


(Received July 9, 1948) 


The electronic transadmittances, transmodulation functions, and other characteristics of a planar, 
multigrid high frequency tube are expressed in terms of four fundamental admittance and modula- 
tion functions which are sufficient for the description of the low level r-f properties of the electron 
beam. A matrix representation is introduced based on the successive velocity moments of the beam 
which leads to the formulation of the characteristics as products involving excitation, grid trans- 
mission, modulation, and modulation admittance matrices. The electronic equations obtained are a 
formulation in which the Llewellyn approximation that the second and all higher velocity moments 
are zero is not made. Thus the theory allows for higher order velocity effects but is dependent on the 
knowledge of the fundamental admittance and modulation functions. 





I. INTRODUCTION 


HE pioneer work of Benham! and Llewellyn? 

on transit effects in electron beams opened 
the field of velocity-dependent characteristics, 
which characteristics are associated with all elec- 
tronic high frequency phenomena. Velocity modu- 
lation soon became a valuable technique at all 
power levels of operation.’ In the low level field 
(small modulation amplitudes) Llewellyn made the 
remarkable achievement of combining density and 
velocity modulation effects in a planar diode in a 
manner completely consistent with electron dy- 
namics and field divergence.‘ This is a single- 
velocity theory which does not include effects of 
velocity spread. 

The electronic characteristics of a planar grid 
tube can be expressed in terms of certain funda- 
mental admittance and modulation functions 
Vm(X,0), Vm'(x), w(x,v;0’), w'(x;v’), which are suffi- 
cient for the description of the r-f components of 
electron beam within a single region of the tube. 
Under conditions of low level operation, which are 
assumed in this paper, the r-f current density at a 
specified plane parallel to the grids resulting from 
electrons of velocity v is the sum of a current 
density proportional to the r-f voltage of the region 
and a current density which arises from modula- 
tions present in the beam as it enters the region. 
Thus, within the region a distance x from its first 
bounding plane, the r-f current density for elec- 
trons of velocity v is 


4 (X,V) = —Vm(x,0) Vo tim(x,2), (1) 


where V, is the r-f gap voltage, ym(x,v) is the 
modulation admittance function at this plane for 





1W. E. Benham, Phil. Mag. 5, 641 (1928); 11, 457 (1934). 
(1936) Llewellyn, Bell Sys. Tech. J. 14, 632 (1935); 15, 575 

?W. C. Hahn, G. E. Review 42, 258 (1939); W. C. Hahn 
and G, E. Metcalf, Proc. I.R.E. 27, 106 (1939). R. H. Varian 
and S. F. Varian, J. App. Phys. 10, 321 (1939). D. L. Webster, 
J. App. Phys. 10, 501 (1939); 10, 864 (1939). 

‘F. B. Llewellyn, Electron Inertia Effects (Cambridge 
University Press, London, 1939). 
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electrons of this velocity, and 7,,(x,v) is the modula- 
tion current density just mentioned.® If electrons 
enter only at the plane x=0, a condition which is 
also assumed throughout this paper, the modula- 
tion current density is an expression of the form 


eo} 


im(X,v) = dv’ u(x,v ;v’)t..(0,v’) 
MO u(x,v0'(0))i(0,0'(0))0.4(0), (2) 


where u(x,v;v’) is the density modulation function at 
the plane x for electrons of velocity v arising from 
modulations due to electrons injected with velocity 
v’; the d.c. current density at the injection plane is 
1(0,v’), and the quantities v'(0) and v,'(0) are the 
d.c. and r-f components, respectively, of the lowest 
velocity with which electrons enter the region (if 
the injection plane is the cathode, v'(0) and v,'(0) 
are zero). 

The first problem of a complete theory is the 
calculation of the admittance and modulation func- 
tions from the dynamics of the electrons and their 
contributions to the fields. Since the charge density 
at any plane depends on the electrons arriving at 
that plane, and hence on the electric fields through 
which they have passed since injection, the electric 
field appears on both sides of the divergence equa- 
tion relating charge density and field. This equa- 
tion is, in fact, an integral equation for the field, 
which in the small signal approximation breaks 
up into an integral equation for the determination 
of the d.c. field and an integral equation linear in 
r-f amplitudes for the determination of the r-f 
field. The admittance and modulation functions 
follow directly from a knowledge of these fields. It 
is assumed that this problem has been solved and 
that these functions are known.® 

In the absence of the r-f fields, the motion of the 
electrons is controlled by the emission characteris- 


5 See chapters 3 and 5 of D. R. Hamilton, J. K. Knipp, and 
J. B. H. Kuper, Klystrons and Microwave Triodes (McGraw- 
Hill Book Company, Inc., New York, 1948). 

6 For a discussion, see the last section of this paper. 
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tics of the cathode and the d.c. voltages on the 
grids and anode. In order to insure uniformity over 
planes parallel to the grids, a uniformly emitting 
cathode, shallow regions, and reasonably fine mesh 
grids are assumed. Because of d.c. field penetration 
of the grid, it is possible, even if the control grid 
is negative, that as far as the motion of the electrons 
is concerned the effective voltage of that grid is 
positive; but nevertheless, since it is negative, 
electron collection by that grid is practically zero. 
It is assumed, however, that there is no r-f pene- 
tration of the grids. The action of grid (k—1) on 
the r-f current density is described by 


i .(0,v) =a(v) “4-2, (d,v) 4 +4,(0) 4, (3) 


where 7,(0,v)“ is the injected r-f current density 
in the region (k), 7,(d,v)“*- is the r-f current den- 
sity at the exit plane of region (kR—1), a(v)“-” is 
the velocity-dependent transmission function of 
grid (k—1), separating the two regions, and 7,(v) “*- 
is the grid modulation current density introduced 
into the beam by that grid. The grid characteristics 
a(v) and i,(v) are determined by the geometry and 
nature of the fields in the neighborhood of the grid 
and by the details of the process of capture and 
the production of secondaries. 


Il. BEAM AND MODULATION CURRENTS 


The r-f electron current at the plane x within a 
particular region has the form 


t(x) = —Ym(x) Vist tm(x), (4) 


where ym(x) is the density modulation admittance 
at the plane x and 7,,(x) is the modulation current 
at that plane arising from modulations present in 
the beam at injection. These quantities can be 
expressed in terms of the admittance and modula- 
tion functions. 

It is first to be noted that Eq. (1) must be supple- 
mented by an expression for the r-f component of 
the lowest (or most negative) velocity at that 
plane. The required relation can be written in 
the form 


v(x) = =ya(a)( 


v'(x) 
~IJe 


in which v!(x) is the d.c. component of this velocity, 
and —TIp, is the d.c. beam current in the region; 
the first term is the modulation in the lowest ve- 
locity proportional to the gap voltage and the 
‘second term is the modulation arising from modula- 
tions present in the beam at injection. The latter 
has the form 





) Vattm'(x), (la) 





v'(x) oo 
Um! (x) = | dv’ p'(x ;v’)t..(0,v’) 
—Lol” »l(o) 
—p'(x ;0'(0))2(0,0'(0))v.'(0) ¢. (2a) 
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Equations (1) and (1a), together with Eqs. (2) and 
(2a), give a complete description of the r-f prop- 


erties of the beam at the plane x in terms of the 


gap voltage V, and 7,(0,v’) and v,'(0), which de- 
scribe the r-f properties at the injection plane. 
Substitution of Eqs. (1) and (1a) in 


= 


ty(x) = f dvi .(x,v) —1(x,v'(x))v.'(x), (5) 
vl (zx) 


and comparison with Eq. (4) yields for the density 
modulation admittance and the modulation current, 
respectively, 


Vm(x) = dvYm(X,v) 
vl(z) ' ; 
—i(x,v'(x))¥m!'(x)(v'(x)/ —Lo), 
im(x) = dvi m(x,v) —1(x,0'(x) vm! (x). (6) 
vl (x) 


Substitution of Eqs. (2) and (2a) in the second of 
Eqs. (6) gives 


3) 


dv’ u(x 3v’)t.(0,v’) 
v!(0) ’ fas 
— (x ;w'(0))7(0,0'(0))v4'(0), (7) 


with the modulation function, 


tn (X) = 


( ) 
p(x ;v’) = dvp(x,v ;v’) 
vl (z) 


—i(x,v'(x))(v'(x)/—Lo)u"(x 50’). (8) 
III. DRIVING AND EXCITING CURRENTS 


The circuits of a high frequency tube are usually 
cavity resonators with the gaps between the cath- 
ode, grids, and anode forming the capacitances. 
Since the regions are shallow, the fields in a par- 
ticular gap in the absence of the beam are nearly 
space-constant and the beam is coupled to the 
resonator and load through the space average of 
the electron current in the gap. This average is 
called the driving current i..; thus 


ie=t/af dxi (x). (9) 


If Ye and Y, are the resonator and load admit- 
tances, respectively, measured at the gap, 74=(Yr 
+ Yz)V.. The space average of Eq. (4) gives 


a) aes (10) 


where 
Ya=1/af dxYm(x), in= t/a f dxim(x), (11) 
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are the beam-loading admittance and exciting cur- 

rent, respectively. If the total circuit admittance is 

Y=YetYrt+ V1, the gap voltage arising from 

modulations in the beam at injection is given by 

the simple relation 7,,= YV.. The exciting current 

is given by 

in= f dv’ u(v’)i..(0,v’) 

v (0) , 

—u(v'(0))2(0,0'(0))v..'(0), 


where u(v’) is the excitation function for modula- 
tions in the injected beam due to electrons having 
velocities v’ : 


(12) 


d 
p(v’) =1 af dxp(x sv’). (13) 
0 


IV. EXCITING CURRENT IN THE CATHODE REGION 


Since the electrons emitted by the cathode have 
all velocities from 0 to «, the exciting current in 
the cathode region is simply 

4) 

inom f dv’ p(v’) ‘2 ..(0,v’). (14) 

0 
The cathode excitation function y(v’)© is the 
space average of the current modulation function 
u(x 37’), which as in Eq. (8) can be expressed in 
terms of the modulation functions u(x,v;v’) and 
u(x 30’), defined by Eqs. (2) and (2a). In the 
cathode region electrons may be reflected from a 
potential minimum as a result of space charge or 
from a retarding field arising from a negative ef- 
fective voltage on the control grid. Hence at a 
given plane in the region, the velocity distribution 
may extend down to a lower negative value; such 
is implied in Eq. (8). The definitions of Eqs. (2) 
and (2a) allow for this possibility. 

The beam-loading admittance of the cathode 
region Ypg) is the space average of ym(x)‘ as 
given by the first of Eqs. (6). The modulation ad- 
mittances Ym(x,v) and ym'(x) are used without 
modification as defined in Eqs. (1) and (1a). 


V. TRANSFER RELATIONS FROM REGIONS 
(k—1) TO (hk) 


Separating regions (k) and (k—1) is grid (k—1). 
The action of the grid is described in the Introduc- 
tion, to which description is added the assumption 
that the r-f components of the velocity of the slow- 
est electrons approaching the grid and the slowest 
electrons leaving it after transmission are identical. 
This limit does not apply to the grid modulation 
current density i,(v)“-», which has the lower 
limit 0. Hence 


14(0,v’)\ = a(v’) D4 (dv) &-Y +i,(v’)-», 


ve'(0) =v,4'(d)-», (15) 
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These relations can be used in Eq. (12) for the 
exciting current, 


tm” = f dv’ u(v’)™a(v’) 94, (d,v’) &-Y 
vl (0) (k) 
— p(v'(0))2(0,0'(0)) av '(d) FY 


+f dv'p(v’)™i,(v’)*-, (16) 


On using Eqs. (1) and (1a) for region (k—1), the 
exciting current of region (%) is transformed into 
the sum of a current proportional to the gap voltage 


_of region. (k—1), an exciting current dependent on 


the exciting current of region (k—1), and an excit- 


ing current arising from grid (k—1). The result has 
the form 


in ® = — VY, HED VY AD4 5, D 4G RD (17) 
in which 


x 


Y,,(4 kD -{ dv’ u(v")ex(v’) Dyn (d,v") 4D 
vl (0) 


— u(0"(0))4(0,04(0))®yqt(d) 
‘ (v'(d) / — Io)», 


L 
tm'* k—1) — dv’ p(v’) a(v’) &-Y4,, (d,v’) FY 
“ vl(o) 


— w(v'(0))4(0,0'(0))om!(d)@—, (18) 


2 
1\* k—1) -{ dv'u(v’)“4,(v’)*-», 
0 


are, respectively, the electronic transadmittance 

from region (k—1) to region (k), the transexciting 

current from region (k—1) in region (k), and the 

grid exciting current from grid (k—1) in region (k). 
The use of Eqs. (2) and (2’) yields 


3) 


Lin * k—1) — f dv’’ u(v"") k—1)4 ,(0,0'’) (-D 
vl (oy (k -1) 


— p(v'(0))* *-Y4(0,0'(0))*—v4(0)4-Y, (19) 
in which 
u(v’’) oe ft dv’ u(v’)a(v’) 
vl (0) (k) 
X w(d,v’ jv’’)*-Y — p(v'(0))“2(0,0'(0)) 
x (0(d)/— To) utd p")* —(20) 


is the transmodulation function from region (k—1) 
to region (k). 
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VI. BEAM CURRENT WITHOUT REFLECTIONS 


If the beam at a given plane contains no elec- 
trons with negative velocities, it is convenient to 
describe it by giving the beam current, instantane- 
ous average velocity, and instantaneous velocity 
moments. Thus if i(x,v,t) is the current density and 
if v'(x,t) is the lower limit of the velocity distribu- 
tion, the beam current i(x,t) and instantaneous aver- 
age velocity v(x,t) are given by 


1(x,t) -{ dvi(x,v,t), 
v!(z,t) 


v(x,t)4(x,t) = f dvvi(x,v,t), (21) 
ol (z, t) 


and the mth instantaneous velocity moment A‘ (x,t) 
by 


eo) 


am(eited = f dv[v(x,t) —v ]"i(x,v,t). (22) 


v!(z,t) 


One notes that the zeroth moment is unity and the 
first moment is zero. 

For the purpose of calculating average values 
the current density can be represented by 


4(x,v,t) =i(x,t)8(x,v,t), (23) 
where 


@o 


1 
5(x,v,4) = >> A” (x,t)6™(v—v(x,t)). (24) 
o n! 


In this expression 6‘")(v—v(x,t)) is the mth deriva- 
tive of the Dirac delta-function. On substituting 
this expression in Eq. (22) it is seen at once that 
all the moments are given correctly. 

It follows that the representation of the d.c. and 
r-f components of the current density are, re- 
spectively, 


i(x,v) = —Iod(x,v), 
4 o(X,v) =74(x) 6(x,v) +I ov.(x) 5’ (x,v) 


o 1 


-Ihh> —A ai (x) 5 (x —0(x)), (25) 
2 Nn. 


where 


1 
5(x,v) =D —A~ (x)d(v—v(x)), (26) 
0 n! 


and 4’(x,v) is the first velocity derivative of 6(x,v). 
In these expressions A‘)(x) and A,‘(x) are the 
d.c. and r-f components of the mth velocity mo- 
ment, respectively. Since Eqs. (25) contain the 
velocity only through the delta-function and its 
derivatives, modulations of the lower velocity limit 
are to be neglected in using this representation for 
the d.c. and r-f current densities. 
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At this point it proves to be very convenient to 
introduce a matrix representation for 7.(x,v) which 
effects a separation of velocity dependence and r-f 
moments. This is done by writing 7,(x,v) as the 
product of two matrices, thus, 


to(x,v) = 8(x,v)i«(x), (27) 
where 
5(x,v) = | 5(x,v), —v(x)8’(x,v), 3A (x) 5 (v—v(x)), 
$4 (x)5@(v—v(x)), ---| (28) 
is a dimensionless one-row matrix, and 
i.(x) - —Io ta(x)/—Io 
0 o(x)/v(x) 
A. (x)/A® (x) (29) 


A(x) /A (x) 








is a one-column current operand. 
The current operand can be written in the form 


i.(x)= —Ym(x) Vatin(x), (30) 


where ym(x) is the modulation admittance matrix 
and i,(x) is the modulation current operand arising 
from the modulations present in the beam at in- 
jection. These two matrices are readily expressed 
in terms Of Ym(x,v), im(x,v), Ym'(x), Um'(x) of Eqs. 
(1) and (1a). In fact 


3) 


Ym(x) = dvym(x,v) 1 (x,v) 


vl(z) 


—4(x,0"(x))¥m'(x) (v'(x)/ — Lo) 1 (x,0"(x)), 


im(x) = Vim (x,v) 1(x,v) 
vl(z) 


—14(x,v'(x))0m'(x)1(x,0'(x)), (31) 
where 


1(x,v) =/1 
— (v(x) —v)/v(x) 
(v(x) —v)?—A® (x) ]/A® (x) 
[(v(x) —v)*—A® (x) 
—3A®(x)(v(x)—v)J/A@(x)| (32) 
[(o(x) —v)*—A™ (x) 
—4A ® (x) (v(x) —v) J/A“ (x) 








is a dimensionless one-column matrix. It has the 
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following properties : 


3(x,v)1(x,v’) = 5(v—v’) , farce 9)6(e9) =1, (33) 


where 1 is the identity matrix. 
The modulation current operand can be written 


in terms of the properties of the beam at the 
entrance plane by using Eqs. (2) and (2a). Thus 


i] 


in(x) = dv'u(x ;v’)i..(0,v’) 

"© 40,0") )04'(Ouew"(0)), (34) 
where 
u(x ;v’) = : dop(x,v ;v’)1(x,v) 


vl(z) 





-itex'ey (wes nese), 65) 


is the modulation matrix function. 


VII. MODULATION AND EXCITATION MATRICES 


It was assumed in the last section that the beam 
at the plane x contained no electrons with negative 
velocities. The r-f current density was then repre- 
sented as the product of two matrices, Eq. (27). 
If there are no negative-velocity electrons at the 
entrance plane, the r-f current density at that 
plane can be given the same matrix representation. 
When introduced into Eq. (33), there results 


im(x) =u(x)i.(0), (36) 
where 


u(x) = f do'y(x 0')8(0,0’) (37) 


is the doubly infinite (square) modulation matrix. 
Similarly, from Eq. (12) for the exciting current, 


tm=vi.(0), (38) 
where 


‘ f dv’ u(v’)8(0,0’) (39) 


is the one-row excitation matrix. 


VIII. TRANSMISSION AND GRID- 
CURRENT MATRICES 


If region (k) has no electrons with negative ve- 
locities, i,(0)‘* =7,,(0). Hence on combining 
Eqs. (34) and (15) and introducing the matrix 
representation of i,(d,v’)“*-) and observing that 
u(0;v’) =1(0,v’), there results 


i.(0) =a Vie +ie-, (40) 
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in which 


al) = f wv) Paiey8a0) (41) 
is the transmission matrix for grid (k—1) and 
ie) = f dv’1(0,0')i,(v')"- (42) 
0 


is the grid-current operand. 

If the transmission function a(v’)“-” is a con- 
stant, 1(0, v’)) =1(d,v’)“-), and the transmission 
matrix is a constant, by the second of Eqs. (33). 


IX. TRANSFER RELATIONS FROM REGIONS 
(k—1) TO (k) 

By Eq. (38), the exciting current for region (k) 
is tm =p i,(0). With the aid of Eq. (40), 
im =y[ ai ,(d)*-) +i,» ], Hence, by us- 
ing Eq. (30), 


in ®) = —ya(tDy,, (d) 4D VD 
+p aj, (d)AY + pi), (43) 


This matrix relation is of the same form as Eq. 
(17). The electronic transadmittance, transexciting 
current, and grid exciting current are, respectively, 


Vn FD = y® gt-Dy,, (gd) @-D, 
in FD = y) qj, (d)A-D, 
i 40h FD = yi e-D, (44) 
That Eqs. (44) are completely equivalent to 
Eqs. (18) can be shown by the substitution into 
the former Eqs. (39), (41), (31), (39), (42) and the 
use of the first of Eqs. (33). Since 7(0,v'(0))™ 
= a(v'(0)) *-Y2(0,0'(0))“-», Eqs. (18) follow as 
written. 


If region (k—1) has no negative-velocity elec- 
trons, the use of Eq. (36) gives 


dim (® FD) = yy, FD}, () FD, (45) 
Um 
(i, BV) = gy (#) gy kD yy (dq) H-D (46) 
Um 


where 


is the transmodulation matrix from region (k—1) to 
region (k). Eq. (45) also follows from Eqs. (19) 
and (20). 

In the cathode region it may prove necessary to 
use the more general Eq. (34), which yields for the 
transexciting current to the second region 


in = f dv’ tm(v’)*%i4(0,0'), (47) 
0 


where the transmodulation function is expressed by 
tin(0') 9 =yPay(d 0’). (48) 
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X. TRANSMODULATION MATRICES, TRANS- 
ADMITTANCES, TRANSMODULATION 
FUNCTIONS 


Repeated use of the relations of the last three 
sections leads to the following general equation for 
the r-f exciting current in region (k): 


k—1 
im'* = > [— Y,,* i) Vo + pm * Dj) 


7=1 


+m Pa i,(d). (49) 


In this expression pp“ ” is the transmodulation 
matrix from grid (j) to region (&). It is given by 


tn * ?) =p aD yy (d) FD gg un(d) (*-2) wie 


Q Vy (d) Grn, (50) 


The transadmittance Y,,‘“?) from region (j) to 
region (k) is simply 


Y,,* DP=p,* Pay, (d)™. (51) 

In the last term of Eq. (49), in(d)‘ can be sub- 
stituted as calculated from Eq. (34). The term can 
then be written in the form 


aD 


f dv’ pm (v’) 07 (0,v') © 


Io) () 
_ 1(0,v'(0))‘®v,,'(0) ( ) um (v'(0)) i) 


where up»(v’)“” is the transmodulation function 
from the entrance plane of region (2) to region (k). 
It is 


(52) 


bm(v’)  =y,,(* Pa n(d y’) 6 (53) 


The fully detailed exciting current of region () is 


k-1 
im (®) = > [- Y,,°* ® Vio +m Dig | 


j=! 


+f diiun(o)* (02). (54) 
0 


It contains currents due to the (k—1) voltages and 
(k—1) grids through which the beam has passed 
and a current due to cathode fluctuations. If there 
are no electrons in the tube with negative velocities, 
the latter current can also be expressed as a matrix 
produce in the form un“ i,(0)°, where un“ is 
the transmodulation matrix from the cathode plane 
to region (k). 


XI. EVALUATION 


In the Llewellyn approximation, which assumes 
that the second and all higher velocity moments 
are zero, the theory takes the form of three equa- 
tions, one for the driving current 7, (which is, in 
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fact, Eq. (10)), an equation for the density modula- 
tion 7,(x) at the plane x (which is Eq. (4)), and an 
equation giving the velocity modulation v,(x) at 
that plane: 


Ly = —T YaVuttm, 
1 w(x) = Tn (x) Ve tim(x), 
Va(X) 


on = 


v(x) 


u(x) 


These are the Llewellyn electronic equations. The 
coefficients ym(x) and f,(x) are the density and 
velocity modulation admittances, respectively, and 
Um(x) is the velocity modulation at the plane x due 
to density and velocity modulations present in 
the beam at injection. 

The second and third electronic equation can be 
combined as in Eq. (30) by setting 





1 »(X) 
i4(x) = | Vu(x) |, 
ee 
Ym(X) 
aad | F(x) 
Im(X) 
in (x) = Um (x) |, (56) 
| "va | 


thus giving them a matrix representation. 


The Llewellyn theory provides explicit expres- 
sions for im, tm(x), and Vm(x) : 


; ’ v4(0) 
in= wel 0) +m —Lo ), 
v(0) 





v0(0) 
dm(X) = wii(x)t.(0) +uits)( Sid —), 








v(0) 

m w(0 
ie ) = Hose) n(0) + Hola) ( —Te \ ’), (57) 

v(x) v(0) 


The exciting current 7, can be written as a matrix 
product as in Eq. (38) by using the excitation matrix 


u=([ui wv], (58) 


Likewise, the second and third equations can be 
represented as in Eq. (36) by means of the modula- 


tion matrix 
u(x) _ en canes 


My i(x) Muv(x) a 


Thus the Llewellyn electronic equations have a 
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matrix representation of the form of Eqs. (60), 
below. 

The Liewellyn calculations give the modulation 
admittance, excitation, and modulation matrices, 
as well as the beam-loading admittance, in terms 
of d.c. field and space-charge parameters and the 
d.c. transit angle.? The results have been used to 
great advantage by Llewellyn and Peterson in the 
discussion of vacuum-tube networks.* The approxi- 
mation has been used by Peterson in the discussion 
of noise suppression in microwave tetrodes;*® the 
contributions to the output noise density and ve- 
modulation currents arising from cathode fluctua- 
tions can be made to cancel each other by the 
proper choice of the grid-screen space-charge 
smoothing parameter. 

The Llewellyn approximation begins to lose its 
validity when the spread in velocities, whether 
thermal or due to secondaries, becomes comparable 
to the average velocity, and also at such frequencies 
that the spread leads to a spread in transit angles of 
an appreciable fraction of a radian. In the cathode 
region the theory is at best a crude approximation. 
However, under most operating conditions, its ac- 
curacy is quite sufficient in other regions. 

Unfortunately the Llewellyn theory has no 
simple extension. In a region having no reflected 
electrons, an extension would require the. calcula- 
tion of the elements of the electronic equations 


y= = Yp Votui.(0), 
i.(d) = —yn(d) V..+u(d)i.(0), (60) 


without assuming the second and all higher velocity 
moments to be zero (the Llewellyn approximation). 
This problem has been reduced to the solution of 
an integral equation the kernel of which is an 
integral over the d.c. velocity distribution at the 
entrance plane of the energy derivative of a d.c. 
phase factor.” In the Llewellyn approximation the 
solution of this equation yields a resolving kernel 
which leads to the Llewellyn equations. The general 
solution of the integral equation would lead to 
Eqs. (60). 

If the cathode region contains reflected electrons, 
the theory requires the calculation of the elements 


7D. R. Hamilton, J. K. Knipp, and J. B. H. Kuper, Klys- 
trons and Microwave Triodes (McGraw-Hill Book Company, 
Inc., New York, 1948). 
(1944) B. Llewellyn and L. C. Peterson, Proc. I.R.E. 32, 144 
*L. C. Peterson, Proc. I.R.E. 35, 1264 (1947). 
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of the somewhat more general electronic equations 


to=— YaVet f dv’ u(v')i.(0,v’), 
0 


iu(d) = —yn(d) Vet f dv'y(dw')ie(0,v’). (61) 


The reduction of this problem to an integral equa- 
tion will be the subject of an investigation by the 
author when he returns to the study of these ques- 
tions in the not too distant future. It is hoped that 
the formalism of this paper will facilitate study in 
this field, both experimental and theoretical. Many 


extensions suggest themselves, in particular, nega- 


tive-velocity secondary currents, double injection 
and two-way transmission, less simple geometries, 
non-linear trajectories, each of which presents its 
own difficulties. 

The Eqs. (61) are a shunt formulation in which 
the driving current and the current operand are 
expressed in terms of the r-f voltage, the exciting 
current, and the modulation current operand. A 
series formulation is also possible. Let C be the 
capacitance of the empty gap, J.=i.—jwCV. the 
total current, and Y,=jwC+Y, the gap admit- 
tance. Then 


Vo= - 216+ f dv'z,(v’)i.(0,v’), 
0 


i,(d) = (dot f dv'a(d;v')i,(0,v’), (62) 


is a formulation in which the r-f voltage and the 
current operand are expressed in terms of the total 
current, an exciting voltage, and a series current 
operand. Equations (61) and (62) are equivalent, 


for 
Z9(v’) =p(v’)Z,, 


a(d;v’) 
= —Yn(d)z,(v’)+u(d;v’). (63) 


It is perhaps well to point out that a matrix 
representation can be introduced for the beam 
even where there are electrons with negative ve- 
locities if care is taken to distinguish between posi- 
tive and negative currents. Thus Eq. (60) can be 
used in a cathode region containing reflected elec- 
trons if only entering electrons be used in calculat- 
ing the current operand i,(0). 


Z,Y,=1, 
a(d) = —Ym(d)Z,, 
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Barkhausen Noise and Magnetic Amplifiers. I. Theory of Magnetic Amplifiers* 


J. A. KruMHANSL** AND R. T. BEYER 
Brown University, Providence, Rhode Island 


(Received July 27, 1948) 


An analysis is made of the operation of magnetic amplifiers or non-linear transformers. Open circuit 
output voltages are calculated for transformers driven by a sinusoidal primary current, and also for 
sinusoidal primary current plus a direct current bias. The case of a loaded secondary is also considered 
and it is shown that there exist conditions for instability or resonance. A subsequent paper will consider 
Barkhausen noise as a limiting condition on the sensitivity of these devices. 





I. INTRODUCTION 


HE term magnetic amplifier or non-linear 
transformer is customarily applied to a trans- 
former containing ferromagnetic core material 
which, under the operating conditions, undergoes 
such magnetization that the flux density, B,' is not 
a linear function of the magnetic field intensity H. 
In the course of some experimental work previously 
reported,’ it was noticed that the background noise 
which ultimately limited the sensitivity of the de- 
vices could not be entirely accounted for by the as- 
sociated circuits. It seemed that this noise might be 
due, at least in part, to Barkhausen effect in the core 
materials. 

In order to discuss this subject this paper (1) 
presents a theory of operation of these magnetic 
amplifiers.* A subsequent paper (11) will present an 
analysis of the Barkhausen noise in terms of a 
simplified physical model. 


Il. METHOD AND ASSUMPTIONS 


The usual mode of operation of the magnetic 
amplifiers can be understood with the aid of Fig. 1. 
Two identical ferromagnetic transformers have their 
primary windings (N,) connected in series. The 
transformers are excited by an audio oscillator which 
drives the cores into the region of saturation during 
each cycle, thus generating a potential across each 
secondary winding which contains the fundamental 
frequency plus odd harmonics. If a given secondary 
winding (2) on each transformer is connected in 
series opposing (relative to the primary), these 
signals will cancel out. However, if a direct current 
is entered on one of the secondaries (N3), even 
harmonics will be produced across each secondary. 


* Part of the research reported in this paper was performed 
under contract NOrd 508. The authors wish to express their 
Seems to Professor T. B. Brown for his encouragement 
of the work done under this contract. 

** Now at Cornell University, Ithaca, New York. 

1 A table of chief symbols used in this paper is included in 
Appendix B. 
asap” Beyer, Miller, and Trischka, Proc. I.R.E. 35, 1375 

* The authors are indebted to the Bell Telephone Labora- 
tories, whose unpublished work leading to results similar to 
those of Part I was made available in the course of the work 
mentioned in reference 2. 
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Because of the method of connection, the contribu- 
tion from each transformer will add, so that the 
voltage across the No. 2 windings will contain only 
even harmonics. , 

The analysis in this paper will be based on one 
such transformer. In the analysis it is assumed that 
the relation B = F(#Z) is known, that there is no flux 
leakage, and that B is uniform over any cross 
section of the core. In most of the discussion, the 
effect of hysteresis is neglected. Finally, only the 
steady state is considered. 


III. EXPRESSIONS FOR INDUCED VOLTAGES 


If the constant cross-sectional area of the core is 
A, the magnetic flux through this area is BA. 
Around this core there may be a number of windings 
(W). Let the jth of these have N; turns. If the cur- 
rent in the jth winding is J;, the total mmf, M, is 


Ww 
M= 7 0.49N,I ;, (1) 


j=l 


where the summation is taken over all windings. 
Also, H= M/l, where / is the length of the core and 
H is the magnetic field intensity. 

The voltage appearing across the terminals of the 
jth coil is 


V;=N;X10-*d¢/dt, (2) 
so 
dB d {M 
v= AN;x10-+(—) — —). (3) 
, dH/ yomjdt\ | 


If in each mesh there is a generator with e.m.f. E; 
and linear impedance Z;(p) (where p is the operator 
d/dt), then the jth equation may be written 


: Ww 
E;=Z(p)I;+ Nigh, Ts, ee Iw) ox Nibli; (4) 





k=1 
where 
0.49rA dB 
g(11, Iz, +++, Iw) = x10-(—) 
! dH/ yoy. 


All voltages and currents may be expressed as 
harmonics of some basic frequency w. Accordingly 
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one can express any of these in terms of a Fourier 
series. Substitution in Eq. (4) and collection of 
coefficients for each harmonic would yield 2W 
simultaneous non-linear algebraic equations for 
these Fourier coefficients. Because of the tedious 
nature of such an operation, the discussion is limited 
to certain special cases, for which the equations can 
be “‘linearized.’” 

Under the assumption that J,>J;+1, the function 
g(1, I2, ---, Jw) can be expanded about the point P: 
I;, 0, 0, ---. That is, 


Ww Og 
gh, In, +++, Iw) =g(h)e+ DL (—) Ij+-+-. (5) 
i=2 \OI;/ p 


In the above expression, g(J,)p and each coefficient 
in the sum are functions of time. 


Then 
- 
gli, Io, +++, Tw) a Niple=elieNiphi 
k=1 
Ww og Ww 
+(Miph) > (—) TKe+g(hi)p & Nipli 
k=2 \OL,/ p k=2 


+higher order terms. (6) 


The first step in approximating solutions to Eq. 
(4) is to solve 


E;=Z;(p)Ij;+NiN gli) pph. (7) 


In this case, one can prescribe either the £; or the J; 
and solve for the other. Since Eq. (7) gives E; ex- 
plicitly, but J; only implicitly, the J; are prescribed. 


IV. FIRST APPROXIMATIONS FOR 
INDUCED VOLTAGES 


It shall be assumed throughout this analysis that 
I,=I 9 sinwt.6 Also, the analysis is carried out for 
one possible approximation to the average mag- 
netization characteristic of the form 


2S Tio 25 
B=— tan-'( =) == tana), (8) 


us T 


S is the saturation value of B and yp is the initial 
permeability. 





‘It is often possible to carry the analysis from this point 
through a Fourier calculation (such as that done in the next 
section) without the linearizing approximation. This is done 
in the work cited in reference 3. 

* Although the present discussion is concerned principally 
with sinusoidal excitation current, this same approximation is 
equally useful for other types of operating conditions. For 
many B—H characteristics, (dg/0I;)p is a sharply peaked 
function of J; in the neighborhood of certain values of J; and 
small elsewhere, which can simplify approximate evaluations 
of Fourier coefficients. 

* While this is not the usual experimental arrangement, it 
serves as a useful first approximation to what is otherwise a 
very involved analysis. 
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Fic. 1. Schematic diagram of windings on a pair of 
non-linear transformers. 


The non-linear term in Eq. (7) is now expanded in 
a Fourier series. After some calculation, it can be 
shown that the coefficients G,(cos) and H,(sin) are 
given by 








G, 2 7? cos  coswtdt 
== f (nus — 
H, T Y_7/2.sin 1+a?D? sin*wt 
where 
0.47A po 
yo———-X 10, 
0.494 
D 7 Nil. (10) 


The solution of integrals of this type is given in 
Appendix A. The results here are 


2y 1 
G,=— exp( —n sinh), n odd 
aD aD 


=0, m even or zero 


H,=0, all n, (11) 
and the final expression for the open circuit voltage 
appearing across the jth winding is 


wy 
e;= 2110NiN;— Zz 


aD n=1, 3,5 


1 
xexp( —n sinh") cosnwt. (12) 
aD 


It can easily be seen that as D->0 (the linear case), 
the higher harmonic voltages vanish. For 7=1, the 
coefficient of I19 coswt becomes the linear inductance 
of the winding. 

The relative magnitudes of the various harmonics 
of the open-circuit voltage appearing across the jth 
winding are plotted in Fig. 2. This graph indicates 
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Fic. 2. Relative amplitudes of odd harmonics produced in non- 


linear transformers. 


that there is an optimum driving condition for the 


production of any given harmonic. 


V. FURTHER APPROXIMATIONS FOR 


INDUCED VOLTAGES 


To consider the effects produced by currents other 
than J;, Eqs. (4) and (6) are combined to give 


E;=Z(p)Ij;+NjiNig(h) eplhitNigth)e 


Ww Ww Og 
x< 2. Nipli t+ Nil phi) > (—) Iyt+-+-. (13) 
O1,/ p 


k=2 k=2 


When the same form for the magnetization curve 
and primary current are used as in the preceding 


section, it can be shown that 








(14) 
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Fic. 3. Relative amplitudes of even harmonics produced in 


non-linear transformers. 


434 


Then Eq. (13) can be rewritten as 
E;=Zi(p)Lj;+NjNig(h) pplitNig(h)p 


Ww Og Ww 
x> wiph+ (pl) (=) > Nile. (15) 


k=2 j/ P k=2 
The various terms of this expression are then de- 


veloped in Fourier Series. If 


I;= $Piot+ y (P in cosnwt +Q jn sinnwt), (16) 
n=1 
and if the linear impedances at each harmonic are 
given by Rjn, X jn, then 


E;= $RyoP jo+ ; >» [(RinP jin tX jnQ jn) cosnwt 
n=1 


+(RjnQin—X jnP jn) sinnwt |+N;N; >» 


n=1, 3,5,++- 
2wy 1 
X—lo exp( —n sinh") cosnwt 
aD 


aD 


7 x 
+ NV — -—-—— + > 
(1+a?D?)} n=2, 4, 6, +++ 


1 
27 exp( —n sinh-—) 


a 





cosnwt 


(1+-a2D*)! 


Ww oo 
x| = > Ninw(— Pen sinnwt 


k=2 n=1 


+ Qin cosnat) | +N, - > 
n=?2, 4, 6, +--+ 


1 
2nwy exp( —n sinh-\—) 
aD 


(1+a2D*)! 


sinnwt 





W 
x| = > Nil Prot Pin cosnwt 


k=2 n=1 


+Qin sinmat) | (17) 


For one particularly simple special case, it is 
assumed that Ry=0, i.e., no d.c. in primary; 
I2=C=}P2 (d.c. on No. 2 winding, P2,=Q2,=0. 
n#0); W=2 (only two windings). In this case, 


substitution in Eq. (17) leads to an expression for 


the even harmonic voltage appearing across the No. 
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Fic. 4. Idealized B—H curve used in estimating effect 
of hysteresis. 


2 winding: 


1 
—2N2*ynw exp( —n sinh"—) 
aD 


Con = 





(1+a2D?)! 


Xsinnwt |C, n=2,4,6,---. (18) 


The amplitude of the e.m.f., éen, is thus directly 
proportional to the input direct current (for small 
C) as has been observed experimentally. The de- 
pendence of these voltages on aD is shown in Fig. 3. 

The effect of hysteresis can be studied by as- 
suming a hysteresis loop as shown in Fig. 4. In such 
a case, dB/dH is represented by a delta-function. If 
the analysis is carried out as above, the first ap- 
proximation indicates that the only effect of 
hysteresis is to introduce a phase shift in the output, 
without diminishing the amplitudes of the various 
harmonics. These amplitudes are affected by 
second-order terms, i.e., large direct currents. 

In the next important case it is supposed that J, 
is prescribed but that J, depends on the loading 
conditions for No. 2 winding. As before, J; = J,9sinwt, 
Ri=0, and W=2. Under these conditions Eq. (13) 
reduces to 


E;=Zj(p)I;+NjNigeplit NjNop(grl2). (19) 


It should be noted that gp is an even function of J; 
and pJ,; is an odd function of J;. In order to elimi- 
nate the effects of this large second term in Eq. (19), 
it is customary to connect the secondaries of two 
identical transformers in series opposing (relative to 
the primaries). The total voltage across the jth set 
of secondaries is then 


E;=Z(p)I;+2NjNop(grl2). 
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Fic. 5. Region of instability for tuned secondary on non- 
linear transformers. 


If the only appreciable currents in the No. 2 sec- 
ondary are the d.c. and the mth harmonic, then 


E2= RooP 20+ (QomR2m— P2mX om) sinmut 
+ (QomX 2m + PomReom) cosmwt + 2N22mw 


1 
— Pay exp( —2m sinh-\—) 4+] 
aD 





x 
(1+a2D*)! 
—_ 2P ay _ 
x0(- m sinh~ —) 

~ (1+e2D? (i+a2D)) 

. —Qomy 
X sinmwt +2N22mw| ——————_ 

(1+a?D%} 


x 





1 
exp( —2m sinh) —1 | cosmwt. (21) 
aD 


If the coefficients of sinmwt, cosmwt, and d.c. term 
are equated, and if the notation 


& 2N22mwy 1 
X m =———————_ ex( —2m sinh"), 


(1+a?D?)! aD 
- 2N2?*mwy 
Aw *—————_, 
(1+a*D*)! 
4N2.*mwy E2 1 
€20m -——~(=) exp( —m sinh-\— ) (22) 
(1+a?D?)!\ Roo aD 
is used, further calculation leads to 
20m = RomQ2m— (Xm+X 2m+X mo) Pan» (23) 


0 = (Xom— Xm+X mo) Qom+ R2P om 
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y their contributions cancel. Hence the value of the 
contour integral is the same as that of Eq. (25). 
Y + The residue at z= sinh~'1/8 is found to be 
1 
exp{| — XK sinh-- 
i ; B 
| Ro= ra ’ (27) 
2(1+87)! 
sinh tl 
| 6 while that at z=7 sinh—'1/8+-7 is 
-Ir or i x 1 
exp{ —XK sinh—'— 
8 a 
Ro=(-—1)* a ; (28) ’ 
° ° . 2(1+8") pe 
The value of the integral is then = 
in 
Fic. 6. Integration contour. 2 1 Sc 
———- ex( —K sinh"), m 
The current amplitudes are then (1+8*) B th 
Di by 
a a K even "% 
— A I(K, B) Gi 
. ™ =0 K odd. (29) ve 
(X om — Xm+Xmo)€20m tri 
Pom = — A ’ (24) For negative K, the integration is performed in the in) 
h lower half plane and Eq. (29) is again obtained, ex- qu 
where , , 
, cept that |K| replaces K. mi 
A = Rom? — Xn2+(X om t+ Xmo)?. ° | | , 
Since A vanishes for certain values of the parame- APPENDIX B co 
ters, it follows that there exists a region of insta- Table of Symbols per 
bility, the locus of which is the semicircle drawn in - 
Fig. 5. B Flux density. te! 
One should note that apparently the behavior in 7 Senge wued oo de 
“hae gg. M Magnetomotive force. of 
the vicinity of the mth harmonic is analogous toa _ 5 Saturation value of flux density. 
voltage generator which supplies the “open circuit” Flux. dis 
voltage €20m to a four-terminal network whose 4 Cross-sectional area of core. Cx 
transfer characteristics are described by Eq. (21). a Length of core. of 
W Number of windings on core. val 
N Number of turns on jth winding. 
APPENDIX A sg Current in jth winding. 7 
T . . . . . E; E.m.f. in jth mesh. 0 
lhe various Fourier expansions involve integrals Z.(p) dane tentiolionne tn Selb wneicls an 
of the form p d/dt. P 
, ‘ite Rin Resistance in jth mesh at mth harmonic. fin; 
. e'Kdz X jn Reactance of load impedance of jth winding 
I(K, B)= Ppa yey (25) at mth harmonic. of 
-r 1+? sin*s w Angular frequency of driving oscillator. loc 
, , Ho Initial permeability of core. cyl 
In the corresponding contour integral, the poles are a rpo/2S. lo 
given by C Direct current input in No. 2 winding. 2 
. Con Even harmonic voltages produced in No. 2 sou 
1 winding. =r 
z=+isinh—--+nr, n=0,+1,+4+2,---. (26) Ji Amplitude of current in No. 1 winding. : 
g(Ih, In, «++, Iw) y(1+e2D? sintwt). ~ 
g(Ii)p gp(Ii, 0, «+, 0). Con 
tetue I i sn Fi : ArA/l)woX10~. on 
The contour for positive K is shown in Fig. 6. Since ; . 40 /1)N r Con 
* ° ° a 218 4V11410- 1 
the integrand —0 as sx +4 «, the upper side of the Pin, Qin Peusler coatiddeate of enpeunian See Zo. (19 j 
contour does not contribute to the integral. Also, Gg, 4, Fourier coefficients of expansion for 2 « 
the opposite sides differ in value of z by 2x so that g(I10 sinwt) coswt. 7 
Vo! 
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On the Infinitely Long Cylindrical Antenna* 


Cwar_es H. Papas 
Cruft Laboratory, Harvard University Cambridge, Massachusetts 


(Received October 13, 1948) 


Using the method of steepest descents the far field and the asymptotic form of the current distribu- 
tion is obtained for an infinitely long, perfectly conducting cylindrical antenna excited by a localized 
electromotive force. The low frequency value of the radiation conductance is determined by inte- 
grating the radiated energy flux over a large sphere. 





I. INTRODUCTION 


HE problem of determining the longitudinal 

distribution of current and the input im- 
pedance of an infinitely long, perfectly conducting 
circular cylinder antenna has interested several 
investigators. Among them are Stratton and Chu,! 
Schelkunoff,? and Hallén.* And recently a new 
method of analyzing the problem was brought to 
the author’s attention through an unpublished note 
by H. Levine in which he elegantly formulates the 
input impedance by means of the free-space dyadic 
Green’s function and a variational principle in- 
volving the current. Using a simple plane wave- 
trial function, he obtains an expression for the 
input impedance which is valid for a broader fre- 
quency range than ‘the expressions of the afore- 
mentioned investigators. 

Unless one uses the variational approach, the 
core of the problem lies in the evaluation of a 
certain contour integral. Once this integral has been 
evaluated the current at any point along the an- 
tenna and hence its input impedance is readily 
determined. In this paper, by means of the method 
of steepest descents we evaluate this integral for 
distances far from the localized electromotive force 
exciting the antenna. This gives the explicit form 
of the far-zone field from which the asymptotic 
value of the antenna current is easily calculated. 
Moreover, the input conductance of the antenna is 
obtained by integrating the radiated energy flux 
over a large sphere. 

To formulate the problem we assume that an in- 
finitely long, perfectly conducting circular cylinder 
of radius a is driven at some cross section by a 
localized electromotive force, V, circumscribing the 
cylinder in a narrow peripheral band. Since only 
longitudinal currents can be generated by such a 
source, the axial and radial components of the mag- 





*The research reported in this paper was made possible 
through support extended Cruft Laboratory, Harvard Uni- 
versity, jointly by the Navy Department (ONR), the Signal 
Corps of the U. S. Army, and the U. S. Air Force, under ONR 
Contract NSori-76, T.O.I. 
sis) Stratton and L. J. Chu, J. App. Phys. 12, 230-235 
(1945). 

?S. A. Schelkunoff, Proc. I.R.E. 29, 493-521 (1941); Proc. 
I.R.E. 33, 872-878 (1945). 

*E. Hallén, J. App. Phys. 19, 1140-1147 (1948). 
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netic vector H and the circumferential component 
of the electric vector E are identically zero. That is, 


E=(E,, 0, E,) (1) 
H=(0, Hz, 0) (2) 


where E,, Hy, E, are the non-vanishing components 
of the field in cylindrical coordinates (Fig. 1). The 
exact expression for E,(p, ¢, 2) is given by 


V f Ho (*-$9)9] 





E.(p, dg, 2) ee 


edt (3) 
2a J _» Ho (k?—£?) 4a] 


where H,(x) is the Hankel function of the first 
kind and zero order, k=22/i, being the free- 
space wave-length, and ¢ is a complex variable. 
The path of integration is along the real axis of 
the complex ¢-plane with an upward indentatioa at 
¢=—k and a downward one at {=+ (Fig. 2). 
This equation can be obtained from previous in- . 
vestigations,/** or from a report by Papas‘ in 
which he derives integral expressions for the field 
components and the antenna current by means of a 
Green’s function used by Levine and Schwinger.® 
Our first task is to find the explicit form of Eq. (3) 
for large values of r. 


II. FAR FIELD 


To find the far field, i.e., find E, when r is large, 
the method of steepest descents® is used. That is, 
the integral 





E.(p, g, 2) eases 


Vi p® Ho [ (R29) 
f o[(k?—$?) Ps (3) 
Qn J_., Ho[(k2—¢2)4a] 


is evaluated for large r. Let u=(k?—{?)!; then for p 
large and “+0, 


2 )3 oa 
Hf (up)e**~| — exp] i( up—+3*) | (4) 


Tup 


4 C. H. Papas, “‘On the Infinitely Long Cylindrical Antenna,” 
Cruft Laboratory Technical Report No. 58, Harvard Uni- 
versity, September 10, 1948. 

° H. Levine and J. Schwinger, Phys. Rev. 73, 383 (1948). 

*P. Debye, Math. Ann. LXVII, 535-558 (1909); “Semi- 
konvergente Entwickelungen fur die Zylinderfunktionen und 
ihre Ausdehung ins Komplexe,’’ Munchener Sitzungsberichte, 
XL (1910). 
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Fic. 1. Coordinate systems used in the formulation of the 
problem. Cylindrical coordinates: p, ¢, z. Spherical coordinates: 
r, 0, o. 


The value of ¢ which makes the exponential 
euetts) an extremum, is a saddlepoint and is de- 
noted by fo. Hence fo is the solution of 


d 
—[(k?—¢?)'!p +o2]=0. 
d¢ 


Carrying out the differentiation, it is easily found 
that 
kz 
fo=———_—_ = k sin0 (5) 
(p?+-2")! 


where @ is measured from the plane z=0. At [=fo 


i(upt+{z) — pik(p cos6+z siné) 
é e 


and in the neighborhood of {> it is approximately 
equal to 


ip 
exp] ik cos@+z sin@) — ———-10] (6) 


cos*@ 


The path of integration through the saddlepoint is 
determined by the constancy of the imaginary part 
of the exponent of Eq. (6); that is, the path is so 
chosen that the imaginary part of the exponent 
remains constant as the path of integration traverses 
the saddlepoint. Since the exponent is a monogenic 
function and its imaginary part remains constant, 
its real part rises from negligibly small values to a 
peak value at the saddlepoint and then down 
again. The requirement for the constancy of the 
imaginary part of the exponent is expressed by 


I.P.[ik(p cosd+z sin@) ] 





ip(¢—$0)? 
=1.P| ik(p cosé+<2z sin@) —— | (7) 
2k cos*@ 


Let ¢{—f)=se**. In this, a is the angle which the 
path makes with the real axis of the ¢-plane. Sub- 
stituting this into Eq. (7), it is found that if (7) is to 
be satisfied, it is necessary that the R.P. (se?**) =0; 
i.e., cosla=0 or a=+7/4. The value a=7/4 is 
excluded because the path through {9 must extend 
from the second to the fourth quadrant of the com- 
plex ¢plane; this leaves a= —2/4 as the angle 
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which the path of steepest descent makes with the 
real axis of the complex ¢-plane. 

We now transform’ from the ¢-plane to the 
complex r-plane by the following transformation 


¢f=ksinr 


7T=60+12y. (8) 


The path of integration in the ¢-plane (Fig. 2) 
transforms into path C in the r-plane of Fig. 3. 
The saddlepoints lie on the real axis of the r-plane, 
i.e., T=8@ is a saddlepoint. We de-form C into Co, 
which starts at —(r/2)+i, goes through @=0, 
and then ends at 7/2+iz. A path such as C, 
passing through a saddlepoint @ would have to 
start at @—(2/2)+i and terminate at 0+(2/2) 
+i. In the r-plane, Eq. (4) transforms into 


2 b} 
(——} e~ i(r/4) pike cosr+z sinr) 
rpk cost] © 


2 j 
= (——) ei! 4) pikr cos(r—6) | (9) 
apk cost 


Let F(k cosr) represent the factor 
(Ho ((k?—§*) 4a }}7 


in the r-plane. Substituting this and Eq. (9) into 
(3), we get 


V 0+(4/2)—ixe 2 } 
E.(p, 2) =—e~i(n/4) f |__| 
T @—(x/2)+i0 lapk cost 


X ei*r cos(t—9) F(R cosr)k costdr. (10) 


This integration is carried out along the path C». 
It is only necessary to carry out the integration on 
a short segment of the path in the neighborhood of 
7r=06. We set r—0= ne~*‘*/® where 7 is the distance 
along the path measured from the saddlepoint. 
When |7—86| is small, 


. (r—6)? ne . 
cos(7 — 6) = 1 —————- = 1 —$ —e*"!?), 
2 
Moreover, dr=dyne~*‘*/®, With these, Eq. (10) 
takes the form 


Vi 
E.(p, 2) = ——e“' F(R cos@) 
T 


2k cos@)' ¢* 
[AEM F cindy 
Tp —e 





Since r is large we can replace the limits, —e to ¢, 


7F. Noether, ‘“Elektromagnetische Wellen an einem Draht, 


_ bei konzentrierter Energiequelle;’’ Physik Zeits. USSR, Band 


8, Heft 1, 1-24 (1935). 
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by infinite ones. That is, 
‘ - 2r\} 
f e-tntindy = f e-wvindy= (—*) . (12) 
—e — kr 
Recalling that r cos@= p, Eq. (11) takes the form 
Vi e*r 1 


E,(r, ¢) = —-— — . (13) 
x r Ho (ka cosé) 





This is the far field of an infinitely long antenna 
driven by a localized e.m.f. V. The 6-component of 
the E-field and the ¢-component of the H-field are 
given respectively by 











E.(r, @) Vi er 1 
E,(r, 0) = = —— (14) 
cos0 x rcosé Hy (ka cosé) 
and 
Ee(r, 0) Vi eter 
H,(r, 6) = =o 
1207 1207? r cos@ 





x . (15) 
H,™ (ka cosé) 


Ill. THE RADIATION CONDUCTANCE 


To find the radiation conductance of the antenna, 
the radial component of Poynting’s vector S must 
be integrated over the surface of a large sphere. 
This integration yields the radiated power P, and 
division by V?/2 gives the radiation conductance. 

Poynting’s vector is given by 


V? 1 


2407? r? cos?6 





1 
a = 


1 
x 
[Ho (ka cos@) | Ho (ka cosé) }* 





16) 


where the asterisk indicates conjugate complex 
values. Thus 





x/2 V2 x/2 1 
P= 2rr?S cosédé = f — 
—r/2 1207?/_x/2 cosé 
dé 





Xx (17) 
[ 79 (ka cos) |[ Ho (ka cosé) |* 


when ka is small, 
21 
H,™ (ka cos@)—1 —— log . 
T yka cos@ 


and the integral in Eq. (17) is approximately 
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Fic. 2. Paths of integration in the ¢-plane. 

















equal to 
w/2 dé x/2 
2 f =2 f 
0 21 2 2 0 
cos6| 1 —— log 
ra yka cos@ 
dé 
(18) 
2% 2 3 
cos6} 1 ——( log-— logka — log cost) 
Tv 7 








where log2/7=0.1159. This integral appears to be 
rather difficult to handle rigorously, and we there- 
fore resort to an approximation. 

Let us examine the integrand 


1 





(19) 


| 24 2 |? 
cos@| 1 — = (tog —logka —log cos# ) 
boy | 


Since ka is small, the term in square brackets can 
be replaced by (4/2*)(logka+log cos@)?. This quan- 
tity varies from (4/2*)(logka)? to (16/2?) (logka)? as 
6 varies from @=0 to @=6,=arc coska. For this 
range we replace (4/2*)(logka+log cos@)? by the 
geometric mean of its end values, i.e. (8/2?) (logka)*. 
Hence for the range @<@<6,, the contribution to 
Eq. (18) is 





1 7? 1 d@- x? 1 
af"s < s 
o 8 [logka]?cos@ 4 [logka]? 


ka 
3? log ctn(—) 
; 2 











T 6; 
X log tan(=+—) = , (20) 
4 2 4[logka |? 
1 \Co\% 
CG 
id 
Fic. 3. he eae 4M \ 
tion in the r-plane. “= 
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since 6; = arc coska = r/2 —ka. Moreover log ctnka/2 
+log2/ka=log2+logi/ka. Neglecting the term 
log2, (20) is approximately equal to 


1 x? 1 do 7 1 
2f cian, —ii 
o 8 [logka]*?cosé 4 1 


log— 
ka 





To complete the range of integration we again 
approximate the integrand Eq. (19) 


v/2 dé Tr: r/2 
2f —_——___+- | 
6 2 2 2 2 61 


cost = log —-—— 
x  ‘ykacosé 


x r 
(F-#)({(3-)]) 
2 2 
This integral can be immediately evaluated and has 
the value (#?/4)[1/log(1/ka) ]. By adding this re- 


sult to Eq. (21) the approximate value of (17) is 
obtained. Hence, 
1 g* 1 


V2 [x 
Re - pm eaten 
“ a. a. 


to 





log— log— 
a ka 


The radiation conductance G is given by 


2P 1 
G=—=—____.. (23) 


Vy? 1 
120 oe(—) 
ka 


On physical grounds G must be positive. Conse- 
quently it is necessary that 1/ka>1 or \>2za. 
(k=2x/d where d is the wave-length of the exciting 
source). This means that Eq. (23) is a low fre- 
quency, i.e., long wave-length, approximation to 
the radiation conductance. When ka = 1, G becomes 
infinite. 


IV. ASYMPTOTIC FORM OF THE CURRENT 


For great distances from the localized e.m.f. the 
carrent distribution can be computed from the far- 
field expression (15). For large values of r we found 
that 

Vi etkr 


1202? r cos@ Hy (ka cosé) 





H,(r, 0) =— (24) 


To obtain the value of H, on the surface of the 
cylinder for large z, it is necessary to restrict r and 6 
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so that rcos@=a; moreover, rsin@é=z. But since 
6=n/2, we have r=z. The argument of the Hankel 
function is given by the approximation: ka cos@ 
+ka*/z. Substituting these into Eq. (24) we obtain 
H, for p=a and 2z large: 


Vi et* 1 
(H,) = —- — ——, (25) 
p=a 1207? a ka? 
oe no() 


o 
~ 





The current for large z is then given by 
Vi eikz 


I(s) = 2xaH, = —— —— (0 


9 


607 ka? 
H0(—) 


ka? ka? 21 22 
—<i, H(—) = —— log . 
Z z ra yka? 





Since 





~ 


~ 


Using this approximation the asymptotic form of 
the current takes the form 








Vs eike 
I(z) =—— : (27) 
120 z 
log 
ka? 


Vv. CONCLUSIONS 


The far field (Eqs. (14), (15)) and the asymptotic 
form of the current at large distances from the 
source of excitation (Eq. (27)) were determined by 
the method of steepest descent. By integrating the 
radiated energy flux an approximate value of the 
radiation conductance was calculated (Eq. (23)). 
This integration, however, yields no information 
about the radiation susceptance since the sphere 
over which the integration is performed lies in the 
far zone where the time average energy flux is 
purely real. There is close agreement for long wave- 
lengths, i.e., for low frequencies, among the values 
of conductance obtained by Schelkunoff,? Hallén,* 
and Eq. (23). For shorter wave-lengths, Eq. (23) 
becomes invalid since the radiation conductance 
becomes infinite when \=2za and negative when 
A< 27a. 
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Note added in proof: Professor S. Silver has also 
derived (3). See Report No. 149, Antenna Lab., 
University of California, Jan. 3, 1949, ‘‘The Field of 
a Slot of Arbitrary Shape In an Infinite Cylinder.” 
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Phase Contrast in Electron Microscope Images* 


E. G. RAMBERG 
RCA Laboratories, Princeton, New Jersey 


(Received November 1, 1948) 


Image contrast in electron microscope images may arise from absorption, scattering, or phase 
changes impressed on the illuminating electron wave. In the present paper the intensity distribution 
in the in-focus and out-of-focus image of an edge of a transparent thin film introducing a prescribed 
phase change in the incident beam is calculated. It is found that the resulting “phase contrast” in- 
creases both with the film thickness (i.e., magnitude of phase change) and with the degree of de- 
focusing and, for thin specimens, exceeds other sources of contrast in magnitude. The model is too 
schematic to permit a quantitative comparison with measured intensity distributions. 





THIN specimen may modify an incident elec- 

tron wave by absorption, scattering, and 
changing the phase of the electron wave. All these 
may give rise to contrasts in electron microscope 
images. With thin specimens absorption is neg- 
ligible. Furthermore, the scattering of electrons 
beyond the effective aperture of the objective may 
yield only slight contrast. At the same time, phase 
changes in the electron waves traversing the speci- 
men may give rise to a pronounced delineation of 
specimen edges. This phase contrast, at the edge of 
a specimen of uniform thickness, will be investi- 
gated in the present paper. For reference, it may be 
recalled that a 100-angstrom-unit film of a material 
with an inner potential of 15 volts produces a phase 
delay of approximately 90° in a 50,000-volt electron 
beam.! 


I. INTENSITY DISTRIBUTION IN THE IN-FOCUS AND 
OUT-OF-FOCUS IMAGE OF AN OPAQUE EDGE 


The first step in determining the phase contrast 
in the image of the edge of a transparent film, 
whose sole effect on the incident electron wave is to 
introduce a prescribed phase change, is the calcu- 
lation of the wave amplitude distribution in the 
in-focus and out-of-focus image of a _ perfectly 
opaque edge. It will be assumed that a plane elec- 
tron wave is incident on an opaque half-plane a 
distance do from the front focal point of the objec- 
tive, which will be regarded as aberration-free 
within the limiting aperture of radius p=6f.? The 
geometrical arrangement is, thus, that pictured in 
Fig. 1. The notation employed is the following : 


f: focal length of objective. 
p= 6f: radius of lens aperture. 
dy: distance of front focal plane from edge. 


* Presented at the E. F. Burton Memorial Meeting of the 
Electron Microscope Society of America in Toronto on 
September 11, 1948. 

‘See, e.g., J. Hillier and E. G. Ramberg, J. App. Phys. 18, 
48-71 (1947). 

* This implies that the physical aperture of the objective is 
made appreciably smaller than the aperture which is optimum 
from the point of view of resolution. 
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d=f?/do: distance of image of edge from rear 
focal plane (d, do positive corresponds to too 
strong objective current). 

£, n, ¢: coordinates in aperture plane. 

x, y, 2: coordinates in reference plane. 

a=x/z 


Cw) +iStw)= f exp[_mix?/2 |dx ; 
0 


Siw) = f sinxdx/x. 
0 


Specifically, the amplitude distribution will be de- 
termined, for small values of a, along the x axis in a 
reference plane an infinite distance from the objec- 
tive (g= «), 

The procedure is to determine the wave ampli- 
tude u,° in the (image-side) aperture plane of the 
lens, assuming that the latter is perfect and un- 
restricted and to find the excitation at the reference 
point (x, y, ) corresponding to the wave amplitude 
distribution in the actual aperture (#-+7?< pp?) by 
means of the formula’ 


e plo —eyt 
up= (1/30) J f (z/r)up%E, 1) 


—(p?—¢)8 


0 
X—(e""/r)dédn, (1) 
of 











Fic. 1. Geometrical arrangement for the determination of 
the intensity distribution in the out-of-focus image of an 
opaque half-plane. 





’R. K. Luneberg, Mathematical Theory of Optics (Brown 
University, Providence, Rhode Island, 1944), p. 354. The 
author is indebted to Dr. M. Herzberger for pointing out 
Luneberg’s formulation of the diffraction problem. 
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with 


r=[(2—s)?+(x—£)?+(y— 0)? }}. 


The excitation in the aperture plane is simply the 
complex conjugate of that produced by a spherical 
wave converging at the focal point F, incident from 
the right on a real opaque half-plane placed at the 
geometric image of the opaque half-plane in object 
space. This is given by* 


up? = (i/2)(P+ 2+ n°)—! exp[ —ik( P+ P+ 1°)*] 


X {S(w) — C(w)+7[1 — C(w) — S(w)]}, (2) 
with 





k d ' , 
w=)— 
rftd(fP+P+n’)! 


ko d \'7) 3845 
a gas) (3) 
r fdtf) 4 f 








(Fy 40 








($) d,=2-1Gem 


4022 Fon $5310" 
4 
_/ 


——INTENSITY DISTRIBUTION IN 
IMAGE SPACE 


--~-INTENSITY DISTRIBUTION IN 
OBJECT SPACE 





($}4,: 5-10 em 





e 3 
d,+0.56 so For rd #310 


“0 “20 “10 ° 





(Has x10"tm 
20 i Av. 





10 20 
wf ron Ygs3-10" 


Fic. 2. The intensity distribution near the image of the edge 
of an opaque half-plane for different degrees of defocusing. 
4M. Born, Optik (Verlag Julius Springer, Berlin, 1933), p. 

191. 
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Equation (1) may be brought into reasonably simple 
form if following approximations are made: 


(i) 1/r and 1/f are neglected in comparison 
with k=27/) in the coefficients. 

(ii) (+ 9*)/(2f) is neglected in comparison with 
f in the coefficients. 

(iii) terms of order higher than (#-+-7?)/(2f) are 
neglected in comparison with f in the ex- 
ponents and arguments of the Fresnel in- 
tegrals. 

(iv) a is neglected in comparison with 1 in the 
coefficients. 


A check of numerical values indicates that all these 
approximations are justified in the present ap- 
plication. Then 


1 
up(x, 0, s) = ——— exp[ik(—f+211+ 7 ]}}) ] 
Anfz 


p (p? —¢2)8 ; 24+ 7? 
x2ik f ae f dn exp| itl at \| 
-» 0 af 


XLS(w) —C(w) +7{1—C(w) —S(w)j J, (3) 


tykf d \3 
He) 
fX\ x f+d 
Integration with respect to 7 simply replaces the 
factor exp[ —kn?/(2f) |dn by 


(2) ({e-0])-((4-0]})] 


which for large argument (£ not too close to p) is 
equivalent to (xf/(2k))e-‘*/4. For convenience in 
calculation it is thus possible to break up the in- 
‘tegral for up in the following manner, making use 
of the fact that C(—w) = —C(w), S(—w) = — S(w): 





with 


up —[4/(42) ](rkf)— exp[tk(—f+s2(1+7)') ] 
X {LitTet+TIs}. (4) 
Here 


I,=2k(1 -if dé exp —tkL —a&+ 2/(2f) J] 


~—i(4rkf)'{1+a(fk/x)! 
X[1 —cos(kp?/2f) —sin(ke*/2f)]} (5) 


JOURNAL OF APPLIED PHYSICS 








I; 


W 


I; 


WwW 


I; 


ti 


Sil 
tu 
sy 
at 
re 


fr 


In 


u} 





(5) 


ICS 





and 


— f " sin(kat) exp[—ik#/(2f)] 


x {[3-S(w) ]-114 — Cw) }dé 


w exp[ —ixw?/2 ] 


~—4nifa f 


0 
x (13 —C(w) +713 — Sw) ]} dw 
(olf (kil) sin[ (rkf)taw | 
4iak 
7 an (rkf)'aw 
afw afw 
aw 








—isin 


2d 2d 





> 4 | cos (6) 


with, e.g., w.=2. Finally, 
wy 
— i) dé exp[ik[ —at+ #/(2f)7] 
0 


x {[2—C(w’) J-i13-S@')]} (7) 
w’ =(k(p?— &)/(rf) }}. 


This has the approximate value 
1,;Si(4akf)' expl —tkp*/(2f) ] 
X {Jo(kap)+1Ho(kap)}. (8) 


Here Jo is Bessel’s function and Ho Struve’s func- 
tion of order zero. The last integral does not de- 
pend on the location of the edge and is clearly 
simply the Fresnel diffraction pattern of the aper- 
ture, vanishing except near the optic axis of the 
system.® It is thus not characteristic of the pattern 
at the edge of an object and is omitted from the 
remaining considerations. , 

Omitting a common factor —(1/z) exp[ik(—f 
+2(1+a?)!)] we thus find, for small deviations 
from focus and for the region close to the geometric 
image of the edge, the approximate relation 


‘ 


with 


m+ (1/x) f “(sinx/z) 


X {cos(x?/2a*dk) —7 sin(x?/2a*dk)}dx. (9) 


For exact focus (d—~«) this becomes 4+(1/7) X 
Si(pka). On the other hand, for large deviations 
from focus (i.e., d=f?/do small) the distribution 
approaches the Fresnel diffraction pattern at a 


_*See G. N. Watson, Theory of Bessel Functions (Cambridge 
University Press, London, 1922), p. 328. 

_ © In the complete absence of an edge the intensity distribu- 
tion near the center of the field is given, to the same order of 
approximation, by 1—[Jo(kap) ? exp[ —ikp?/(2f)]. 
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distance dy behind an edge illuminated by a parallel 
beam.’ 

Figure 2 shows the calculated intensity variation 
|up|? as function of af-p/f for an electron wave- 
length 5.10-'A (~60 kv) for perfect focus and for 
different distances dy between the edge and the 
focal plane; changing the sign of the defocusing 
merely changes the sign of the imaginary term in 
Eq. (9) and does not affect the intensity distribution 
in the pattern for the opaque edge. Furthermore, 
if the objective aperture is reduced from po/f to 
pi/f and the defocusing separation is, simultane- 
ously increased from do to d;=(po/pi)"do, the in- 
tensity pattern is simply expanded in the ratio 


45° PHASE SHIFT 








90° PHASE SHIFT 





3 
— G ds0 . 
aneee (F)-4,+2:10 om \ 
e? =o \ 
— (Fraes- 10 cm \ 
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Fic. 3. The intensity distribution near the image of the edge 


of a transparent half-plane imparting a prescribed phase delay 
on the incident electron wave. 





_ 7 Equation (9) is not an adequate approximation under these 
circumstances. 
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po/pi1 without changing relative intensities at 
corresponding points. The curves in Fig. 2 indicate 
the transition from the diffraction pattern arising 
from the restricted lens opening to the simple 
Fresnel pattern of an edge. At large distances from 
the edge the intensity oscillations are attenuated 
most. This is to be expected since here the maxima 
and minima lie closest together and hence are most 
affected by the limited resolution of the system. 


The intensity patterns expected for infinite lens” 


aperture are indicated for comparison. 


Il. PHASE CONTRAST AT THE EDGE OF A 
TRANSPARENT OBJECT 


The phase contrast at the image of the edge of a 
transparent object film is obtained by adding the 
amplitudes for two complementary edge diffraction 
patterns, one of which is shifted in phase by the 
phase delay 6 imparted by the transparent film: 


Intensity = | up(a)+e*up(—a) |”. (10) 


The phase delay 6 is here given by Ak-t=k-t 
*-AV/(2V), where AV is the inner potential of the 
film (~15 volts) and ¢ is the film thickness. As 
mentioned above, a film 100A in thickness is seen 
to result in a phase shift of the order of 2/2. 
Figure 3 shows the phase contrast for phase shift 
x/4, x/2, and w. The contrast is seen to increase 
with phase shift, reaching a maximum for z. For 
the last value, and generally at perfect forcus, the 
intensity distribution is symmetrical about the 
edge. Otherwise, changing the sign of the defocusing 
distance simply has the effect of changing the sign 
of a (i.e., the abscissas) on the intensity plot. 
One intensity distribution is the mirror image of 
the other through the edge. 

The edge contrast is seen to increase, also, with 
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the degree of defocusing. Even for perfect focus, 
however, the phase contrast at the edge may be con- 
siderably greater than the contrast arising from 
scattering beyond the effective objective aperture: 
For a 100A film of a substance with the effective 
atomic number 6 (carbon) and a specific gravity 1.5, 
the contrast arising from the scattering of 50-kv 
electrons beyond an aperture of 3.10-* radian may 
be estimated to be approximately 5:4,° while the 
phase contrast (for 90° phase shift) is seen to be of 
the order of 5:2. 

Caution is, of course, required in comparing the 
results here obtained with measured intensity dis- 
tributions at the images of film edges. Qualitatively, 
the shift of the first minimum from the shadow to 
the clear side of the edge in passing from too strong 
to too weak objective currents is, in fact, observed.°® 
Other features, such as the difference in contrast for 
the overfocused and underfocused image of the 
edge, which may be ascribed to the prismatic 
action of an edge of non-uniform thickness,!® are, 
of course, not reproduced. The same applies to the 
contrast arising from scattering and absorption 
losses. Finally, most high resolution images of 
film edges have been obtained with large’ physical 
objective aperture, so that spherical aberration 
influences the intensity distribution. Even so, the 
calculations give some indication of the importance 
which phase contrast may assume in revealing 
specimen contours in normal electron microscope 
operation. 


8See Zworykin, Morton, Ramberg, Hillier, and Vance, 
Electron Optics and the Electron Microscope (John Wiley and 
Sons, Inc., New York, 1945), p. 691. 

* See, e.g., J. Hillier and E. G. Ramberg, J. App. Phys. 18, 
48-71 (1947), Fig. 10. 

10 See H. Boersch, Physik. Zeits. 44, 32-38 (1943) and E. 
Kinder and A. Rocknagel, Optik 2, 346-363 (1947). 
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Flow through a Rocket Nozzle with and without Vibrational Equilibrium 


S. S. PENNER 
Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California 


(Received November 8, 1948) 


A procedure is outlined for calculating the maximum changes in the exhaust velocity of gases 
which can be introduced during flow through a rocket nozzle as the result of vibrational time lags, 
Representative calculations have been carried out for a number of propellant systems. These calcu- 
lations indicate that the vibrational time lags cannot produce changes in exhaust velocity which 
exceed a few percent. The numerical values of the changes are found to depend not only on the 
nature of the gases but also on the occurrence or non-occurrence of chemical reaction during flow. 





NOMENCLATURE 


c=exhaust velocity if complete equilibrium is reached 
with respect to internal vibration. 

c’=exhaust velocity if the vibrational heat capacity is 
frozen at the chamber temperature. 

C,=total heat capacity per mole at constant pressure 
including the heat capacity corresponding to dis- 
sociation. 

C,’=sum of molar rotational and translational heat 
capacities at constant pressure, including the heat 
capacity corresponding to dissociation. 

C,=average molar total heat capacity (including the 
heat capacity corresponding to dissociation) at 
constant pressure, between 7. and 7,. 

C,’=C,—C.i»=average molar rotational and transla- 
tional heat capacity (including the heat capacity 
corresponding to dissociation) at constant pressure, 
between T,’ and T.. 

Cyi» = vibrational heat capacity per mole. 

C.vi»=average molar vibrational heat capacity between 
T,. and T.. 

Cyi»’=average molar vibrational heat capacity between 
T.’ and T.. 

E=internal energy per mole. 

H=enthalpy per mole. 

AHT,.7¢=molar enthalpy change during flow through the 
rocket nozzle under conditions where complete 
equilibrium is reached with respect to internal 
vibration. 

AHT,’’«-=molar enthalpy change during flow through the 
rocket nozzle under conditions where the internal 
vibrations remain frozen at the chamber tempera- 
ture. 

p= pressure. 

p-=exit pressure. 

« = chamber pressure. 

R=gas constant per mole. 

T =absolute temperature. 

T.=chamber temperature. 

T.=exhaust temperature if complete equilibrium is 
reached with respect to internal vibration. 

T.’ =exhaust temperature if the vibrational heat capacity 
is frozen at the chamber temperature. 

T.i»= vibrational temperature defined through the actual 
energy residing in the vibrational degrees of freedom. 

u=linear flow velocity of the gases through a rocket 
nozzle, if the flow occurs parallel to the nozzle axis. 
v=volume. 


DIABATIC expansion of the combustion prod- 
ucts of a rocket motor through the nozzle! 
may occur so rapidly that there is not enough time 


1F. J. Malina, J. Franklin Inst. 230, 433 (1940). 





VOLUME 20, MAY, 1949 


to establish equilibrium with respect to the internal 
vibrations of the molecules. It is necessary to make 
a complete kinetic and flow analysis in order to de- 
termine the approach to flow with equilibrium 
partition of vibrational energy. The analysis of the 
flow process can be carried through, in a manner 
similar to that previously used to determine the 
presence or absence of equilibrium with respect to 
chemical reactions.” Some of the required kinetics 
data for calculations involving vibrational equi- 
librium may be found in the literature.*-® 

Before proceeding with the laborious analysis of 
the history of vibrational energy distribution during 
flow, it is desirable to determine the effect of various 
vibrational distributions on the performance of a 
rocket motor. Accordingly, the maximum perform- 
ance differences associated with the extremes of 
vibrational energy distribution which can be intro- 
duced during flow through the rocket nozzle have 
been estimated. Calculations have been carried out 
only for the case where complete equilibrium with 
respect to internal vibration is reached in the rocket 
chamber. Flow calculations with complete vibra- 
tional equilibrium and with internal vibrations 
frozen at the chamber temperature permit the 
determination of the maximum difference in exhaust 
velocity associated with the extremes of vibrational 
distribution introduced through zero or infinite 
vibrational time lag during expansion through the 
nozzle. The vibrational correction is assumed to be 
superimposed on the usual type of calculations for 
equilibrium flow or constant-composition flow.*? 

The method of calculation employed for the 
present study may be readily applied in those cases 
where the vibrational states are frozen at some tem- 
perature intermediate between the chamber and 
exhaust temperatures. 


2S. S. Penner and D. Altman, J. Franklin Inst. 245, 421 
(1948); D. Altman and S. S. Penner, J. Chem. Phys. 17, 56 
(1949). 

3W. Jost, Explosions- und Verbrennungsvorgdnge in Gasen 
(Verlag Julius Springer, Berlin, 1939), 

h., Lasies and E. Teller, Physik. Zeits. Sowjetunion 1, 34 
(1936). 


5 A. Kantrowitz, J. Chem. Phys. 14, 150 (1946). 


6 P. W. Huber and A. Kantrowitz, J. Chem. Phys. 15, 275 
(1947). 
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OUTLINE OF CALCULATIONS 


Thermodynamic considerations lead to two alter- 
nate methods of calculation, both of which are 
described below. 


A. Calculation of Exhaust Velocities from the 
Law of Conservative Energy 


Flow through the rocket nozzle is assumed to be 
adiabatic. Conservation of energy leads to the 
relations 


3c?=constant AHr,"*, (1) 

$c’*=constant AHr7,’". (2) 
From Eqs. (1) and (2) it follows that 

c/c’ = (AHr,7:/AHr7,'"*). (3) 


The exhaust temperatures 7, and T,’ are given by 
the relations!” 


(pe/ pe) BICr = T/T 0 
(pe/ pe) B/Cr’ = we | T- 


Equations (3) and (4) permit the determination 
of the exhaust velocity ratio corresponding to flow 
with and without equilibrium with respect to in- 
ternal vibration. 


(4) 


B. Calculation of Exhaust Velocities 
from the Enthalpy Theorem 


The condition of conservation of energy may be 
expressed by the enthalpy theorem 


H+ }3u?=E+pv+ }u?=constant, 
or 


. 


T 
f C,dT+4u*=constant. (5) 
0 


Equation (5) holds for equilibrium partition of 
energy. If equilibrium is not maintained with re- 
spect to internal vibrations, then Eq. (5) should be 
replaced by the relation 


T Tvib 
f caT+ f CyindT+ 3u?=constant. (6) 
. Sr 


A relation between ¢ and c’ may be obtained by 
applying Eqs. (5) and (6) to the two types of flow 
under consideration. The result is 


” eats Te’ 
f CdT+3c?= f C,dT 
0 0 


Te 
+ f CundT+4e"2, (7) 
Te’ 
or 
Te T.’ 
a—cr=a f CyindT+ ca | (8) 
- , 


. Te 
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Reference to Eqs. (3) and (8) shows that Eq. (3) 
determines the velocity ratio, whereas Eq. (8) per- 
mits calculation of the difference between c and c’. 
The two expressions together, therefore, permit 
calculation of the absolute values of both ¢ and c’. 
This fact will not be particularly useful in practice, 
however, since the numerical value of ¢ can be 
calculated in the usual manner! from the data re- 
quired to determine C,, C,’, AH7.7*, and AH7,’T¢. 


APPLICATION TO ROCKET MOTORS 


Representative calculations of c and c’ have been 
carried out and are described in this section for the 
hydrogen-fluorine and hydrogen-oxygen motors 
employing the reactants in stoichiometric propor- 
tions. Under these conditions it is justified to equate 
Cyiv, as a first approximation, to the mean vibra- 
tional heat capacities of hydrogen fluoride and 
water, respectively. 

The calculations are illustrated in detail for the 
hydrogen-fluorine rocket motor if chemical equi- 
librium exists during flow through the rocket nozzle. 
The required thermodynamic data were taken from 
unpublished work by M. Farber and are reproduced 
here: T,.=4915°K; T.=3640°K; c=11344 ft./sec. 
= 3.46 X 10° cm/sec. ; C, = 19.82 cal./°K/mole; Cyi» 
for HF =2.23 cal./°K/mole at 4915°K; Cyi» for 
HF = 1.90 cal./°K/mole at 3640°K. From the pre- 
ceding data it can be seen that Cy j,~2.06 cal./°K/ 
mole and C,’~19.82—2.06=17.76 cal./°K/mole. 
From Eq. (4) it is found that T,’~3515°K. There- 
fore AHr,7-~19.82 X (4915 —3640) =25,270  cal./ 
mole and AHr,’T:~17.76X (4915 —3515) = 24,864 
cal./mole. Substitution of the calculated enthalpy 
changes in Eq. (3) leads to the result c/c’ = 1.008, or 
c’=11,253 ft./sec. Therefore, the maximum differ- 
ence between c and c’ is seen to be less than one 
percent in the present case. Application of Eq. (8) 
instead of Eq. (3) is found to lead to the same result. 
Calculations for the hydrogen-fluorine motor have 
also been carried through for the case where flow 
through the nozzle occurs without composition 
change. In this case C.yi» represents nearly 20 per- 
cent of C,. The ratio c/c’ is now found to be in- 
creased to 1.034. 

Calculations similar to those outlined for the 
hydrogen-fluorine rocket with chemical equilibrium 
during flow through the nozzle have been carried out 
for the hydrogen-oxygen system using unpublished 
calculations by R. B. Canright. The result is 
c/c’ = 1.006. The conclusion may therefore be drawn 
that the state of vibrational energy distribution 
during flow through the nozzle is of no practical im- 
portance as regards the performance of the hydro- 
gen-fluorine and hydrogen-oxygen motors if chemical 
equilibrium is maintained during flow. For both of 
the rocket motors considered, the exhaust velocity 


is lowered by less than one percent if the vibrational 
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states remain frozen at the chamber temperature. 
On the other hand, if constant-composition flow 
occurs through the rocket nozzle, then the exhaust 
velocity may be lowered by as much as 3.3 percent 
for the hydrogen-fluorine motor. On the basis of 
previous work?’ it appears probable that equi- 
librium is maintained during flow with respect to 
atomic reactions but that composition changes as 
the result of molecular reactions are less likely to 
occur. The actual correction to the exhaust velocity, 
corresponding to completely frozen vibrational 
energy states, may therefore be intermediate be- 
tween the small values obtained for equilibrium flow 
and the larger values determined for constant- 
composition flow. 

Calculations for other rocket motors indicate that 
the changes in exhaust velocity associated with 
frozen vibrational states are generally of the same 
order of magnitude as the results calculated for the 
hydrogen-fluorine rocket motor. For example, for 
the red fuming nitric acid-aniline rocket, at a weight 
ratio of nitric acid to aniline of 2.80, c/c’ is found to 
have a value of 1.039 for equilibrium flow and 1.056 
for constant-composition flow. 

Reference to Eqs. (3) and (8) as well as to the 
calculated results of c/c’ allows some qualitative 
conclusions regarding the effect of frozen vibrational 
energy states on the exhaust velocity. In general, 
c/c’ will be large for systems in which Cy, consti- 
tutes a large fraction of C,. The numerical value of 


7S. S. Penner, J. Am. Chem. Soc. 71, 788 (1949). 


VOLUME 20, MAY, 1949 


Cyi» is nearly the same for flow with and without 
chemical equilibrium. On the other hand, C, is 
always larger for equilibrium flow than for constant- 
composition flow. It would therefore be expected 
that ¢/c’ is larger for constant-composition flow 
than for equilibrium flow. This conclusion is borne 
out by the numerical results given earlier in this 
discussion. 

The ratio c/c’ depends on the chamber tempera- 
ture, chamber pressure, and exhaust pressure. This 
can be shown by writing Eq. (3) in the form 


c/c' ={C,(T.—T.)/Cy'T [1 — (pe/pe)®/*>" J}. 


The vibrational heat capacity Cy, is relatively 


‘large for gaseous mixtures containing polyatomic 


molecules. Thus the red fuming nitric acid-aniline 
system contains CO, and H,O with sizable vibra- 
tional heat capacities. As the result, c/c’ is found to 
be greater for the red fuming nitric acid-aniline 
rocket than for the hydrogen-fluorine rocket. In 
conclusion, it should be pointed out that the heavy 
polyatomic molecules such as COz, which have large 
vibrational heat capacities, may also be the mole- 
cules with relatively long relaxation times which are 
most likely to lead to flow with frozen vibrational 
energy states.*® In other words, in some of the 
systems for which c/c’ is largest the exhaust velocity 
is likely to decrease from c to c’. Just how large this 
decrease actually is can be determined only after a 
complete study of the flow process and of the 
relaxation times in gaseous mixtures has been made. 
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A Note on the Substitution Principle for Steady Gas Flow 


R. C. Prim III 
Naval Ordnance Laboratory, White Oak, Maryland and Princteon University, Princeton, New Jersey 


(Received November 8, 1948) 


The Substitution Principle established by Munk and Prim for flows of a gas having constant 
specific heats is extended to the wider class of fluids having a state equation of the form p= P(p)S(s). 
The equations of motion for this class of fluids is reduced to a canonical form involving only the re- 
duced velocity vector and the pressure. It is shown that explicit elimination of the pressure from these 
equations is possible only for those fluids whose state equation has the form p= p*S(s) with k<1. 
Relations between the Mach number and stagnation pressure and the reduced velocity field are dis- 


cussed for the class of fluids under investigation. 





INTRODUCTION 


HE following Substitution Principle was es- 

tablished! by Munk and Prim for steady flow 

in the absence of extraneous force fields of a non- 

viscous, thermally non-conducting gas with con- 
stant specific heats: 


If in any flow field the densities are multiplied by m and 
the velocities by 1/(m)*, where m is a positive function which 
is a constant along each individual streamline, the new flow 
field will be a possible flow of the same gas having the same 
pressure distribution and Mach number distribution as the 
original flow. 


This Substitution Principle is particularly useful 
in clarifying the study of rotational gas flows—that 
is, flows in which the entropy and/or the total flow 
energy vary from streamline to streamline. For 
example, the function m can always be so chosen 
that the total flow energy of the derived (or “‘sub- 
stitute” flow) will be uniform throughout the flow 
field. Or, alternatively, in a region between two 
successive shock fronts an isentropic flow can be 
substituted for a non-isentropic one. Hence ques- 
tions involving flows which are both non-isoener- 
getic and non-isentropic can be resolved by the 
study of isoenergetic flows alone. 

Furthermore, by the choice m =a’, where a is the 
ultimate velocity magnitude along a streamline, we 
obtain the canonical formulation?’ of the equations 
of gas dynamics in terms of the reduced velocity 
vector (w=i/a) and the pressure alone: 


W-gradw+(y—1/2y)(1—w”’) grad Inp=0 (1) 
and 
div[ (1 —w*)!/7"9]=0, (2) 


where p denotes the pressure and y the (constant) 
ratio of specific heats. The condition of integrability 


1M. Munk and R. Prim, “On the multiplicity of steady 
gas flows having the same streamline pattern,’’ Proc. Nat. 
Acad. Sci. 33, 137-141 (1947). 

2M. Munk and R. Prim, “On the Canonical form of the 
equations of steady motion of a perfect gas,” Naval Ord- 
nance Laboratory Memorandum 9169 (June 1947). J. App. 
Phys. 19, 957 (1948). 

+B. Hicks, P. Guenther, and R. Wasserman, ‘‘New formu- 
lations of the equations for compressible flow,’”’ Quar. App. 
Math. V, 357-361 (1947). 
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for the pressure field is derived from (1) in terms 
of the w field alone: 


ww Xcurlw 
curl [| =0. (3) 


1—w" 


The set of Eqs. (2) and (3) then represent necessary 
and sufficient conditions involving the w field alone. 
They have proved particularly useful for the pur- 
poses of discovering families of exact solutions and 
establishing general geometric and kinematic prop- 
erties of rotational flow fields.*~* 

The purpose of the present note is to determine 
to what extent these results can be extended to 
flows of fluids having a more general state equation 
than that of a perfect gas 


p=Kp'!" exp(—s/Ccp). (4) 


(p denotes the density, s the specific entropy, and 
Cc» the (constant) specific heat at constant pressure. 
K is a dimensional constant.) In particular, it will 
be determined (a) for what equations of state 
p=p(p, s) the Substitution Principle is valid and a 
consequent canonical formulation of the equations 
of flow in terms of the pressure and reduced velocity 
fields alone can be obtained; and (b) for which of 
these equations of state the pressure can be elimi- 
nated from these canonical equations. 


THE SUBSTITUTION PRINCIPLE 


The equations governing the steady flow in 
absence of extraneous force fields of a compressible 
fluid without viscosity or thermal conductivity can 
be written in the form 


divpi =0 (5) 


4P. Neményi and R. Prim, ‘“‘Some geometric properties of 
plane gas flow,”’ J. of Math. and Phys. XXVII, 130-135 (1948). 

5R. Prim, “On the existence of steady gas flow in plane 
isothermal streamline patterns,” Bull. Am. Math. Soc. 54, 
1009-1012 (1948). 

6 P, Neményi and R. Prim, ‘Some Properties of Rotational! 
Flow of a Perfect Gas,” Proc. Nat. Acad. Sci. 34, 119-124 
(1948). 

7 R. Prim, ‘Extension of Crocco’s theorems to flows having 
non-uniform stagnation enthalpy,’’ Phys. Rev. 73, 186 (1948). 

8 R. Prim, “On a family of rotational gas flows,” Quar. App. 
Math. VI, 319-325 (1948). 
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d-gradd+(1/p) gradp =0 (6) 
d-grads =0 (7) 
p= pl, S). (8) 


In order to establish the Substitution Principle 
we shall introduce the “‘kinetic flux’’ vector g de- 


fined by 
q=(p)4o. (9) 
In terms of this vector, (5) and (6) can be written 
divg+3(q-grad Inp) =0 (10) 
and 
q:gradg—3(q-grad Inp)g+gradp=0. (11) 


Also, since from (8) 


0 Inp 0 Inp 
grad inp=( ) gradp+(—*) grads, (12) 
Op /. Os J» 


(7) can be written 





0 Inp 
(q-grad inp) -(—*) gG-gradp=0. (13) 
Op /, 

Equations (10), (11), and (13) are clearly invariant 
under any change of dependent variables leaving 
q, p, and (@-gradInp) unchanged provided that 
(0 Inp/dp), is a function of p alone. This permitted 
change of dependent variables is of the form 


p*=p; p*=mp; o*=(1/(m)')d, (14) 
where m is a constant along each individual stream- 
line (6-gradm=0). The provision that (0 Inp/dp), 
be a function of p only limits the validity of this 
substitution to fluids having a “‘product”’ equation 
of state; that is, an equation of state (8) of the form 


p=P(p)S(s). (15) 


The invariance of the Mach number under the 
substitution (14) follows from the identity 


0 Inp cin M? 
- ¢ 
for fluids having a celia equation of state. 
The validity of this substitution is not affected 
by the presence of shock discontinuities in the flow. 
The familiar equations relating the flow quantities 


on the two sides (1 and 2) of a shock front can be 
written in the form 





(16) 


(p1) 4g, sina = (p2) go Sina 


pitqi sin’a: = p2+ 2’ sin*a2 

qi° Pl dp q2” P2 dp 
=i 42 f 

Pi 0 p(P, $1) pe 0 p(p, Se) 


(a denoting the angle between flow direction and 
shock) in which it is immediately evident that the 
substitution (14) is permissible. 
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Hence, the Substitution Principle stated in the 
Introduction for flows of a perfect gas is valid in the 
same form for flows of any fluid having an equation 
of state of the form p=P(p)S(s). 

The Principle is also valid for flows of a non- 
homogeneous incompressible fluid. (In place of (13) 
we have g-grad Inp=0. The Mach number is then 
of course zero throughout.) 


THE CANONICAL EQUATIONS 





From Eq. (6) we have the Bernoulli or energy 
equation 
v’+2h(p, s) =a? =2h(po, s) (17) 
wherein the enthalpy h is defined by 
p dp 
h(p, s)= ; (18) 
0 p(P, s) 


and the ultimate velocity @ and the stagnation 
pressure po are constants along each individual 
streamline in a region free from shocks. (The 
quantities a and h(po, s) are unchanged by a shock 
front; but po and s are changed discontinuously.) 

We now set about the elimination of the redun- 
dancy of dependent variables which the Substitu- 
tion Principle reveals by introducing the reduced 
velocity vector # through 


w=i/a (19) 


(equivalent to the choice m =a’). 
The dynamic Eq. (6) then becomes 


1 
w-gradw+— gradp=0 (20) 
pa? 
and (17) 
h(p, s) =3a*(1—w*) =h(po, s)(1—w*). (21) 
By the use of (21), (20) can be written 
gradp 
w-gradw+ 3(1—w*) =0. (22) 





p(p, s)h(p, s) 
Now for the product equation of state (15) the term 


p(p, s)h(p, 





(23) 


is a function of p alone. (It can further be shown 
that only for the product equation of state is this 
the case.) Hence, letting 


1 
(=f ——; Pe)=-———_._ (aa) 
0 P(p) (dH(p)/dp] 
we can now write (22) in the canonical form 


wW-gradd+3(1—w*) grad InH(p)=0 (25) 
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involving W and p alone. Also (21) becomes 


H(p) = H(po)(1—w"), (26) 


and from (25) and (26) we obtain the important 
relation 


w Xcurlw 
grad In (po) = ——e, (27) 


connecting the W and pp fields. The integrability 
condition for (25) is the same as that (3) for a 
perfect gas. 

The introduction of the reduced velocity vector 
into the continuity equation (5) yields 


divpw =0 


which for the product equation of state takes on 
the canonical form 


divP(p)i=0. (28) 


In (25) and (28), then, we have a formulation of the 
equations (Eqs. (5)—(8)) of gas dynamics in terms 
of w and p alone for any fluid with a product 
equation of state. (For a non-homogeneous in- 
compressible fluid, (10) and (11) immediately yield 
such a formulation in terms of 9 and p). Given a p 
field and a w field satisfying (25) and (28), the 
stagnation pressure po is determined by (26). The 
ultimate velocity a can then be assigned arbitrarily 
for each streamline to obtain a velocity field 3. The 
entropy is then determined by 


S(s) =(2H(po)/a?*], (17a) 


and the density by (15). As will be shown later, the 
Mach number is given by (34). For a shock-free 
region, the entropy s may be assigned on each 
streamline instead of the ultimate velocity a, and 
the latter computed from (17a). 


ELIMINATION OF THE PRESSURE FROM 
THE CONTINUITY EQUATION 


We have shown that any product equation of 
state leads to an integrability condition (3) in 
terms of w only. We now ask which product equa- 
tions of state permit the elimination of P(p) from 
the continuity equation (28) (by the use of (25)) 
to obtain a form of this equation involving w# only 
(such as (2) for perfect gases). 

To answer this question we rewrite (28) in the 
form 
. d inP(p) 

divw+————# - grad Inf7(p) =0. (29) 
d |InH(p) 


The factor (wW-grad InH(p)) can be expressed in 
terms of w by the use of (25). If the resulting form 
of (29) is not to involve the pressure, it is then 
necessary that 


[d InP(p)/d InH(p) ]=constant. 
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It can be readily shown from (24) that this implies 
that the function P(p) is of the form 


P(p)=Cp*, (k<1), (30) 


where C and & are constants. 

Hence, the most general product equation of state 
permitting the elimination of the pressure from the 
equations is 


p=p*S(s). (k<\1). (31) 


The equations obtained are the same as those 
((1) and (2)) for a perfect gas with y set equal to 
1/k. Given a W field satisfying (2) and (3) (with y 
replaced by 1/k), a consistent pressure field can be 
obtained by integration from (1). The other flow 
quantities can then be determined as described 
above for the more general product equation of 
state. It will be observed that for the exponential 
state equation (31), the pressure field can be multi- 
plied by any positive constant without violating 
the differential equations, but that this is not the 
case for the general product equation of state (15). 

It should be noted that the variety of physically 
meaningful state equations represented by (31) 
is not substantially greater than the family (4) of 
state equations for perfect gases with 1<y<5/3, 
because such considerations as the positiveness of 
temperature and of internal energy preclude the 
existence of substances with state equations of the 
form (31) with k<0, S(s) <0, or S’(s)>0. 


RELATION BETWEEN MACH NUMBER AND 
REDUCED VELOCITY 


It is well known that for a perfect gas the Mach 
number and the magnitude of the reduced velocity 
vector are substantially equivalent quantities con- 
nected by the relation 


M? = (2/7—1)(w?/1—w*). (32) 


For the general equation of state p=p(p, s) the 
Mach number can by use of (19) and (21) be ex- 
pressed as follows 


M?=[4p(p, s)/dp }e” 
=2[dp(p, s)/dp jh(p, s)(w?/1—w*)] (33) 
= M*(w’, p, s). 


For the product equation of state, (33) can be 
written 


M?=2[d InP(p)/d InH(p) ](w?/1—w®) 
= M*(w*, p), 


and we see that only for the power law state equa- 
tion (31) does (33) take on the form 


M?= M*(w*). 


(34) 


(For this case 
M?=[2k/1—k ][w?/1—w*].) (35) 
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Stabilization of Breakdown Voltage and Measurement of Some 
Critical Currents in a Glow Tube 


A. B. PAL 
Depariment of Electrical Engineering, Princeton University, Princeton, New Jersey 


(Received January 14, 1949) 


This paper describes a special technique by which it is possible to stabilize the breakdown voltage 
of a glow tube having zirconium electrodes in pure argon at a pressure of 8.5 mm Hg. A brief account 
is also given of some experimental methods by which it is possible to find the order of magnitude of 
two critical currents, whose concepts have been used in the theory of spark discharge at low pressures. 
One of these critical currents is the “threshold current,”’ which is defined as the maximum steady- 
state current in a glow tube just before a spark. The other critical current is the maximum electronic 
current (attained during the development of a spark) at cathode before space-charge distortion 
becomes appreciable. The methods described below were confined to two tubes having zirconium 
cathodes of area 0.3 sq. cm, the electrode spacing being 3.5 mm. One of these two tubes had a zir- 
conium anode, and the other one had a thorium anode. The tubes contained pure argon at a pressure 
of 8.5 mm Hg. For both the tubes, the order of the: ‘threshold current” was found to be 10-* amp., 
and the order of the other critical current (electronic current at cathode) was estimated as 10-* amp. 





STABILIZATION OF BREAKDOWN VOLTAGE 


N order to produce a spark discharge in a glow 

tube, the presence of some initial ionization 
current 79 is necessary. Experiment shows that the 
“time lag’’ may be several seconds, if the initial 
ionization current ip <10-'4 amp. If the value of io 
is fixed and kept much above the value 10~'* amp., 
then the breakdown voltage will have a definite 
value, and the value of the time lag will be almost 
zero. It is found that a 100-w lamp placed near the 
tube with a zirconium cathode produces an initial 
ionization current whose value is of the order of 
10-" amp., and so the breakdown voltage assumes 
a definite value. We may define the breakdown 
voltage as the value of applied voltage which pro- 
duces a spark discharge almost immediately, when 
the initial ionization current has a given value, say, 
10-"! amp. In defining the breakdown voltage, the 
arbitrary value of 7) must be stated, as the value of 
the breakdown voltage depends upon the value of io. 

The value of the breakdown voltage E, at low 
pressures of the order of a few mm. Hg is mainly 
governed by the following three factors: 


(1) Nature and purity of the enclosed gas. 
(2) Nature of the cathode surface. 
(3) Initial ionization current io. 


If these three factors are kept unchanged, the value 
of the breakdown voltage E, will also remain un- 
changed. In actual practice, it is very difficult to 
keep the first two factors unchanged. A heavy dis- 
charge may produce impurities and also change the 
nature of the cathode surface, resulting in a different 
value of E, for the next spark. Continuous glowing 
of the cathode is also undesirable, as it may heat 
the cathode and so change its work function. Such 
a heating may lower the breakdown voltage by 
several volts. Thus, in an actual test, a tube with 
zirconium cathode had an Ey equal to 157 volts at 
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room temperature. On glowing the cathode with a 
d.c. current of about 50 ma for a few minutes, the 
cathode became red hot, and the value of EZ fell 
to 112 volts, i.e., the fall in E, was 45 volts. On 
cooling the tube to room temperature, the value 
of E, was restored to the original value 157 volts. 

In order to control the first two factors, we must 
select a suitable cathode material and a suitable gas 
for filling the tube. A suitable material for the 
cathode is pure zirconium, and a suitable gas is pure 
argon. Impurities must be removed by repeated 
induction heating of the electrodes, while the tube 
is connected to a pump. The tube must be baked 
for a long time. There must be long traps between 
the tube and the pump cooled by liquid air, in 
order to prevent mercury vapor and other im- 
purities from entering the tube. When the tube is 
ready, its cathode surface requires further cleaning 
in order to make the surface stable for a long time. 
This requires intermittent bombardment of the 
cathode surface with positive ions. The bombard- 
ment must occur in almost instantaneous blows, 
repeated after sufficiently long intervals, so that 
heat produced by bombardment radiates away 
during the vacant intervals and does not penetrate 
into the body of the cathode. The surface cleaning 
may last for several hours before it is satisfactorily 
complete. The writer got very good results by 
continuing this process for as long a time as 100 
hours! The breakdown voltage then becomes con- 
stant and no change is detected in its value even 
after two months. 

Figure 1 shows a very efficient circuit for cleaning 
the cathode surface. It enables us to control the 
discharge current, keep its value almost constant 
during a discharge, control the duration of a dis- 
charge, and also control the interval between two 
successive discharges. This circuit thus meets all 
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the requirements of the cleaning process. Its action 
may be explained as follows: 

The condenser C acquires a d.c. voltage (about 
500 volts) much above the breakdown value £, for 
the tube T and so a permanent glow current flows 
through the tube. When S is open, the glow current 
is very small due to the presence of the resistor R, 
which has a resistance equal to 0.1 megohm. When 
S is closed, the value of the glow current rises to a 
high value, which is controlled by the tapping 
point on the resistor 7, which has a value of 1000 
ohms. The gap S is mechanically opened and closed 
by the vibrating lever L (attached to an ordinary 
electric bell) and the rugged terminals m and m. 
The period of vibration of the lever L controls the 
interval between two successive heavy discharges. 
The duration of a heavy discharge is controlled by 
adjusting the position of m, which changes the time 
of contact between m and m. It can be shown that 
the value of the heavy discharge current during the 
time of contact between m and m remains approxi- 
mately constant, if the value of C is about 100 uf 
and the time of contact is about 1/1000 second. 

When the cathode surface becomes clean, the 
glow on its surface shows uniform brightness at 
all points. Unclean surfaces give glows whose 
brightness varies from point to point on the cathode 
surface. Thus, the appearance of the glow also 
helps detect any impurities in the surface of the 
cathode. 


MEASUREMENT OF THE FIRST CRITICAL CURRENT 


The concept of this critical current (also called 
the “threshold current’’) has been used by Loeb.! 
We shall denote this current by 7, (Loeb denotes 
it by 2). The writer has measured the threshold 
current 7, by several methods, and all these methods 
agree in giving the order of magnitude of this 
eritical current as 10-* amp. for the two tubes used 
by the writer. We shall briefly describe these 
methods one by one. 


First Method 


Suppose 7, is initial ionization current produced 
by some light illuminating the zirconium cathode. 


1L. B. Loeb, Fundamental Processes of Electrical Discharge 
in Gases (John Wiley & Sons, Inc., New York, 1939), pp. 420-423. 
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There is also some current 7, due to initial ionization 
produced by stray radiations, such as cosmic rays. 
The value of i, changes at random. The total 
current due to initial ionization, i.e., 7%, is thus 
given by 


19 =1p+1>. (1) 


Experiment shows that 7 has large random fluctua- 
tions, if i, < 10-'* amp. These fluctuations are due to 
the existence of the random component 7, As we 
raise the value of i,, the percent random variation 
in to decreases. The random changes in 7 produce 
random changes in £). Experiment shows that if 
we illuminate the zirconium cathode by direct sun- 
light at noon, then the value of 7, is of the order of 
10-® amp. With such a high value of 7,, the percent 
random fluctuation in 79 becomes almost nil, and 
the value of E, remains constant. The value of time 
lag for such a large value of io is also almost nil. 

Now, consider the circuit shown in Fig. 2, in 
which T is a tube having a zirconium cathode 
illuminated by direct sunlight at noon, and J is a 
known current, which can be set at any desired 
value in the range 10~-* to 3X10-* amp. The current 
I can be produced by a device using a photo-cell 
exposed to light of controlled intensity. 

The current J charges the condenser C. Let the 
voltage built up in C be denoted by E. When E is 
well below the breakdown value E,, then the tube 
current i<J, and the rise of E with time ¢ is almost 
linear. As E-—E,, the tube current approaches the 
critical value 7,, which we want to measure. Sup- 
pose J>>i,. As the tube current 7 (before breakdown) 
can rise only to the maximum value 7,, the charging 
current (I—i) remains almost equal to J, since 
I>i,. The voltage E thus rises almost linearly with 
time ¢ from E; (extinction voltage) to E,, and we 
get a series of discharges, which are heard as a 
series of clicks in the sounder. The interval T 
(between two successive clicks) is given by 


T=C(E,—E,)/I, since >i. (2) 


As we lower the value of J, so that ultimately 
I—i,, then (2) cannot hold, since the value of the 
charging current (I—7) drops to a negligible value 
as i—%, near the breakdown condition. Suppose E 
rises to Eo (which is slightly below £,) in a time To, 
the value Ey being such that for E< Eo, the tube 
current 7 remains negligibly small. Let the time 
taken by E to rise from E» to E, be ta. The total 
time T (between two consecutive clicks) is thus 
given by 

T = To+ta=C(Eo—E,)/I+ta. (3) 


If I>i., then To>ta, and T~T>. As I—%,, the 
value of ta. Thus, if we gradually lower the 
value of J from a large initial value, and notice 
the value of I when T becomes very large, then this 
value of J will be roughly equal to 7,. In this way, 
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it is possible to find an approximate value of %,. 
Experiment shows that the interval between two 
clicks becomes very large for values of J around 
2X10-* amp. This shows that the order of magni- 
tude of the current 7, is 10-* amp. 


Second Method 


If the zirconium cathode of a glow tube be 
illuminated by a lamp of 100W instead of direct 
sunlight at noon, we notice appreciable random 
fluctuations in 7. As the value of E, depends upon 
the value of zo, there will also be some noticeable 
fluctuations in E,. In this case, the value of ¢g in 
Eq. (3) will depend upon chance. In other words, 
the value of tz will be erratic. For values of [>>z,, 
the value of ta&T , and so the total interval T will 
be approximately constant. But if we lower the 
value of J, so that J remains only slightly above 
the critical current 7,, then the value of the erratic 
component tg will be very large, and this will make 
the total interval T also erratic in value. Measure- 
ments of the interval ¢ for different values of the 
capacity C have shown that for values of J around 
2X10-* amp. the interval tg becomes very large 
and erratic. Hence, the ‘value of 7, is roughly equal 
to 2X10-* amp., and its order of magnitude is 
definitely 10-* amp. 


Third Method 


The principle used in this method will be clear 
from the circuit shown in Fig. 3. 

Suppose the d.c. supply voltage Eo is kept several 
volts above the breakdown value FE of the tube T. 
When S, is closed and S; is open, a glow current 
passes through the tube 7. The value of this glow 
current depends upon the resistance R. If we open 
the switch S:, the glow current will charge the 
condenser C. In this way, a voltage is built up in C 
and as the value of this voltage rises, the value of 
the glow current falls, and ultimately the current 
drops to a value which is almost zero. Suppose the 
voltage on C rises from 0 to x volts in a time f¢. 
The value of x can be measured as follows. After a 
time ¢, when the voltage on C acquires the value x, 
we suddenly lower Ey to 0. If we again raise the 
value of Eo, we shall get a spark discharge in T 
(heard as a loud click in the sounder) for some 
value E,’ of Eo, such that 


E,’ =Eyxtx or x = FE,’ —E,. 
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As we know E£,’ and E&, we can at once find the 
value of x. We can measure x for different values 
of ¢, and then plot x vs. ¢t. From this plot we can 
find the value of dx/dt, and hence of the tube 
current (Cdx/dt), for any value of t. As ¢ increases, 
the value of the tube current falls, and ultimately 
drops to a value which is almost zero. As the value 
of the tube current for tube voltage nearly equal 
to E, is 7, and the tube current drops to almost 
zero for a slight fall in value from E,, the minimum 
measurable value of the tube current in the plot of 
tube current vs. time must be of the order of %,. 
The writer measured x for several values of ¢, and 
plotted x vs. ¢ for different values of R, such as 22, 
44, 66, 88, and 110 megohms. The value of C was 
0.5 wf. From these plots, the values of tube current 
were computed for different values of t. These 
values showed that the tube current had a minimum 
value of the order of 10-* amp. before suddenly 
dropping to a value which was almost zero. Hence, 
the order of magnitude of the critical current 7, is 
also 10-* amp. 


Fourth Method 


For this method, we use the circuit arrangement 
shown in Fig. 4. 

Here Risa very high resistance—several thousand 
megohms! The writer used the wooden spacing be- 
tween two terminals fixed on a wooden block as a 
high resistance, and found its value quite suitable. 
The value of R changes with time, but remains 
pretty constant for the short interval during which 
we use R in our experiment. C is simply the stray 
capacity. As we raise the supply voltage above E&,, 
the tube may be in the breakdown condition. But 
the tube current is kept low by the high resistance 
R, and unless it rises to the critical value 7,, a break- 
down cannot occur. We can increase the value of 
tube current, however, by raising the supply volt- 
age. So we increase the supply voltage to a value E,, 
when a spark discharge is produced and is heard as 
a click in the sounder. The energy is supplied by 
the capacity C, which is discharged. Just before the 
breakdown, the voltage across the tube is nearly 
equal to E,, and the tube current is of order 7,. 
Hence, we have the simple relation: 
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In an actual experiment, the following values were 
observed for a tube. 


E.=263 volts, 


E, = 154 volts. 
Putting these values in Eq. (4), we get 
i,=109/R. (5) 


To find R, we set C=5000 uyf, and keep the supply 
voltage equal to E.. We note the time T taken by C 
to build up its voltage from 0 to E,, when a dis- 
charge will pass through the tube. The following 
values were observed for 7. T=17.2, 17.0, 16.8, 
17.0, 17.1. Mean value of T=17.0 sec. These 
values clearly indicate that R remains pretty con- 
stant for several minutes, and its value is quite 
reliable for our purpose. Now, we use the relation: 


E=E(i—e-"/®°), 
At t=T, E becomes equal to E;. Thus 
E,y=E(i1—e-T!/®°), 


In this equation, we know the values of all the 
quantities, except R. We can, therefore, compute 
the value of R. Calculation shows that R=3.8X10° 
ohms. Putting this value of R in Eq. (5), we get 
i,=2.9X10-* amp. The order of magnitude of 7, is 
therefore definitely 10-* amp. 


Fifth Method 


Consider the circuit shown in Fig. 5. As before, 
I is a current which can be set at any value in the 
range 10-* to 3X10-* amp. by means of a suitable 
device. When the switch S is closed, the voltage 
across C2 remains 0, and the tube breaks down each 
time the voltage built up in C, exceeds the value E 
for the tube 7. We thus get a series of discharges, 
which produce clicks in the sounder at regular 
intervals. The interval between two successive 
clicks depends upon the values of J, Ci, E», and E; 
(the extinction voltage). If J>>7,, then this interval 
is C\(E,—E,)/I. 

What happens to the discharges, if S is opened? 
Experiment shows that if S is opened just after a 
discharge, the next discharge may not occur under 
certain conditions for values of J below a certain 
critical value J. 

Let us see how Jp is related to 7,, the threshold 
current. As the voltage across T reaches F,, the 

‘ tube current rapidly builds up to 7., just before the 
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spark. Thus, when the voltage across T is a little 
below E,, the value of the tube current is an 
appreciable fraction of z,. The tube current indi- 
cated in Fig. 5 by 72 is thus of the same order as %,. 
Hence, the determination of 72 gives a value whose 
order is the same as that of 7,, in which we are 
interested. 

Now, let us see how it is possible to measure 1%., 
As the resistance in series with the tube T is 1K 
(or 1000 ohms), the voltage drop across it by a 
current of the order of 10-* amp. is almost nil, and 
so we shall ignore it. Let the instantaneous value of 
the voltage built up in C,; be E,, and that in C, 
be E2. The voltage acting across the tube is then 
(E,—E:). If the circuit Components are so chosen 
that (for i2<%,) 


dE,/dt =dE:/dt, (6) 


then the value of (Z,—£,) will remain constant. 
Notice dE,/dt is almost nil, when (£,— £2) is well 
below the breakdown value £&, since the tube 
current 72 becomes appreciable (and of the order 
of i.) only when the applied voltage (E:— £2) be- 
comes close to the breakdown value EF. Hence, 
condition (6) is satisfied only when 7 is of the order 
of z,. When this condition is satisfied, then the tube 
voltage (E,—E:2) remains permanently below &, 
and so we cannot get another discharge. 

The value of i2 can be determined in terms of Jo, 
C,, and Cs, if condition (6) is satisfied. It is obvious 
that 


i:=CidE,/dt, and 
Thus, by Eq. (6), we get 


12 = C.dE,/dt. 


41/C1=12/C2. (7) 
But [9 =71+72, and so we get the value of 72 given by 
42 = IpC2/(Ci+ C2). (8) 


The writer performed an experiment to measure 
iz by Eq. (8). Several values of C,; and C2 were 
tried. The value of J was gradually reduced till on 
opening S the next discharge was delayed by a long 
time. The value of J, when the discharge time be- 
comes very large, is then very nearly equal to the 
critical value J) in Eq. (8). (As mentioned before, 
we open S just after a discharge, which is one of 
the series of discharges which appear when S is 
kept closed.) In performing this experiment, the 
zirconium cathode must be illuminated by light, 
as 72 cannot flow continuously in dark. The writer 
used a lamp of 100w for this purpose. In dark, 
io (current due to initial ionization) fluctuates at 
random, and this makes iz also unsteady, so that 
Eq. (7) cannot be satisfied. Thus, with the tube 
having a thorium anode in dark, no critical value 
of I could be observed. On illuminating the zir- 
conium cathode with a lamp of 100w, however, it 
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its value to be of the order of 10-5 amp. According 

to the writer’s estimate, however, the critical 

current 7, seems to be of the order of 10-* amp. This 

may be shown as follows. As each electron starting 

from the cathode produces several electrons on its 

way to the anode, the current 7, assumes a higher 
value J, at the anode. We may set 


I,=1,p, (11) 


where p denotes a multiplying factor. Theoretically, 

it is hard to find the value of the factor p. It is 

possible, however, to estimate roughly the factor 

pb from some experimental values of currents at the 

cathode and at the anode which are not very differ- 

ent in order of magnitude from the currents i, 

and J, Experiment shows that the value of 1% 

(initial ionization current) for the tube with zir- 

conium electrodes exposed to direct sunlight at 

noon is of the order of 10° amp. When the applied 

voltage is just below &;, the current at anode has 

the value of order 10-* amp., as shown by the above 

five methods of measuring this current. If we assume 

the positive ions are not liberating appreciable 

number of electrons from the cathode, then the 

order of p is seen to be 10, since the electronic 

current at cathode is 10-® amp., and the electronic 

current at anode is 10-* amp. As a rough approxi- 

mation, we may set p in Eq. (11) equal to 10, 

since the space-charge distortion is assumed to be 

negligible upto the value 7,. In order to compute 4, 

from Eq. (11), we still need the value of the second 

unknown current J,, which is the minimum glow 

current flowing in the tube just after a spark dis- 

charge gets started. The value of this current may 

be measured experimentally by the circuit shown 
in Fig. 6. First, we close the switch S, and lower E, 
to 0 volt. Then slowly increase the voltage EF, so 
that voltage built up in C is nearly the same as Fj, 
the value of the time constant Cr being very small. 
When £,=E,, we get a spark discharge in the 
tube 7. Then we lower £, to 0, in order to restore 
the tube to normal condition, and then raise E, to 
a value slightly below the breakdown value &. 
Now, open S, and gradually raise the voltage Ep. 
As (E,+ E:) exceeds the value E,, the tube breaks 
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was possible to detect and measure the critical 
value Jo. 
As a result of several observations, the values of 

ig computed by using Eq. (8) were seen to be of 

the order of 10-* amp. 

It will be seen from the above discussion that 

all the above five methods agree in giving the order 

of magnitude of the threshold current 7, equal to 

10-§ amp. for the two tubes used by the writer. 

These tubes had zirconium cathodes of surface 

area 0.3 sq. cm, the anode-cathode gap being about 

0.4 cm, and the enclosed gas was argon at a pressure 

of 8.5 mm Hg. The anode of one of these tubes was 

of zirconium, and that of the other one of thorium. 

The electrodes were properly seasoned so that each 

tube had a definite and stabilized breakdown 

voltage. 
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MEASUREMENT OF THE SECOND 
To, CRITICAL CURRENT 


on The second critical current, which we shall denote 


by 7,, was introduced by Schade? in his theory of 
the formative time lag. According to this theory, 
the main part of the formative time lag is spent in 
raising the cathode electronic current from its 

initial value 7%» to the critical value z,. The space- 

charge distortion is assumed to be serious only 

when the cathode current exceeds the critical 

value i,. According to Schade’s theory, the forma- 

tive time lag T; is given by the equation: 


T;=(t;/€) Inf 1+ €(i,/%o) ], (9) 


where ¢; means the time taken by the positive ions 
to cross the anode-cathode gap, which is usually of 
the order of one centimeter. The value of ¢t; is of 

the order of 10-* sec. The quantity ¢ is a function 

of the over-voltage (E—E,), where E is the applied 
voltage, and E; is the threshold voltage. As E-E;, 

the quantity e—0. Hence, in case E is set at the 
threshold value, Eq. (9) reduces to the simple form 


T's = (t t2/t0), 1, = (Tfo/ti). (10) 


Schade verified Eq. (9) in a rough manner. His 
experiments can give no value for the critical 
current 1,. Schade regards the critical current 4,, 
however, as the minimum glow current, and assumes 


*L. B. Loeb and J. M. Meek, The Mechanism of the Electric 
Spark (Stanford University Press, California, 1940), pp. 12-20. 
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down. But the tube current (a weak glow discharge) 
is controlled by the high resistance R. As we in- 
crease E; still further, the glow current is increased. 
With the rise of the glow current, the voltage 
across the tube falls below £;, the diode circuit is 
at once closed, the condenser C sends a heavy dis- 
charge through the tube 7, and this discharge is 
heard as a loud click in the sounder. At this mo- 
ment, the voltage across T is roughly equal to &,, 
and as all the current flowing in the tube (just before 
the condenser C discharges) also flows through the 
resistance R, we have the simple relation : 


Glow current just after breakdown 


=(E,\+E:—E,)/R. (12) 


The current given by Eq. (12) is expected to be 
quite small, since the voltage across T is a little 
less than £,, and the glow current becomes appreci- 
ably high only when the voltage across T has an 
appreciable fall below Ey. As the exact value of the 
voltage across T just before the moment of hearing 
the loud click is not known, the value given by 
this method is only an approximation to the lowest 
glow current J,, which can possibly flow just after 
a breakdown. The following values refer to an 
actual experiment: 


(1) R=66 megohms, E,=225 volts, (E:+E2)=260 volts. 


Putting these values in Eq. (12), we get J,=5.3 
X10-7 amp. 


(2) R=110 megohms, E,=225 volts, (E, + £2) =280 volts. 


Putting these values in Eq. (12), we get J,=5.0 
X10-7 amp. This experiment shows that the order 
of magnitude of the minimum glow current J, in 
Eq. (11) may be roughly taken as 10-7 amp. As the 
order of the factor p in Eq. (11) has been shown to 
be 10, we see that the order of magnitude of the 
critical current 7, is 10-* amp. 

It will be noted that the value of 7, assumed by 
Schade, namely, 10-° amp., differs very widely from 
the writer’s estimate which is 10-* amp. In order 
to decide which estimate is the correct one, we may 
resort to experiment and measure the order of 7, by 
using Eq. (10), which is based on Schade’s theory. 
For this purpose, we have to measure the value of 
T; and to. The value of t; may be taken as 10~ sec., 


so far as the order of magnitude is concerned. The 


writer produced ip by illuminating the zirconium 
cathode of a tube with a lamp of power 15W, 
placed at a distance of about four inches from the 
cathode. In order to control the intensity of light 
(and hence the value of io) given by the lamp, the 
voltage across the lamp terminals was supplied 
through a variac, having 130 scale divisions. When 
the variac was set at 0 division, the cathode of the 
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tube was in perfect darkness. In this case, the 
value of JT; was several hundred secs. and also 
highly erratic, since the only source of 7) was stray 
radiations which change at random. On gradually 
increasing the intensity of light, the value of JT; was 
found to fall, and its different values for a given 
intensity of light (and hence a given value of 4p) 
became less erratic. With the variac set at 70 
divisions, the value of 7; became of the order of 
50 sec. When the variac was set at its maximum 
position, the value of 7; reduced to about 2 sec. 
The value of %) (corresponding to the maximum 
position of the variac) was still too small to be 
measured by a direct and convenient method. On 
reducing the distance of lamp from the cathode to 
1/4 of its previous value, the value of the intensity 
of light falling on the cathode was increased 16 
times, and the new value of 7) was thus made 16 
times greater than its previous value. This is, 
however, an approximation, as the lamp was not a 
point source of light. Now, the value of ip was 
5X10-" amp., which was determined by using the 
“second method” described in an article already 
published by the writer.* The value of iz» for the 
previous position of the lamp may, therefore, be 
taken as (1/16) K(5X10-"*) or 3X10~-™ amp. ap- 
proximately. The corresponding value of 7} was 
2 sec., as already stated. We are, now, in a position 
to compute the order of magnitude of 7, by using 
Eq. (10). Taking only the orders of magnitude of 
the quantities 7;, io, and ¢t;, we get at once the 
order of magnitude of i, equal to 10-* amp., which 
agrees with the writer’s estimate. Hence, the order 
of magnitude of the second critical current 7, is 


. 10-° amp., and not 10-° amp., as assumed by 


Schade. The validity of Schade’s estimate could 
not be checked by the rough manner in which 
Schade tested his theory. 

It will be noticed that the currents mentioned in 
this paper are total currents, and not current 
densities. As the electrode area is 0.3 sq. cm, the 
order of the current density is not very different 
from that of the total current. 
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High Frequency Total Emission Loading in Diodes* 


NicHoLas A. BEGOVICH 
Hughes Aircraft Company, Culver City, California 


(Received August 16, 1948) 


An expression is derived for the total emission conductance produced by the electrons returned to 
the cathode of a diode operating with a retarding off-cathode field. For the retarding voltage and . 
diode spacing used, the assumption of a linear potential distribution is justified. Examples of the 
variation of the total emission conductance with retarding voltage and frequency are illustrated. 





I. INTRODUCTION 


AN interesting high frequency loading phenom- 
enon occurring in very close-spaced diodes is 
reported by a number of investigators.'~* The ex- 
perimental set-up is briefly as follows. The ‘‘Q”’ of 
the cathode-anode cavity of the diode is measured 
as a function of the d.c. retarding voltage on the 
anode.‘ The frequency used for the measurements is 
3X10° c.p.s.5 At a very large d.c. retarding voltage, 
the “Q” of the cavity is the same whether the 
cathode is heated or not. With the cathode heated, 
the ‘‘Q”’ is measured as the retarding voltage is de- 
creased in magnitude. The ‘‘Q”’ sharply decreases 
and reaches a minimum value before any measur- 
able value of d.c. current is collected at the anode 
(i.e., before any appreciable number of electrons 
reach the anode). 

Clearly, the loading of the cavity must be pro- 
duced by the electrons that are returned to the 
cathode by the retarding anode field. On an average, 
because of the finite to and fro transit time, the 
electrons returning to the cathode must have an 
energy greater than their emission energy. This 
extra energy is supplied by the high frequency field 
and is converted into heat when the returning elec- 
trons strike the cathode. This loading has been 
named ‘‘total emission loading”’ since it is produced 
by the total emission of the cathode.' ® 

In the normal operation of vacuum tubes as 
amplifiers, the anode current is always only a small 
fraction of the total emission current. Consequently, 
there exists in front of the cathode a potential mini- 
mum and a resulting retarding off-cathode field. 
This retarding field drives most of the electrons 
(~90 percent) back into the cathode. If these 
electrons return to the cathode with an energy 


* The results of this calculation were reported at Pasadena, 
June 24-26, 1948, meeting of the American Physical Society. 

'C. N. Smyth, Nature 157, 841 (1946). 

2 A. van der Ziel and A. Versnel, Nature 159, 640 (1947). 

*D. R. Hamilton, J. K. Knipp, and J. B. H. Kuper, 
Klystrons and Microwave Triodes (McGraw-Hill Book Com- 
pany, Inc., New York, 1948), p. 156. 

* This experiment was also performed on the input cavity 
of a triode. (See reference 3.) 

*An exception to this is the van der Ziel measurements 
which were made at 4107 and 5.2 107 c.p.s. 

* A. van der Ziel and A. Versnel, Philips Research Reports 
3, 13 (1948), 
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greater than their emission energy, the input circuit 
of the tube will require a greater input driving power 


for a given power output. 


II. THEORY 


For the very close cathode-anode spacing used in 
these experiments and with retarding voltages of the 
order of 1 or 2 volts, the potential distribution be- 
tween the cathode and anode is very close to linear. 
This assumption is suggested by the experimental 
evidence of maximum loading before an appreciable 
number of electrons strike the anode of the diode. 
An additional justification for this assumption will 
be given in Section III. In the equations that follow, 
c.g.s.-practical units will be used. - 

Figure 1 shows the configuration of the diode 
consisting of two parallel infinite planes a (the 
cathode) and 6 (the anode). The spatial dependence 
for all quantities will be given by the x variable 
measured from the cathode, i.e., the problem is 
assumed one-dimensional. To calculate the energy 
of the returning electrons, the force equation must 
be integrated.’ 


mé = —107’eE 
= —107(e/d)(Vso— Vor cos{wit@}). (1) 


The stationary (d.c.) retarding voltage at the b 
plane is Vso; the time dependent voltage is V;. Note 
that the sign of the electron charge is written into 
the equation so that when e appears it will always 
mean the absolute value of the electron charge. At 
the cathode (x=0), the electron’s initial condition 
is V=¥Vo at time ¢=0. Integrating (1) and using the 
cited initial condition, the velocity at time ¢ is 


; Vou 
b=m—(/d)(Vot+—“fsind—sinfw+4I}), 2) 
Ww 
where /=10%e/m. The position at time ¢ is 
Van 
x =Vot — wa(s Voot? +-— 
w? 


X {wt sind+cos(wi+¢) —cos¢} ). (3) 


7 This method of calculating the transit time has been used 
by a number of writers. Compare C. J. Bakker and G. deVries, 
Physica 1, 1045 (1934); and K. S. Knol, M. J. O. Strutt, and 
A. van der Ziel, Physica 5, 325 (1938). 
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— 7 Fic. 1. Basic diode configuration. 
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The stationary transit time ft) (transit time when 
Vs1=0) for the to and fro path of the electron is 
given by one solution of (3) with x set equal to zero, 
the other solution being the trivial one t)=0. 


to = 209d /1 Vio. (4) 


With V,, different from zero, the transit time T can 
be written as 


T=tot+piitpri2t:--, (S) 


where p= V:/ Vio. Clearly, p is the proper parame- 
ter for the transit time expansion since the ratio of 

- Voi to Veo will have a direct influence on the pertur- 
bation of fo. Using (4) and (5) in (3) with x set equal 
to zero, expanding the resultant expression in powers 
of », and separately setting the coefficient of each 
one equal to zero, give 


5; = }(to/@) { (cos20—1) cosd 
—(20—sin26) sing} (6) 
and 


52 = — (5:°/to) + (61/6) { (sin26 cos 
—(1—cos26) sing}, (7) 


where 0=wt)/2 (the half-way transit angle). 

The velocity v, of the returning electron at the 
cathode will be given by substituting (4) to (7) in 
(2). Squaring the resultant expression, and multi- 
plying by m/2, will give the difference in the emitted 
and returning kinetic energy for an electron leaving 
the cathode with a velocity vp at a phase angle ¢ for 
the potential V,:. The average energy interchange 
per electron for electrons leaving the cathode with a 
velocity v» will be the integration of this expression 
over a complete period of the phase angle ¢. 


$m(v,? —v9") = m(vop/8)?(1—@ sin2@—cos26). (8) 


-The number of electrons per second per cm? (An/At) 

leaving the cathode in the velocity range v9 to 
VYotdvy is given by the Maxwell-Boltzmann ex- 
pression 


(An/At) = (Io/leVr) exp(—v0?/21Vr)vedv, (9) 


where J, is the emission current amps./cm*, and 
Vr=kT./10*’e k being Boltzmann’s constant 1.380 
X10-"* ergs/deg., and T, being the absolute cathode 
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temperature. The energy interchange per second per 
cm? between the cathode and the emitted and re- 
turned electrons in the velocity range vp to v9 +d» is 
given by multiplying (8) by (9). Dividing this 
product by 10’, the energy interchanged for the 
velocity packet dv is 


AP =(Io/Vr)(vop/10)2(1 —@ sin20 —cos26) 


watts 
Xexp(—v9?/21 V7) vedv9 





(10) 
cm 


Since the transit angle @ is a function of vo, it will be 
more convenient to write (10) completely in terms 
of 0. . 

Introducing the following parameters 
h = Voo"l(2w*d? Vr)—! = 1.021 - 10'9(wd)-? TT Veo”, (11) 


and 


a= wd (1 V5)", (12) 
(10) can be written as 
AP =2Iop*h(la?)“"(1 — 0 sin20—cos26)e—""6d@. (13) 


The critical emission velocity, the velocity of emis- 
sion which separates the electrons reaching the 
anode from those returning to the cathode, is 


v= 21| Veol, (14) 
or in terms of the transit angle 
6,=av,;=2.115- 10-7 fd(Vio)-3, (15) 


where f is the frequency of Vs:. When V,, is different 
from zero, (14) should be modified. However, as will 
be seen shortly this will be an unnecessary correc- 
tion. Using (15) for the upper limit in the integration 
of (10) (see Appendix A), and defining a total emis- 
sion stream conductance G, as 


P=3V,,°G, watts/cm’, (16) 
give 
41 kT, 


:= pean Caled (h+1) cos26, 
107e Vio? 


h 
+h(0, sin20;— 1)} —_ i+ 





h 
Xe-"" Fy(1/2; 3/2; 1/0) mhos/cm*, (17) 


where ;F; is the confluent hypergeometric function. 
A curve is given on page 279 for ,F,(1/2; 3/2; 1/h) 
in Jahnke and Emde,® which in their notation is 


8E. Jahnke and F. Emde, Tables of Functions (Dover 
Publications, New York, 1943). 
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denoted as M(1/2; 3/2; 1/h). Since 
1/Vh 
1Fu(1/2; 3/2; 1/8) = f edt, (18) 
0 


the table given on page 32 of Jahnke and Emde for 
the exponential integral can be used since it is more 
convenient for numerical work. For the examples to 
be considered, the exponential term in (17) can be 
completely neglected since, to justify the constant 
field assumption, Vio/Vr will always be greater 
than 10. Consequently, 


G,= (41 oRT./107%eh Vio?) { (h+-2)e71/* 
iF (1/2; 3/2; 1/h)—h} mhos/cm?, (19) 


where h is given by (11). For a very low frequency 
Vi1 coswt (@~0), the first term in the bracket ex- 
pression in (19) is 4/3 (see Appendix B), so 

G, = (32]omw*d?k?T,,”) /(3e* Voot1 0?!) 


w-0 


= (32m wd? Vr?) /(3e107 Viot) mhos/cm?. (20) 
The frequency square dependence of the total emis- 
sion conductance has been experimentally obtained 
by van der Ziel. For an extremely high frequency 
Vs1 coswt (see Appendix C), 


G, = (2e107Jo)/(mw*d?) mhos/cm?. (21) 


Ill. EXAMPLES 


(1) To illustrate (17), the constants given by 
Kuper* for the frequency and diode spacing will be 
used. f=3X10° c.p.s.; d=8.89X10- cm; T,=10' 
deg.; Jo =0.1 amps./cm*. Since Voo/ Vr = 11.605 Vio 
for the above given value of the cathode tempera- 


ture, the exponential term in (17) can be neglected * 


and (19) can be used for the calculation. To justify 
the assumed field distribution, the stationary po- 
tential distribution was calculated for the above 
given values of d, T., and J». The results are shown 
in Fig. 2. For retarding potentials greater in magni- 
tude than 1 or 2 volts, the assumption of a linear 
potential distribution is certainly justified. In Fig. 3, 
the total emission conductance and its reciprocal, 
the resistance, is plotted as a function of the re- 
tarding voltage Vso. The rapid increase in loading as 
the retarding voltage is decreased is clearly ex- 
hibited. The resistance curve is in good agreement 
with that given by Kuper,’ who cites the experi- 
mental results obtained at the M.I.T. Radiation 
Laboratory by Dr. H. V. Neher, when the shunt 
resistance of the cavity used in Dr. Neher’s experi- 
ment is considered. Since the stationary anode cur- 


rent is related by the Boltzmann expression 
I,=Io exp(— | Voo|/ Vr) (22) 
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‘ 


to the emission current, the value of the anode cur- 
rent for the range of Vso shown in Figs. 3 and 4 is of 
the order of 10-* amps./cm? or less. Consequently, 
the electrons that cross to the anode will have a 
completely negligible effect on the conductance. 

(2) Another interesting example is the variation 
of loading with frequency for a constant value of 
retarding voltage. The values of d, T., and Jo used, 
were those of example 1. Figure 4 shows the varia- 
tion of the total emission conduction as a function of 
frequency for retarding voltages of 3, 5, and 7 volts. 
Kuper® states that the loading should decrease as 
the frequency increases, though experimental re- 
sults seem to indicate the opposite. As Fig. 4 shows, 


. the behavior at 3X 10° c.p.s. is an increase in loading 


as the frequency is increased. The decay of the 
loading for very high frequencies is also very inter- 
esting and can be explained as follows. For a very 
high frequency V5; coswt, and consequently large 
transit angles, the energy interchange between the 
field and electrons decreases (note Eq. (8)). In addi- 
tion, a small change in the emission velocity will 
produce a large change in the transit angle and a 
greater averaging out among the different electron 
velocity groups. The asymptotic value of the con- 
ductance is given by (21). In the derivation of (21) 
the assumptions used in Appendix A must be noted. 
Namely, though (21) appears to be independent of 
Vio, it was assumed that Vio>V 7 so that the first 
term in (17) could be neglected. 

(3) Another interesting calculation is on the dis- 
tance the electrons, that produce the maximum 
loading, penetrate in the retarding field. This can be 
obtained by differentiating (13) with respect to @ 
and setting the resultant equal to zero. The follow- 


o 






“a 


Vo/ts 


Distance 1-107 cm 


Fic. 2. Stationary potential distribution for different 
values of anode retarding voltage. 
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ing transcendental equation is obtained for the 
halfway transit angle (@) for the electrons pro- 
ducing maximum loading, 


cos26y{1+2047(1—h)} 
== | + 2hO4?(O14 sin26 4 —_ 1). (23) 


Using the data given in example (1), 84 was calcu- 
lated and is shown in Fig. 5 plotted against the 
retarding voltage Vo. In addition, the ratio of the 
maximum distance (xy) the electrons penetrate in 
the field, to the diode spacing d, is plotted. 


CONCLUSIONS 


The large magnitude of loading calculated above 
suggests that, for close spaced high frequency tubes, 
the cathode-first grid region should not be in the 
radio frequency circuit. Though the above calcula- 
tion was made assuming a linear potential distribu- 
tion, while a triode in normal operation would have 
a potential distribution different from linear outside 
the cathode (see Fig. 2), the order of magnitude of 
the loading obtained will be very nearly the same. 
This can be clearly seen from the behavior of the 
experimental curves.2* An obvious, though me- 
chanically difficult, way of excluding the potential 
minimum region from the radio frequency circuit is 
to operate a tetrode as a triode with the first and 
second grids as the input cavity region. By proper 
spacing the potential minimum can be made to 
reside in the cathode-first grid region, thus obtaining 
the benefit of shot noise reduction by space charge 
without undesirable loading. 
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APPENDIX A 


By repeated integration by parts, the integral in (13) be- 
comes 


6 
Ps e*(1 —cos20—6 sind) ad0 
=e Lin 1) cos26;+h(6; sin2@;—1)] 


—1+(h+2) sindoc-*ao (a) 


Because of the rapid decay of the exponential in the sin26, 
integral, the upper limit can be replaced by infinity. The 
maximum possible error introduced by this change of the 
upper limit can be evaluated by setting sin2@ equal to one 
and integrating from 6; to infinity. The value of this integral 
is given in terms of the error function as 


cruel} 
Sie de i(: 1 —erf V;, (b) 


The exp(—6?) sin2@ integral with infinity as the upper limit 
can be evaluated in the following manner. Let 


u =f. e~** sin2pxdx. (c) 
Differentiating with respect to the parameter p gives 
du/dp=2 [ xe-** cos2prdx. (d) 


Using integration by parts, the integral can be written in 
terms of u as 
du 1 
—=—(1—2pu). e 
dp hs p (e) 


The solution of this first order differential equation is 


unter & ff er*rap, (f) 


where the constant of integration is zero since u=0, when 
p=0. Changing the integration variable in (f) and setting 


6&0? mhos-em* 





) I 70 Ss 20 2s 
Frequency 0:10 cps 


Fic. 4. G, versus frequency of Vo. 
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p=1, give 


1 1nd 
—h2z? <3 enum le 12 . 
I, e sin2xdx ae fi edt. (g) 


Using the expansion for e and integrating term by term, 


. 2 T'(2m+1) (1\" 
t — 
Wf, dtm 2 aor FIN th) 





By the use of Legendre’s duplication formula® and the defini- 
tion of the confluent hypergeometric function,” 


3 
wf ‘edt=sFi(3; 53 1). (i) 


Substituting (i) in (g) and the resulting expression in (a) 
gives (17). 


APPENDIX B 


To obtain the small transit angle expression for (19), the 
iF, and e~* can be written in series form and multiplied 
together. However, an easier method is to make use of 
Kummer'’s first transformation.!! 


’2f 3 1 : 
* a Gt 2) = r(155; -;). () 


Using the expansion shown in (h) 
3 1) = 2 4 8 16 
(1; 2) = —3R Tse ~ 105K" * 945k 5. 


°E. T. Copson, An Introduction to the Theory of Functions 
of a Complex Variable (Oxford University Press, New York, 
1935), p. 225. 

10 See reference 9, p. 260. 

1 See reference 9, p. 261. . 
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Substituting (k) in the bracket expression in (19) 
131 } 
ih 
{ (h+-2)e- iFi\ = 22°} —h 


4 12 32 
={1-540-se 





+ +f (I) 
APPENDIX C 
Using the asymptotic expansion for ,Fi(1/2; 3/2; 1/h) for 
large (1/h) gives 


1 3 1 athe += sit hb } 


FNS 55 7, (m) 


Substituting (m) in the bracket expression in (19), 


{ (h-+2)¢~ ash A383 : :t)- nben|1+h+2 h?+-- |. 
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A Visual Method for Demonstrating the Path of Ultrasonic 
Waves through Thin Plates of Material 


CHARLES J. BurToN* AND R. BowLinG BARNEsS* 
Stamford Research Laboratories, American Cyanamid Company, Stamford, Connecticut 


(Received December 1, 1948) 


A visual method for studying the path of ultrasonic waves through thin plates of material has been 
devised, and a set of photographs showing the complexity of ultrasonic reflection and transmission by 
metal plates has been obtained. An attempt has also been made to estimate the critical angles of total 
reflection for dilatation and shear waves by visual examination of the reflection of a supersonic beam. 
Using these data, Young’s modulus and the shear modulus for a ?#g-inch aluminum sheet have been 
determined and the values are compared with those obtained using a previously described electrical 


method and with published data. 





INTRODUCTION 


ULTRASONIC techniques for determining the 
mechanical properties of metals, plastics, and 
glasses have been the subject of considerable study 
during recent years. Details concerning various 
phases of this work have been described in a number 
of papers.'~* Schneider and Burton’ have shown 
that a rotating plate technique in which ultrasonic 
transmission is plotted electrically as a function of 
the angle of incidence of the waves allows determi- 
nation of the velocities of dilatation and shear waves 
in the plate. From these data, Poisson’s ratio and 
the mechanical moduli may be determined. Elastic 
constants of materials measured in this way are in 
good agreement with previously published values. 
Another paper® described an improved optical tech- 
nique for visualizing ultrasonic waves, the method 
having been developed in an effort to clarify some 
of the data which appeared in the earlier work. 

It is the purpose of the present paper to describe 
briefly a series of experiments showing the nature of 
the reflection and transmission of an ultrasonic beam 

















SCREEN 


Fic. 1. Schematic diagram of optical arrangement used for 
making sound waves visible. 


* Present Address: American Optical Company, Southbridge, 
Massachusetts. 
1L. Bergmann, Ultrasonics (John Wiley & Sons, Inc., New 
. York, 1938). 
27'W. P. Mason and T. G. Kinsley, Bell Lab. Record (May 
18, 1948). 
oui” Ballou and S. Silverman, J. Acous. Soc. Am. 16, 113 
1944). 
*D.L. Arenberg, J. Acous. Soc. Am. 20, 1 (1948). 
5A. W. Nolle, J. App. Phys. 19, 753 (1948). 
6S. Siegel, J. Acous. Soc. Am. 16, 26 (1944). 
049) C. Schneider and C. J. Burton, J. App. Phys. 20, 48 
1949). 
( —_ Bowling Barnes and C. J. Burton, J. App. Phys. 20, 286 
1949). 


462 


incident upon an aluminum plate at various angles 
of incidence. It will be shown that the critical angles 
of total reflection for longitudinal and shear waves 
can be rapidly determined by inspection of the 
optical pattern and that these data can be used to 
obtain with fair accuracy the mechanical moduli 
and Poisson’s ratio. Although not capable of high 
precision in its present stage of development, the 
method does have the advantage that the me- 
chanical properties may be determined very rapidly 
from a single set of measurements. 


EXPERIMENTAL 


The apparatus used in the present set of experi- 
ments is the same as that described in detail in a 
previous paper.*® It comprises an optical system of 
the schlieren type which makes possible the direct 
visualization of an acoustic wave in a liquid medium 
as a result of’ refractive index gradients which are 
set up by passage of the compressional wave 
through the medium. An X-cut quartz crystal is 
employed as the electromechanical transducer, pro- 
viding ultrasonic vibrations of 5-megacycle fre- 
quency, or at odd harmonics thereof. An electronic 
oscillator capable of supplying a high frequency po- 
tential of 300 volts is used to drive the quartz 
crystal. A device for supporting thin plates of vari- 











Fic. 2. Photograph of the experimental equipment for 
making sound waves visible: oscillator and point light source 
at left, ultrasonic tank in center and camera at right. 
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Fic. 3. Ultrasonic beam incident upon an aluminum plate at 0°. 





Fic. 4. Ultrasonic beam incident on an aluminum plate at 5°. 


ous materials at discrete angles with respect to the 
ultrasonic beam was constructed and used during 
the course of this investigation. A schematic dia- 
gram and a photograph of the equipment taken 
directly from the preceding paper® are included in 
Figs. 1 and 2, respectively. 

For precise visual measurement of angles and 
careful observation of the phenomena of transmis- 
sion and reflection of ultrasonic waves it is highly 
desirable to employ a narrow beam or “‘pencil’’ of 
ultrasonic energy. Since, however, the wave-length 
of 5-megacycle ultrasonic vibrations in water is 0.3 
mm, narrow beams of energy as employed in the 
present investigation are not very wide in terms of 
wave-lengths (ca. 25—30d) and as a result diffraction 
effects are pronounced. This situation presents some 
difficulty in visual ultrasonic studies, and there 
appears to be no simple solution to the problem. To 
a limited extent the problem can be alleviated by 
using as narrow a beam as possible since when the 
beam width is halved the diffraction effects are 
condensed into one quarter the distance from the 
radiator. In the present work, an effort has been 
made to use a small source (3?” diameter), to mask 
with a suitable set of acoustically absorbent baffles 
the diffracted lobes which appear and to work as far 
as possible from the source consistent with the 
absorbing properties of the tank liquid. Under these 
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conditions, diffraction effects are minimized and 
present no great problem in practical work. It 
should be noted in passing that diffraction effects of 
this nature are present in similar optical experiments 
although in this case the dimensions of the appa- 
ratus are so large compared with the wave-lengths 
of light that the effects become negligible. 

Another problem of some difficulty is the con- 
struction of a tank of small dimensions which will at 
the same time be anechoic. In the past, various 
workers have gone to great lengths in an effort to 
circumvent this problem: experiments have been 
carried out in lakes which approximate an “‘infinite”’ 
medium ; tanks of considerable size have been con- 
structed with double walls, the inner wall being 


‘constructed of some material (e.g., rho-c rubber) 


whose acoustic impedance closely matches that of 
water and with the space between the walls filled 
with some highly absorbing material. In the present 
study, no extreme measures of this sort have been 
taken. The tank has been lined with rubber and 
several layers of copper window screen have been 
placed on the bottom and ends of the tank to pre- 
vent the formation of strong reflections. Some 





Fic. 5. Ultrasonic beam incident on an aluminum plate at 10°. 





Fic. 6. Ultrasonic beam incident on an aluminum plate at 15°, 
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Fic. 7, Ultrasonic beam incident on an aluminum plate at 20°. 


spurious reflections still exist, particularly those re- 
sulting from water wave reflection, although these 
have been relatively easy to identify and have not 
interfered with study of the primary plate transmis- 
sion and reflection. 

A series of photographs showing the nature of the 
transmission and reflection of an ultrasonic beam by 
a js-inch aluminum plate for angles of incidence 
between 0° and 45° at intervals of 5° is included as 
Figs. 3-12. A careful study of these photographs 
brings out several interesting features. The initial 
transmission is relatively small, decreases somewhat 
at 5° angle, increases through 10° and 15°, increases 
evén more at 20°, decreases through 25° and be- 
comes zero when the angle increases beyond ap- 
proximately 30°. The reflections of the acoustic 
waves are also of interest: noticeable specular re- 
flection is evident at 5° and 10°, with a second 
reflected beam making its appearance at 10°. Calcu- 
lation indicates that the second reflection can be 
readily explained as a result of a partial reflection of 
the ultrasonic beam at the aluminum-water inter- 
face. The phenomenon is not unfamiliar in optics, 





Fic. 8, Ultrasonic beam incident on an aluminum plate at 25°. 
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Fic. 9. Ultrasonic beam incident on an aluminum plate at 30°. 


witness the mechanism of the Fabry-Perot inter- 
ferometer, the Lummer-Gehrcke plate, etc. Actually, 
at intermediate angles as many as 7 reflected and an 
equal number of transmitted beams can be seen and 
photographed. It is extremely interesting to watch 
the “parade” of interference bands moving up and 
down the face of the plate as it is moved slowly 
through these regions. One should also note that the 
refraction in the case of an ultrasonic beam passing 
from a liquid to a metal is away from the normal, as 
opposed to the similar case in optics where light 
under these conditions is bent toward the normal. 
This situation arises, of course, because of the 
greater velocity of sound in a metal than in a liquid. 
The multiple-beam reflection disappears after reach- 
ing an angle of approximately 30° and strictly 
specular reflection then obtains. 

As previously pointed out, these phenomena can 
be partially explained by analogy with familiar 
optical effects. On the other hand, many of the 
observed details do not appear to be explainable on 
this basis. This is understandable: light is a trans- 
verse wave phenomenon, whereas sound in fluids is 
propagated as a compressional wave and in a solid as 
either a compressional or shear wave, or both. It has 
been stated by other workers? in the field that as 
many as six different transmitted waves may result 
from the incidence of a single ultrasonic wave upon 
a thin plate of material, and there is considerable 
difficulty in separating observationally the longi- 
tudinal and transverse wave systems. Furthermore, 
there are at certain angles of incidence so-called ‘‘ex- 
change waves” which are produced by the trans- 
verse wave, generated at the farther wall by a longi- 
tudinal wave in the plate, interfering in phase, on 
meeting the front wall and being reflected again, 
with a second entering ray. The extreme complexity 
of this situation, in addition to the fact that certain 
of the theoretical points are still highly speculative, 
makes it virtually impossible to analyze the problem 
completely with any degree of certainty. The photo- 


*R. Bar and A. Walti, Helv. Phys. Acta 7, 658 (1934). 
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Fic. 10. Ultrasonic beam incident on an aluminum plate at 35°. 


graphs themselves, simple as they appear at first 
glance, are in reality quite complex and do not 
clarify the problem of ultrasonic transmission 
through plates to the extent which we had initially 
hoped. 

It is pertinent at this point to consider the 
mathematical relations which are required for calcu- 
lation of the mechanical moduli from ultrasonic 
transmission data. The indices of refraction for 
dilatation and shear waves are given by 


sind V, 
ng =—_ = — 
sin Ya Va 
sin@ V, 
n,=—_ => — 
sing, V, 


where 
na4=index of refraction for dilatation waves, 
n,= index of refraction for shear waves, 
V.=velocity of sound in the liquid, 
Va=velocity of dilatation waves, 
V.=velocity of shear waves. 


The angles of total reflection 6; and 62 are given 
by the preceding equation when gg and yg, become 
90° and their sines become 1. Under these conditions 








Vi 
Va= 
sin6, 
Vi 
V, = 
sin 65 


The velocity of propagation of the shear waves is 


(0) 


where G=modulus of rigidity, p=density of plate. 
The velocity of the dilatation waves is 


(1—y)E ‘ 
Va= 
e +u)(1 ore 
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_ Fic. 11. Ultrasonic beam incident on an aluminum plate at 40°. 


where E=Young’s modulus, p=Poisson’s ratio. 
Poisson’s ratio is expressed as 


k?—2 
.~=———— 
2(k?—1) 
where 
k= Va/ Vs. 


And the three constants are inter-related in the 
following way : 


E=2(1+y)G. 


The foregoing equations are of course strictly cor- 
rect only when applied to homogeneous, isotropic 
media of infinite extent. These conditions are obvi- 
ously not met in the case of a rolled sheet of metal, 
nor are they met in practically any experimental 
arrangement which can be set up. Since, however, it 
is impossible to specify the degree of inhomogeneity 
or the extent of non-isotropy, it appears legitimate 
to use these expressions as a first-order approxima- 
tion to the “true” values of the dynamic moduli. 

It has been shown! that characteristic transmis- 
sion patterns are obtained when the acoustic energy 
transmitted through a plate is plotted as a function 





Fic. 12. Ultrasonic beam incident on an aluminum plate at 45°. 
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of the angle of incidence of the energy. Figure 13 is 
an example of such a curve, showing the transmis- 
sion through a ;;-inch sheet of aluminum. One notes 
that at 12° and 30° there is approximately zero 
transmission of acoustic energy. This is interpreted 
to mean that total reflection of the longitudinal 
waves which are excited in the plate occurs at ap- 
proximately 12°. Beyond this point strong shear 
waves are set up in the plate, their intensity building 
up to a maximum and then receding to zero as the 
point of total shear reflection is reached. The 
periodic structure in both the longitudinal and shear 
portions of the curve is, according to one interpre- 
tation, a manifestation of the interference between 
the ‘‘exchange’”’ waves and the incident longitudinal 
waves. There appear to be conflicting opinions’® in 
this regard although reference to our earlier work’ 
will show that these interference maxima can be 
used to determine the indices of refraction for 
dilatation and shear waves. 

The angles of total reflection for the shear and 
dilatation waves can also be estimated with a fair 
degree of precision by inspection of the photographs 
in Figs. 14-16. It can be seen in Fig. 14 that there is 
zero transmission through the plate at this specific 
angle (11.5°). In Fig. 5, on the other hand, definite 
transmission is visible showing clearly the rapid 
building up of the shear wave intensity in the plate 
at an angle of 15°. From our previous discussion of 
the transmission curve (Fig. 13) it will be recalled 
that this corresponds to the region of total reflection 
for longitudinal waves. Thus, it is possible to esti- 
mate that this angle of total reflection for aluminum 
occurs at 11.5°. Similarly, the angle of total reflec- 
tion for the shear waves is that angle beyond which 
all transmission ceases. From Figs. 15 and 16 it can 
be seen that this effect occurs at approximately 31.5°. 

Using these data and substituting in the mathe- 
matical expressions recorded previously, the values 
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Fic. 13. Ultrasonic transmission pattern for 7-inch aluminum 
plate, plotted electrically. 








© J."Gotz, Akustische Zeits. 8, 145 (1943). 
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TABLE I.f 

Electrical* Visual Published** 

recording recording value 
Dilatational velocity (cm/sec.) Va 7.05105 7.0105 
Shear velocity (cm/sec.) V; 2.82105 2.9105 
Poisson's ratio 0.40 0.4 0.33 
Shear modulus (dynes/cm?) G 2.1410" 2.310" 2.3-2.7 X10” 
Young's modulus (dynes/cm?) E 6.73 X10” 64X10" 6.3-7.5 X10” 








+ Material—Aluminum sheet (¢ inch); specific gravity—2.70. 
* See reference 6. 
** Landolt-Boernstein tables. 


of the longitudinal and shear velocities, the me- 
chanical moduli and Poisson’s ratio are readily ob- 
tained. The calculated values are recorded in 
Table I and are compared with those previously 
obtained using the electrical method of recording 
ultrasonic transmission data and with published 
handbook values. There is good agreement with 
handbook data save in the case of Poisson’s 
ratio which is somewhat higher than the accepted 
value. However, in considering these data, it must 
be recognized that the calculations of mechanical 
properties have been made merely to show that rapid 
measurement of these properties is possible with 
simple equipment and with reasonable accuracy. 
Moreover, virtually all handbook data on me- 
chanical properties are a result of static determina- 
tions, and dynamic measurements of the type shown 
here cannot be expected to be in perfect agreement 
with these values. 

In conclusion, it is felt that the foregoing experi- 
ments are instructive for illustrating the complexity 
of the phenomena of ultrasonic plate transmission 
and reflection. In addition, the technique may be 
used to determine rapidly and with reasonable pre- 
cision (+1°) the angles of tctal reflection for dila- 
tation and shear waves. These values in turn may be 
used to calculate the dynamical values of Young’s 
and the shear modulus with a precision limited 
mainly by the determination of the dilatational 
wave velocity (ca. +8 percent). It is of course obvi- 





Fic. 14. Ultrasonic beam incident on an aluminum plate at 
11.5°, just at the critical angle for total reflection of dilatation 
waves. Note relatively high transmission in Figs. 4 and 5, for 
10° and 15° respectively. 
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Fic. 15. Ultrasonic beam incident on an aluminum plate 


at 29.5°, just before reaching the angle of total reflection for 
shear waves. 


ous that the precision of the method as a technique 
for measuring mechanical moduli can be improved 
to almost any desired extent by devising either 
electrical or mechanical methods for improving the 
precision of angle measurement. However, these 
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Fic. 16. Ultrasonic beam incident on an aluminum plate at 
31.5°, just at the critical angle for shear wave reflections. 


experiments have been carried out primarily to ob- 


serve ultrasonic reflection phenomena by metal 
plates and the calculations of mechanical moduli 
have been made merely to illustrate the fact that 
the angles of total reflection for dilatational and 


shear waves can be observed easily and measured 
rapidly. 
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The High Current Carbon Arc and Its Mechanism 


WOLFGANG FINKELNBURG 
Engineer Research and Development Laboratories, Fort Belvoir, Virginia 


(Received October 26, 1948) 


The high current carbon arc with maximum crater brightness of 2000 candles per mm, axial in- 
tensity of 1.2 10* candles, power inputs up to 200 kw, light efficiencies exceeding 90 lumens per watt, 
radiation efficiency of 73 percent, crater temperatures up to 8000°K, and arc stream temperatures of 
the order of 12,000°K is the most powerful and one of the most interesting radiation sources known. 
The rising voltage characteristic, the anodic vapor stream which causes its excellent radiation proper- 
ties, and the contracted arc stream distinguish the high current carbon arc from the well-known 
normal low current carbon arc. The properties of the contracted arc stream are discussed and a 
general theory of arc streams is presented. Starting from the anodic mechanism for the low current 
carbon arc, the anodic mechanism of the high current arc is developed. The rising voltage charac- 
teristic, the anodic vapor stream, and the high crater brightness are explained as a consequence of an 
abnormal voltage drop in front of the anode which increases with increased current density and is 
caused by a very rapid evaporation of the anodic material. The important role which the magnetic 
field of the arc current plays in the stabilization of all high current arcs is pointed out. The author 
believes that this unique high temperature arc has bearing on other arc and spark discharges, and also 


will play an important role in future developments of high temperature physics and chemistry. 





INTRODUCTION 


HE high current carbon arc, especially in the 
most common form of the high intensity arc 
or Beck arc, has, for more than three decades, pro- 
vided the most powerful and brilliant light source in 
existence.! Only recently is this unique high temper- 
ature arc discharge finding a more general interest 
among physicists and chemists. In the first place, it 
is the only easily accessible source where the be- 
havior of matter and its reactions can be studied at 
temperatures between 6000 and 12,000°K (plasma 
physics, high temperature chemistry, etc.). In the 
second place, this arc displays certain features 
unique in gas discharge physics and hence an under- 
standing of its mechanism helps our general knowl- 
edge of the more complicated arc and spark 
processes. In this connection, the high current car- 
bon arc exhibits (a) the first example of a discharge 
whose properties are essentially determined by its 
own magnetic field; (b) an arc stream where like in 
a metal the electric and heat conductivity are due to 
free electrons; (c) the first well studied case of a 
discharge whose mechanism is fundamentally in- 
fluenced by mechanical vapor streams which are 
produced by the discharge itself. All these phe- 
nomena seem to play an important role in spark 
discharges.” 


THE HIGH CURRENT CARBON ARC AS A NEW CLASS 
‘ OF ARC DISCHARGE 


Figure 1 is a schematic sketch and Fig. 2 a photo- 
graph of a fully developed high current carbon arc. 


1 Wolfgang Finkelnburg, Der Hochstromkohlebogen (Verlag 
Julius Springer, Heidelberg, 1948). English presentation (not 
— as complete) : The =~ Current Carbon Arc (FIAT Final 

eport 1052—PB-81644, Dept. of Commerce, Washington, 
D. C.,- 1947). 

2 Wolfgang Finkelnburg, Phys. Rev. 74, 222 and 1475 (1948), 
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Both figures display two important distinctions of 
the high current carbon arc from the normal low 
current carbon arc: (a) the deeply molded positive 
crater is filled with a streaming, brilliant, high tem- 
perature vapor ejected from the base of the crater 
which accounts for the high crater brightness, and 
(b) for currents above 100 amperes the arc stream 
appears to be contracted into a slender flamelike 
column (Fig. 3) with a current density more than 
ten times that of the low current arc stream. Further 
distinctive properties will be discussed later. 

Ever since its discovery by H. Beck 35 years ago, 
the high current carbon arc has been operated with 
a cored positive carbon containing primarily cerium 
fluoride and oxide at anodic current densities of 100 
to 400 amp./cm?. This exceeds the current density 
of the usual low current carbon arc by a factor of 
3 to 10. 

It was generally assumed that this cerium com- 
pound core material was a necessary condition for 
the operation of the high intensity arc. The author 
found that an essentially similar phenomena could 
be obtained with many kinds of positive electrodes 
if operated at a sufficiently high anodic current 
density. Indeed, we now believe the high current 
carbon arc should be recognized as a new form of an 
arc discharge which appears whenever sufficiently 
high current density causes vapor eruptions from 
the anode and, in certain cases, also from the 
cathode. The term high current arc has been intro- 
duced to distinguish this general type of discharge 
from the included special high intensity arc using 
cerium compounds in the core of the positive carbon. 
The vapor eruption from the positive carbon influ- 
ences the arc mechanism in a way to be described 
later. It should be mentioned that the general shape 


’ Wolfgang Finkelnburg, Zeits. f. Physik 112, 305 (1939). 
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of the arc, especially the direction of the anodic 
vapor flame and the arc stream as shown in Figs. 1 
and 2, is independent of its orientation in space and 
does not result from the convection of hot vapors or 
gases. It is, rather, a general property of all high 
current arcs of more than 100 amperes that their 
specific shape is greatly influenced by their own 
magnetic fields.* We call them, therefore, field- 
stabilized in contrast to the wall-, electrode-, or 
convection-stabilized arcs.‘ 


THE ANODIC PROPERTIES OF THE HIGH CURRENT 
CARBON ARC 


Before we begin with the discussion of the dis- 
charge mechanism, it seems convenient to list the 
most important properties of high current carbon 
arcs, specifically those of the high intensity arc 
insofar as radiation properties are concerned. As has 
been mentioned, the high current carbon arc is 
operated with anodic current densities between 100 
and 400 amp./cm? compared with approximately 40 
amp./cm? for the low current arc. Evidence for a 
fundamental difference in the arc mechanism is the 
rising voltage characteristic of the high current arc, 
which is caused by an increase of the anode voltage 
drop with increased current. Partly because of this 
larger current density and partly because of the larger 
anode drop, there results a comparatively large con- 
sumption of the positive carbon ranging from 10 to 
about 50 inches per hour compared with the average 
of 3 to 4 inches per hour for the low current arc. The 
author was able to prove that the consumption in 
the high current carbon arc is due to evaporation, 
whereas in the low current arc it is predominantly 
due to chemical combustion.’ This rapid evapora- 
tion of the positive carbon affects the core more than 
the harder carbon shell and thus causes the deeply 
molded crater. As a consequence of this evaporation 
and an additional heating effect, to be discussed 
later, the positive crater of the high current carbon 
arc in contrast to the low current arc is filled with 
brilliantly radiating vapors. These vapors, ejected 
from the carbon, form the anodic vapor ‘‘flame”’ 
mentioned above and indicated in Figs. 1 and 2. 
Whereas in the low current carbon arc only the 
incandescent electrode surface at a temperature of 
4000°K or below is responsible for the radiation, in 


* At currents above 500 amperes the magnetic field often 
Causes a spiral curvature of the contracted arc stream and a 
very fast rotation of it round its foot point. This motion dis- 
turbs the arc considerably and prevents its steady operation. 
British searchlight investigators, not knowing what caused 
this unsteadiness, designated this arc state as the “exploded 
arc” (Symposium on Searchlights, Illuminating Eng. Soc., 
London 1948). All details of this phenomena, its theory, and 
methods of its prevention may be found on pages 11 and 127 of 
The High Current Carbon Arc (reference 1). 

* Wolfgang Finkelnburg.and K. H. Hoecker, Naturwiss. 33, 
55 (1946). 

* Wolfgang Finkelnburg, The High Current Carbon Arc 
(reference 1), p. 77. 
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the high current arc the crater vapors at tempera- 
tures of 6000 to 8000°K are the main source of 
radiation. If these anodic vapors arise from com- 
pounds containing elements like cerium, which has 
closely spaced spectral lines throughout the visi- 
ble region, copious amounts of white light are 
emitted. Crater brilliancies of arcs in practical 
application range from 500 to 1600 candles per mm? 
and brilliancies exceeding 2000 candles/mm? have 
been measured with experimental arcs. This bril- 
liancy exceeds that of the surface of the sun (1600 
candles/mm?). A total radiation as high as 6 kw per 
cm? crater cross section has been obtained. This 
total radiation is equal to that of the sun’s surface 


_ and that of a blackbody at 5800°K. In further con- 


trast to the low current carbon arc, the light effi- 
ciency is amazingly high. More than 70 percent of 
the total power input of high intensity arcs is 
emitted as radiation, and the light efficiencies (in- 
cluding the radiation of the anodic vapor flame) 
range from 55 to 90 lumens per watt compared with 
23 lumens per watt for one of the best low current 
carbon arcs. Table I presents a summary of inter- 
esting data for a number of representative high cur- 
rent arcs in comparison with a low current arc.® 


THE CONTRACTED HIGH CURRENT ARC STREAM 
AND ITS THEORY 


Although all the properties of the high current 
carbon arc discussed so far are essentially anodic 
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_ Gantrected Arc Stream 


Fic. 1. Schematic sketch of the high current carbon arc with 
the different regions of theoretical importance. 


6 Wolfgang Finkelnburg, Optik 2, 149 (1947). 
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Fic. 2. Photograph of a high current carbon arc, burning 
with 185 amperes between water-cooled carbons (positive 
carbon 11 mm, negative 9 mm diameter). 


phenomena, there exists a property common to all 
stationary arcs operated at sufficiently high current. 
At currents above 80 amperes the normal diffuse 
and expanded arc stream seems to contract into the 
above-mentioned bluish-white narrow column which 
we called contracted high current arc stream (Fig. 3), 
and which is surrounded by well defined zones of 
different color. Our investigations’ indicate that 
practically all the current passes through the con- 
tracted column at the high current density of 1000 
to°3000 amp./cm?*. It causes there a temperature in 
excess of 10,000°K which manifests itself by a strong 
continuous spectrum with superimposed widely 
broadened spectral lines of atomic nitrogen, oxygen, 
and carbon, but without any band spectra. 

’ The theory of the contracted arc stream has been 
developed recently by Hoecker and the author® and 
corrected in details by Hoecker and Schulz.’ In 
an arc stream of given diameter, according to 





7 Wolfgang Finkelnburg, Zeits. f. Physik 116, 214 (1940). 

®K. H. Hoecker and W. Finkelnburg, Zeits. f. Naturfor- 
schung. 1, 305 and 382 (1946). 

*K, H. Hoecker and P. Schulz, Zeits. f. Naturforschung. 
(1948), in press. 
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Elenbaas,'® the electric field strength E adjusts 
itself to such a value that the power input per cm 
arc length JE is compensated by the radial heat 
conduction and by radiation, S: 


IE=2rek(dT/dr)+S, 


where J is the current and & the coefficient of heat 
conduction of the arc stream plasma. In a molecular 
plasma like that of air, there are four processes 
contributing to the heat conductivity: the classical 
energy transfer by diffusing atoms (or molecules) ka, 


_ that by free electrons k,, that by molecules dis- 


sociating in the arc stream and recombining outside 
of it, thus conducting away the dissociation energy 
(ka), and finally that by atoms ionized in the arc 
stream and carrying away the ionization energy. 
The contributions of each of these four factors have 
been computed for air as a function of the tempera- 
ture by Hoecker.* His most interesting results are 
that below 7000°K the contribution of the free 
electrons, k,., as well as that of the ionization, R;, can 
be neglected, and in the region of 6800°K there 
exists a minimum of heat conductivity because the 
dissociation of the air is already complete and the 
ionization does not yet play an essential role. It is 
interesting that the measured axial temperatures of 
the low current arc stream in air agree fairly well 
with this theoretically computed temperature of 
6800°K. At temperatures above 10,000°K, on the 
other hand, according to Hoecker’s computation, 
the heat conductivity of the free electrons accounts 
for nearly all the heat conduction ; all other processes 
can be neglected in comparison. ‘The computation, 
carried through for a 200-ampere arc stream, re- 
sulted in an axial temperature of 11,000°K. Still 
higher temperatures are to be expected for higher 
currents. Consequently, the distinction between 
high and low current arc streams, introduced origi- 
nally by the author for phenomenological reasons,’ 
has found its clear physical meaning: The high cur- 
rent arc stream is distinguished from the low 
current arc stream by such a high a degree of 
ionization that the free electrons are predominantly 
responsible for the heat conductivity as well as for 
the optical behavior. The continuous spectrum re- 
sults from the stopping of the free electrons in the 
electric field of the positive ions, and the few 
emitted lines of N, O, and C are widely broadened 
as a result of the electric fields of the surrounding 
electrons and ions. 

Also the transition from low to high current arc 
stream seems understandable now. With increased 
current, the arc stream can either increase its diame- 
ter at the constant temperature of 6800°K, corre- 
sponding to the minimum of the total heat con- 
ductivity, or it can remain at constant diameter 


10 W. Elenbaas, Physica 1, 673 (1934). 
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with increased temperature and consequently with 
increased, electric conductivity. Above currents of 
approximately 100 amperes the second possibility 
seems to be more favorable for the energy balance 
because of the exponential increase of the electric 
conductivity with temperature. 


THE ANODIC MECHANISM OF THE LOW CURRENT 
CARBON ARC 


In order to understand the anodic mechanism of 
the high current carbon arc, we have to begin by 
discussing that of the low current arc. We will 
discuss several limiting cases and assume at first 
that the positive carbon does not emit positive ions 


even at its high temperature of the order of 4000°K. . 


Immediately in front of the positive electrode the 
total current then is a pure electron current. The 
resulting negative space charge causes a potential 
drop immediately in front of the positive electrode 
surface. It is called the anode drop. By it, the 
electrons are accelerated and, by ionizing collisions 
with atmospheric gas molecules, produce the posi- 
tive ions which diffuse in the electric field toward the 
cathode and which are necessary in order to compen- 
sate the negative space charge in the quasi-neutral 
arc stream. 

Let it now be assumed that the positive carbon 
not only is able to emit positive ions at high temper- 
ature, as is the case to a certain extent for carbons 
containing metal compounds, but that its ion 
emissivity at 4000°K is sufficiently high to provide 
all the positive ions necessary for the quasineutral 
arc stream. Then no anode drop would be necessary 
for the production of ions, and we would find only a 
very small anode drop caused by the higher resist- 
ance of the ‘‘cooler’”’ gas layer of 4000°K near the 
anode surface compared with the normal arc stream 
with a temperature of about 6000°K. But this state 
of the discharge is not a stable one. The incident 
energy of the bombarding electrons would not com- 
pensate for the cooling of the anode by radiation 
and heat conduction and hence its ion emission 
would fall below the necessary amount. But as this 
occurs the excess of electrons would cause an in- 
creasing space charge in front of the anode surface 
and thus increase the anode drop until a state of 
equilibrium is reached. Hence it follows that a cer- 
tain part of the positive ions necessary is produced 
by ionization in the anode drop region while the rest 
is emitted by the anode surface, which is kept at 
sufficiently high temperature by. the incident elec- 
trons accelerated by the anode drop. Thus a certain 
anode drop is a necessary consequence of the arc 
mechanism. The anode drop of the low current 
carbon arc using an ion-emitting, metal-cored posi- 
tive carbon actually is only about 10 volts, whereas 
the anode drop for a pure homogeneous positive 
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carbon without ion emission is approximately 30 
volts. 


THE ANODIC MECHANISM OF THE HIGH CURRENT 
CARBON ARC 


Our interpretation of the anodic mechanism of 
the high current carbon arc, as distinguished from 
that of the low current arc, is based on one premise: 
the anodic vapor stream blows away the positive 
ions from the anode and thus, by an increased nega- 
tive space charge, causes the increase of the voltage 
drop in front of the anode with all its consequences, 
and it is these consequences which are characteristic 
of the high current carbon arc. 

To prove this premise, let us assume a low current 
carbon arc with the anodic mechanism described in 
the previous section. Now imagine a vapor stream 
of velocity u ejected from the anode surface through 
the anode drop into the arc stream. This velocity u 
has been found to be of the order of 10* cm/sec. and 
is vectorially added to the previously existing 
velocity of the positive ions and electrons in the arc 
stream. Since the velocity u is roughly 100 times 
that of the velocity of the positive ions in the arc 
stream due to the electric field, but only of the order 
of a few percent of the corresponding field velocity 
of the electrons, it is clearly evident that the positive 
ions are swept away from the region of the anode 
with a velocity essentially equal to uw. A corre- 
spondingly larger number of ions, therefore, must be 
produced by ionization in front of the anode to 
compensate for those carried away by the vapor 
stream. This increased ionization requires an in- 
creased anodic voltage drop. How this is established, 
may readily be seen: While the positive ions are 
swept away out of the anode drop region with a 
greatly increased velocity, the electrons still travel 
out of the arc stream into the anode drop region and 
to the anode with a velocity which is not changed 
considerably by the vapor stream. Consequently, 
the negative excess space charge in the anode drop 
region, and with it the anode drop voltage, increases 
until sufficient ionization by the electrons, acceler- 





Fic. 3. Photograph of a contracted arc stream at 180 
amperes, with positive carbon of so large diameter that current 
density is below the critical point; therefore no positive flame. 
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ated to greater velocity in the larger anode drop, 
establishes an equilibrium. In this stationary state 
not only enough ionization takes place in the anode 
drop region to compensate for the ions carried away 
by the vapor stream, but the electrons accelerated 
by the increased anode drop also carry to the anode 
enough energy in order to sustain by rapid vaporiza- 
tion the anodic vapor stream which is the dominant 
feature of the high current carbon arc and its anodic 
mechanism. 

We can now explain the transition from the low 
current arc to the high current carbon arc as follows: 
By increasing the anodic current density we transfer 
more energy to the positive crater by the electrons 
accelerated in the normal anode drop than can be 
dissipated by heat radiation and conduction. Conse- 
quently, the anode surface begins to evaporate. As 
each vapor layer is pushed away normal to the 
anode surface by the new vapor being formed, as if 
ejected from a nozzle in the anode surface, the 
evaporation results in a vapor stream. The velocity 
of this stream is given by the amount of electrode 
material evaporating per second multiplied by the 
increase of the volume due to vaporization of the 
solid material and heating it to the temperature of 
6000 to 8000°K."' The vapor stream causes an in- 
crease of the anodic voltage drop and thus of the 
electron energy transferred to the anode surface, so 
that in turn the evaporation rate increases until a 
state of equilibrium is réached. The energy trans- 


ferred to the anode surface by the electrons is given 
by 


W,=1( V.t+¢), 


4=electron current 
V,=anode drop 
¢=work function 


where 


and is equal to the energy dissipated by heat con- 
duction and radiation from the anode plus the 
amount consumed for the production of the high 
temperature vapor stream. The rising voltage char- 
acteristic of the high current arc in contrast to the 
falling characteristic of the low current arc is now 
explainable by our mechanism: an increase of the 
arc current causes an increase of the vapor stream 
velocity, and does so through increased vaporiza- 
tion. This in turn causes the increase of the voltage 
drop in front of the anode which has been measured 
by the author’ and causes the rising voltage charac- 
teristic. 

This anodic vapor stream is affected by a mag- 
netic field because near the anode it forms a part of 
the arc stream (in Fig. 1 designated as the ‘‘turbu- 
lent arc stream’’) in which the electrons are mainly 
carrying the arc current. Any magnetic deflection of 
the electron stream in this region will be imparted to 
the vapor stream and cause an over-all change of its 
direction because the anode flame is a continuation 
of the vapor stream from the region of the turbulent 


TABLE I. Characteristic data of high intensity arcs (see reference 6). 








Low current 





carbon arc High intensity arcs in practical application High load experimental arcs 

Diameter of 
positive carbon mm 10 7 7 9 16 24 20X30 7 10 13 
Current amp. 30 50 75 75 200 450 1200 75 200 500 
Arc voltage _—-volts 55 40 55 50 78 90 100 65 85 ~100 
Arc wattage | kw 1.65 2 4.1 3.75 15.6 40 120 4.9 17 ~S50 
Consumption mm/h 
of itive 
pet ove 100 320 1000 400 360 800 2000 2000 3000 4600 
Brightness 

candles/mm? 188 550 1000 700 850 1100 1600 1400 1700 2000 
Candle- 
power candles 9500 22,000 45,000 35,000 140,000 400,000 1,200,000 — — — 
Total light 
output lumens 37,000 110,000 280,000 225,000 ~900,000 ~3,400,000 ~10,000,000 — — — 
Light lumens per 
efficiency watt 23 55 68 60 ~60 ~85 ~90 — — _ 
Blackbody 
temperature 7, °K 
of pos. crater 3800° 4300° 5000° 4800° 4500° 4700° 4800° 5500° 5400° 5800° 











” Wolfgang Finkelnburg, Phys. Rev. 74, 1475 (1948). 
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arc stream. In the absence of external magnetic 
fields the anodic vapor stream is deflected in the 
direction of the bisector of the obtuse angle formed 
by the positive and negative carbons as shown in 
Figs. 1 and 2, by the magnetic field of the arc itself. 
This field, in trying to expand the arc current path, 
exerts an outwardly directed force on the arc 
stream, thus accounting for the curvature of the 
turbulent as well as the contracted arc stream as 
visible in Figs. 1 and 2. This outward curvature of 
both parts of the arc stream is greatest at their point 
of contact because here the inherent stiffness of both 
so-called ‘‘flames”’ is smallest, as has been discussed 
in detail elsewhere.''* We believe that the magnetic 
field of the arc also helps in the contraction of the 
high current arc stream which may be regarded as 
consisting of parallel threads conducting the current 
in the same direction and therefore attracting each 
other. Because at currents greater than approxi- 
mately 100 amp. these effects of the magnetic field 
predominate, as compared with other factors like 
convection currents of hot gases which influence any 
arc, Hoecker and the author‘ have designated high 
current arcs quite generally as ‘‘field-stabilized.”’ 


THE MECHANISM OF LIGHT EMISSION AND ITS 
DEPENDENCE ON THE ARC PARAMETERS 


The influence of the arc parameters on the light 
emission by the high current carbon arc crater 
follows from the anodic mechanism discussed above. 
We have two ways of regarding the problem, though 
both in effect are equivalent. We may rely on the 
fact that in most arcs at higher pressures we have 
thermodynamic equilibrium throughout the arc 
except for a space extending a few electron mean- 
free-paths from each electrode. The emission of the 
crater vapors is then considered as thermal emission 
and the effect of the arc parameters on the vapor 
temperature in the crater must be determined. On 
the other hand, we may disregard the temperature 
concept and assume that excitation of the radiating 
crater vapors is entirely caused by electron col- 
lisions. The radiation depends in both cases on the 
density, average velocity, and collision frequency of 
the electrons, which in the first model are mainly 
responsible for the energy transfer and the heating 
of the crater vapors, while in the second model they 
directly excite the observed line emission of these 
vapors. We clearly distinguish two different regions 
of energy exchange (Fig. 1), the anode drop region, 
confined according to our calculations to about a 
few tenths of a millimeter, and the vapors in and in 
front of the crater through which the current passes. 
These vapors have been called the turbulent arc 
stream because of their rapid internal motion in 
contrast to the adjacent contracted arc stream in air 
(Figs. 1-3). 

In the small anode drop region of a practical high 
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current arc, the high power of approximately 6 kw is 
dissipated per cm? of crater surface. This is sufficient 
to transfer to each atom of the vapor stream an 
average energy of 8 ev, i.e., more than the ionization 
energy of the vapor atoms. The anode drop space is 
thus the main region where excitation and heating 
to a very high temperature take place. Conse- 
quently, the higher the anode drop for a given 
carbon the higher the crater brilliancy. High anode 
drop, according to our discussion of the anodic 
mechanism, is caused, however, by the high velocity 
of the anode vapor stream resulting from the very 
rapid evaporation of the positive carbon. Our theory 
then accounts for the empirical relation that high 
crater brilliancy requires high anode drop (i.e., high 


‘arc voltage) and high consumption of the positive 


carbons. 

In addition to this energy exchange in the anode 
drop space, there is heating and excitation due to 
electron collisions in the turbulent arc stream within 
and immediately in front of the positive crater. As 
much as 2 to 3 kw is dissipated here per cm? cross 
section of the vapor stream. The influence of this arc 
region on the crater brilliancy, as measured from a 
point on the axis of the positive carbon, evidently 
depends on the current density in this region and on 
its thickness. This last point explains why the 
crater brilliancy depends on the crater depth. The 
same change, however, can be attained by in- 
fluencing the thickness of the vapor layer in front 
of the crater, for instance by magnetic means.” 
There is an interesting possibility of distinguishing 
between the influence of the anode drop and that of 
the current density on the resulting crater brilliancy. 
Normally, by increasing the current the anode drop 
is increased together with the current density. We 
can influence the anode drop alone, however, either 
by changing the chemical composition of the posi- 
tive carbon or, according to a recent investigation of 
the author," by cooling the positive carbon im- 
mediately behind the burning end. By comparing 
carbons of different composition, or a cooled one 
with an uncooled one, at the same current density, 
we thus can study the influence of the anode drop on 
the radiation characteristics. The influence of the 
current density can be studied by adjusting the 
current of arcs with different carbon composition 
(or with a cooled and an uncooled positive carbon) 
so that they show the same anode drop. However, 
no systematic studies of the last sort have been 
made. 

The effect of cooling the positive carbon, which 
has proved to be of great technical importance," is 

2 Wolfgang Finkelnburg, The High Current Carbon Arc 
(reference 1), p. 88. 

_ ® Wolfgang Finkelnburg, J. Soc. Mot. Pict. Eng. 52, (1949), 
in press. 


“ M. T. Jones and F. T. Bowditch, J. Soc. Mot. Pict. Eng. 
52, (1949), in press. 
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of considerable theoretical interest. The same posi- 
tive carbon, when cooled up to the burning end, 
behaves quite differently from an uncooled carbon 
which is above incandescent temperature for at least 
one inch behind the burning end. This change in 
behavior seems to be caused by structural changes 
in the carbon due to the heating which are prevented 
by effective cooling. Such structural changes have 
been actually found by the author by x-ray diffrac- 
tion studies which are as yet unpublished. The ex- 
planation of the cooling effect derived from these 
studies is as follows: Although the decrease of the 
anode drop which results from cooling causes a de- 
crease of consumption and brilliancy for a given 
current density, the unsoftened solid structure of 
the cooled carbon allows operation at much higher 
load with consequent higher brilliancy. The over- 
compensating effect of an increased current density 
as against the moderate anode drop accounts for the 
surprisingly low positive carbon consumption per 
unit brilliancy. 


THEORY OF THE ANODE DROP FOR LOW AND HIGH 
CURRENT CARBON ARCS 


The foregoing sections have made evident the 
important role which the anode drop seems to play 
in the mechanism of the high current carbon arc. It 
must be possible, in principle at least, to develop a 
mathematical theory of the anode drop based on the 
anodic mechanism for the low and high current 
carbon arc described above. The difficulties en- 
countered in this task seem to be of general interest 
in the field of arc physics. 

An attempt at setting up a quantitative anode 
drop theory has been made by Heinzmann and the 
author.'® Under certain simplifying conditions (mon- 
otonic change of the potential, velocity of electrons 
and ions given by the product of mobility and local 
field strength, field strength in the adjacent arc 
stream negligible compared with that in the anode 
drop, diffusion negligible) the distribution of po- 
tential, electric field strength, positive and negative 
space charge, and positive and negative current 
density could be calculated. In addition, a formula 
has been derived for the thickness of the anode drop 
region ; 

d=(12b-V.?/rj)}, 


in which b~ is the electron mobility, V, the anode 
drop, and j the anodic current density. 

The only difference between the systems of equa- 
tions for the low current arc and the high current 


16 W. Finkelnburg and G. Heinzmann, unpublished descrip- 
tion in reference 1. 
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carbon arc is the introduction of the vapor stream 
velocity into the equations of the latter. We believe, 
by the way, that our anode drop theory for the low 
current arc, and with it our formula for its thickness, 
should apply for discharges quite generally, e.g., for 
glow discharges where it could be checked rather 
simply. 

The anode drop itself could not be calculated but 
entered the equations as an empirical boundary 
condition. The reason for this deficiency of the 
theory is that we do not have an expression for the 
ionization in a high temperature discharge as a 
function of potential or local field strength in the 
regions of anode and cathode drop where at one end 
the thermal velocity of the electrons is larger than 
their field velocity, whereas the case is reversed at 
its other end. 


OPEN PROBLEMS AND CONCLUSIONS FOR HIGH 
TEMPERATURE PHYSICS 


In the last section we mentioned the unsolved 
problem of the ionization which is essentially a 
problem of high temperature physics. The same 
applies to the problem of how much the ionization 
potential of any atom or molecule is lowered owing 
to the presence of the large density of electrons and 
ions in the surrounding plasma, which, furthermore, 
influence the mobility of electrons and ions in a way 
not yet accurately known. A number of other details 
in the behavior of the contracted high current arc 
stream are still unexplained. How drastically the 
properties of an ordinary gas like air are changed by 
temperatures above 10,000°K becomes evident 
when it is realized that the electrical conductivity in 
such gases is of the same order as that of some 
metals and that the optical behavior (predominantly 
continuous emission and absorption) has little in 
common with that of ‘‘normal’’ gases. 

These few remarks may suffice to indicate what 
interesting new insights may result from further 
study of the high temperature behavior of matter.'® 
We think that the high current carbon arc is ideally 
suited for such studies because (a) it can be operated 
in any atmosphere at any desired pressure; (b) any 
kind of solid material can be introduced into the 
positive carbon core and thus vaporized and heated 
to plasma temperatures; (c) the arc in all its parts is 
easily accessible for detail studies. 

The author is greatly indebted to Messrs. M. K. 
Laufer and G. W. Franks for their help in the final 
presentation of this paper. 


16 W. Finkelnburg, Naturwiss. 32, 105 (1944) ; Die Chemische 
Technik 15, 141 (1942). 
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Radio Doppler Investigation of Meteoric Heights and Velocities* 


L. A. MANNING, O. G. VILLARD, Jr., AND A. M. PETERSON 
Electronics Research Laboratory of the Department of Electrical Engineering, Stanford University, Stanford, California 
(Received October 1, 1948) 


A continuous-wave Doppler method of determining meteoric velocities and heights was tested 
during the 1948 Perseid meteor shower. The accuracy of the technique was found to compare very 
favorably with optical and other radio methods. Doppler measurements have the advantage of 


relatively great sensitivity. 





T the present time three methods are known for 
obtaining meteoric velocities directly from 
radio reflections returned by meteor ionization 
trails. The first to be described was that of Hey.! 


Working with the slow meteors of the Giacobini- ° 


Zinner shower of October 9, 1946, he obtained an 
average experimental velocity of 22.9 km/sec. with 
a standard deviation of 1.3 km/sec. The theoretical 
velocity is 23.7 km/sec. Hey’s method consisted of a 
measurement of the change in range of the pulse 
echo as the meteor came in. He used a wave-length 
of five meters, and based his results upon 12 obser- 
vations. Ellyett and Davies? have determined 
meteoric velocities by a second method based upon 
a measurement of the amplitude changes in pulse 
echoes from a filamentary meteor trail as the trail 
becomes extended in length. From observations of 
the Geminids, they deduced a mean velocity of 34.4 
km/sec. with standard deviation of 2.8 km/sec., 
which they compared with Whipple’s* photographic 
mean geocentric velocity of 34.7 km/sec. The third 
published method of determining meteoric velocities 
utilizes the Doppler frequency shift imparted to a 
continuous-wave signal by the motion of the re- 
flecting cylinder of ionization.‘ The present paper 
describes a test of this technique. Velocities of the 
meteors of the August, 1948 Perseid shower were 
measured and the results found to be in good 
agreement with those obtained by optical means. 


EXPERIMENTAL ARRANGEMENTS 


For this test continuous-wave transmitters were 
operated on 23.1 and 30.66 Mc, as well as a pulse 
transmitter on 17 Mc for range determinations. The 
radiated power of the continuous-wave transmissions 
was of the order of magnitude of a kilowatt, as was 
the peak power of the 150-microsecond pulses. Both 
transmitting and receiving antennas consisted of 
half-wave dipoles, one-quarter wave high. The re- 


* The research described in this paper was supported by the 
ONR and the US Army Signal Corps. : 

1 J. S. Hey, ‘““The Giacobinid meteor shower, 1946,” Nature 
159, 119-121 (1947). ; 

?C. P. Ellyett and J. G. Davies, ‘Velocity of meteors 
measured by diffraction of radio waves from trails during 
formation,” Nature 161, 596-597 (1948). 

’F. L. Whipple, Proc. Am. Phil. Soc. 91, 189 (1947). 

*L. A. Manning, “The theory of the radio detection of 
meteors,” J. App. Phys. 19, 689-699 (1948). 
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ceivers were located in a small valley about four 
miles from the transmitters. For continuous-wave 
reception, the strength of the ground wave was kept 
within bounds by the sheltered nature of the valley, 
the orientation of the receiving dipoles, and by out- 
of-phase cancellation of a portion of the ground 
wave by means of a signal from directional antennas 
pointed towards the transmitters and so arranged as 
to be insensitive to signals arriving at angles much 
above the horizontal. Although a _ near-perfect 
ground-wave null is obtainable by the above 
methods, some ground-wave signal was allowed to 
leak through in order to beat with the reflected 
meteor signals. 

The output of the continuous-wave receivers was 
recorded in two ways. A pair of d.c. amplifiers fed 
directly from the receiver second detectors were used 
to drive a two-channel direct-inking magnetic oscil- 
lograph. This oscillograph, employing Brush ‘‘pen- 
motor’ units, has good response from zero to 120 
c.p.s., and was operated at a paper speed of 125 
mm/sec. It thus indicated all slow- and medium- 
speed variations in received field strength, and 
enabled a side-by-side comparison of the perform- 
ance of the two channels. Doppler beats of higher 
frequency were recorded by means of a pair of 
magnetic tape recorders. For later analysis, these 
‘whistles’ were spliced out and played back into a 
standard two-channel oscillograph having a film 
speed of 48 inches per second. Thus, all Doppler 
shifts from zero to about 1500 c.p.s., were recorded 
for study. 

The output of the pulse receiver was displayed on 
an “A ’scope”’ and slant ranges were read visually 
with the aid of ten-kilometer timing marks. The 
accuracy of range determination was of the order of 
plus or minus five kilometers for typical echoes. 
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Fic. 1. Typical low frequency Doppler oscillogram. 
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Fic. 2. Typical high frequency Doppler oscillogram. 


Visual meteor observations were carried out on 
the night of the test with the aid of a team of ob- 
servers facing the four quadrants of the sky. In 
order to obtain clear seeing conditions, the visual 
observations were made at an elevation of about 
2000 feet in the Santa Cruz mountain range about 
three miles from the receivers. The bearings of the 
observed meteors were relayed to the receiving site 
by radio, for correlation with received echoes. 


NATURE OF THE REFLECTIONS 


The general nature of the phenomenon giving rise 
to Doppler meteor records has been described in a 
previous paper.‘ In Fig. 1 is shown a well-behaved 
penmotor field-strength record illustrating the prin- 
cipal effects of present interest. The ‘down whistle”’ 
with which the echo commences is associated with 
the component of meteoric velocity towards the ob- 
server, and falls below the audible range when the 
meteor arrives at that part of its straightline path 
nearest the observatory. It is this down whistle 
which is analyzed for velocity determination. In the 
case of large meteors, a non-linearity of the whistle- 
frequency versus time curve is sometimes observed 
at the lower whistle frequencies as a result of ex- 
pansion of the ionization trail. The method of 
analysis of this non-linearity to determine trail size 
is a topic for a separate paper; in determining 
velocities, however, it is only necessary to restrict 
one-self to the high pitched part of the whistle in 
order to avoid inaccuracies due to the rate of trail 
formation. 

A large increase in signal strength follows the 
down whistle because the incident wave is then re- 
flected normal to the cylinder of ionization. It is 
this maximum signal, or ‘“‘burst,’’ that is seen on the 
A.'scope when determining range. The “body 
Doppler” which exists during the burst, or “‘body”’ 
of the echo, is of very much lower frequency than 
the down whistle and is of opposite sign, indicating 
a receding reflector. It is apparently caused by 
radial contraction of the fully formed column, again 
a subject for a separate report. The signal finally 
dies out when the ionization is dissipated by recom- 
bination and/or diffusion. 

Figure 2 is a good example of a whistle oscillogram 
as recorded via the magnetic tape. This represents 
that portion of the down whistle which is too high in 
frequency to be recorded by the penmotor of Fig. 1. 
Unfortunately, the wave form is somewhat distorted 
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by the magnetic recording process, but for the 
purposes of frequency analysis, wave form is non- 
essential. Although velocities of slow meteors can 
frequently be obtained from the portion of the 
whistle appearing on direct-inking oscillograms as 
in Fig. 1, fast meteors such as the Perseids have so 
high a rate of change of Doppler frequency with 
time, that an insufficient number of cycles is ob- 
tained within the frequency response of the instru- 
ment. It is then necessary to use magnetic tape 
oscillograms which show the higher whistle fre- 
quencies, as in Fig. 2. 


THE PERSEID DATA 


The measurements described here were made on 
the night of August 11-12, 1948, starting at ap- 
proximately 10 p.m., P.S.T., and concluding at 
about 3 a.m. Visual observations were made in the 
interval from 11:50 p.m. to 2:30 a.m. That an 
abundant shower was in progress is apparent from 
the rate of visual detection presented in Fig. 3. A 
very large number of meteors was recorded by the 
radio equipment, but only those for which either an 
audible whistle or a visual azimuth were obtained, 
which could be definitely correlated with a pulse 
range, will be discussed here. A total of 138 whistles 
were catalogued, and visual azimuths were assigned 
to 36 of that number. Of the 138 meteors with 
recorded whistles, it was possible to determine 
velocities for 60; the remaining whistles were too 
weak to reduce or were otherwise inadequate. Only 
14 of the 138 were analyzed on both 23 and 30 Mc, 
because of an unfortunate accident in the processing 
of the 23-Mc records. Velocities determined at the 
two frequencies were in satisfactory agreement. In 
the group of 36 whistles of meteors with azimuths, 
20 were reducible to velocities. 
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Fic. 3. Rate of visual detection—August 11, 1948. 
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Fic. 4. Typical whistle-pitch versus time plots. 


THE VELOCITIES 


In a previous paper‘ one of the authors has shown 
that if the frequency, f, of a down whistle can be de- 
termined as a function of time, then a relation be- 
tween the velocity and range of the meteor is 
specified. Denoting the slope with which the whistle 
pitch curve approaches zero frequency by fo’, the 
velocity of the meteor is 


V =(—fo'AR/2)!, 


where \ is the wave-length in the medium, and R is 
the pulse range. In Fig. 4 are shown two experi- 
mental curves of whistles pitch versus time. The first 
curve is linear and presents no questions in inter- 
pretation. The second curve is linear at the higher 
frequencies, but does not drop to zero as quickly as 
must the radial velocity of a meteor moving in a 
straight line with constant speed. The additional 
component of Doppler shift present at low fre- 
quencies is caused by radial expansion of the ioniza- 
tion column. The relation between this expansion 
and the Doppler shift has been studied theoretically 
and it has been found that the slope of a straight line 
drawn as shown will give the correct meteoric 
velocity despite diffusive expansion. 

Velocities were determined for 60 meteors using 
the technique just described. In Fig. 5 is shown a 
plot of the resulting data, less six meteors for which 
either the range or the slope was in doubt. From the 
complete grouping of Fig. 5 it is possible to form 
two sub-groups with the aid of the visually-observed 
meteoric azimuths. Since we are now dealing ex- 
clusively with meteors producing a Doppler whistle, 
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we know that the ray from the observatory is per- 
pendicular to the direction of meteoric travel. The 
geometry of Fig. 6 must then apply to the reflection 
process. For each meteor it is possible to construct a 
line, L, seen in the views, which represents the 
intersection of a plane, P, perpendicular to the 
meteor path and containing the observatory, with a 
plane, M, at the height, 4, of meteoric reflection. All 
shower meteors must produce their reflection when 
passing through L, so the difference between their 
observed azimuths A and the azimuth of the shower 
radiant Z must always be between 90 and 270°. In 
addition, a relation is specified between (A —Z), h, 
and the pulse range R, if a meteor is to be along the 
line L determined by a radiant of elevation, H. 

If we segregate those whistles having azimuth 
data into two groups as definitely Perseids or non- 
Perseids depending upon whether they satisfy the 
condition of the preceding paragraph, the distribu- 
tion curves of Fig. 7 are obtained. It is strikingly 
evident in the case of these identified meteors that 
the Doppler velocities truly represent the shower 
speed. For the definite Perseids, the average velocity 
determined by the Doppler technique is 62.3 
km/sec., with a standard deviation of 1.6 km/sec. 
There is one photographically-determined Perseid 
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Fic. 5. Meteor velocities—August 11, 1948 (average of 
Perseid meteors: 62.3 km/sec.). 
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velocity available for comparison.’ Whipple found a 
velocity of 61.19 km/sec. with a deceleration of ap- 
proximately 10 km/sec.?. The accuracy of the 
photographic velocities is probably one percent. 

In plotting Fig. 5 the definitely non-Perseid 
meteors of Fig. 7 have been cross-hatched with 
negative slope. It is seen that there remain at 
velocities between 56 and 69 km/sec. a single well- 
defined cluster. It is not at all certain that this 
cluster contains no non-Perseid meteors, but if we 
assume them to be Perseids we may compute an 
average velocity of 63.2 km/sec. with a standard 
deviation of 2.5 km/sec. 


ACCURACY OF THE VELOCITIES 


The accuracy of the Doppler velocities determined 
for the Perseids may be expressed in terms of 
standard deviation as a percentage of the average 
velocity. For the select group of definite Perseids, 
the accuracy figure becomes 2.6 percent, while for 
the group of presumed Perseids, the accuracy figure 
is 4.0 percent. By way of comparison, the same 
figure for Hey’s determination of Giacobinid ve- 
locities' is 5.5 percent, and for Ellyett and Davies’ 
Geminid observations, 8.1 percent. It should be 
pointed out that the Perseid meteors, because of 
their high speed, provide an especially severe test of 
radio methods of velocity determination in the 
following respect. If errors are to arise because of 
such effects as an ‘ionic bow wave’’® lagging behind 
the particle, the effect should be most pronounced 
for fast meteors, because for them the ionization 
contours must have the greatest difficulty in keeping 
up. 

Consideration of the cause of the observed scatter 
must begin with the accuracy of range determina- 
tion. Ten-kilometer range marks were used, so that 
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Fic. 6. Geometry of reflections from shower meteors. 


5 F. L. Whipple, ‘Meteors and the earth’s upper atmos- 
phere,” Rev. Mod. Phys. 15, 246-264 (1943). 

6J. A. Pierce, “Ionization by meteoric bombardment,” 
Phys. Rev. 71, 88-92 (1947). 
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a perfect observer, reading range to within 10 km, 
would have an average error of 2.5 km. Assuming 
the actual accuracy of range determination to aver- 
age 5 km, and assuming the average range as 130 
km, the average error in the velocities would be 
about 1.2 km. As much or more error was probably 
introduced in finding the slope of the whistle pitch 
curves. Error of the latter sort can be reduced by 
careful attention to the quality of the whistle wave 
forms. 


HEIGHT OF THE PERSEIDS 


Use can be made of the perpendicular nature of 
the reflection for meteors with whistles, to determine 
altitudes with good precision. Consideration of 
Fig. 6, or preferable its three-dimensional counter- 
part, shows that the pulse range may be represented 
in terms of its rectangular components as: 


®=(h tanH)?+[h tanH tan(A—Z) +h’, (2) 
so that solving for the height h, 
R 


h= , 
(1+tan?*H-sec?(A —Z))! 





(3) 


The azimuth Z and elevation angle H of the meteoric 
radiant are very accurately known for shower 
meteors. The azimuth A can be observed visually, 
and if those meteors are chosen for which sec(A —Z) 
is near 1, the accuracy of the visual bearings be- 
comes almost immaterial. Here is a simple method 
for determining the heights of shower meteors, in 
which the accuracy of height determination can 
approach that of pulse range measurement. As an 
illustration, suppose a meteor is chosen with an 
indicated angle (A —Z) of 0°, and a radiant eleva- 
tion H of 30°. Then a 5° error in azimuth determina- 
tion introduces a mere 1/20 of one percent error in 
altitude. | 

In contrast, if altitude is determined by noting 
the pulse range and visual elevation of the meteor, 
errors in angle measurement are accentuated in 
computing heights, and uncertainty as to the part of 
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Fic. 7. Velocities of identified Perseid and 
non-Perseid meteors. 
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the trail responsible for radio reflection introduces 
additional inaccuracy. In the present method the 
existence of the whistle is, of course, essential in that 
it guarantees the correctness of the geometrical 
construction. 

From the group of meteors having both whistles 
and visually determined azimuths, thirteen Perseids 
were selected, and their heights as computed from 
Eq. (3) are presented in Fig. 8. The range is princi- 
pally in the 90-110-km region, but a tendency for 
the altitudes to decrease with time is evident. Such a 
behavior might be anticipated from the behavior of 
the radiant, since its elevation angle was steadily 
increasing from a low value during the period of 
observation. 


CONCLUSION 


In this test the Doppler method of determining 
meteoric velocities has demonstrated itself to be 
remarkably accurate and reliable, deviations in 
measured velocities of only a few percent being ob- 
served for the Perseid meteors. The method has the 
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Fic. 8. Altitudes of identified Perseid shower meteors. 


advantage of relatively great sensitivity in compari- 
son with other means of finding meteoric velocities. 
Study of shower meteors producing Doppler whistles 
has also been shown to be an exceptionally accurate 
way to obtain the altitudes of meteoric ionization 
columns. 
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Improved Sectioning Technique 
for Electron Microscopy 
RICHARD F. BAKER AND DANIEL C. PEASE 


Departments of Experimental Medicine and Anatomy, School of Medicine, 
University of Southern California, Los Angeles, California 


(April 27, 1949) 


METHOD of cutting sections with an ordinary micro- 
tome for electron microscopy was published about a year 
ago.! Sections 0.2-micron thick were described originally. 
Since that time, however, refinements have been made which 
make possible cutting and mounting 0.05-micron sections. It 
seems worth while to comment upon these improvements now, 
pending the publication of a fully detailed account. 
Successful ultra-thin sectioning with a conventional micro- 
tome depends fundamentally upon three factors. 1. Suitable 
support must be given the specimen. 2. The knife must be 
sharp. 3. The adjustment of the knife tilt is critical. Also, 
since no commercial microtome is geared to cut fractional 
micron sections, an adujstment of the lever system is neces- 
sary. A very simply constructed adaptor for the Spencer No. 
820 rotary microtome was described originally, designed so 
that the block advanced 0.1 micron with each cycle. The angle. 
of lever action has been further reduced so that 0.05-micron 
sections now can be cut. 








Fic. Electron micrograph of transverse section of doubly embedded 
hair. A overlapping fibers are shown. Insert at higher magnification 
shows complex structure around edge. Section thickness apprx. 0.2 micron; 
50 kv EMU microscope, 
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The original method of embedding specimens has not been 
significantly improved. After suitable fixation, biological 
material is dehydrated in a graded series of alcohols and then 
infiltrated with Parlodion (Mallinckrodt) dissolved in equal 
parts of ether and alcohol. The maximum concentration used 
is 12 percent. The Parlodion is then hardened in chloroform. 
It actually is precipitated as a very fine three-dimensional 
network. The block is then transferred to molten paraffin 
(m.p. 65°C). This replaces the chloroform and fills the inter- 
stices. When cooled the block is ready for sectioning. The 
Parlodion has sufficient elasticity and toughness to hold the 
thin sections together; the paraffin supplies the requisite 
hardness for cutting. 

The knives used are standard microtome blades, and satis- 
factory edges can be obtained with hand sharpening. Our ex- 
perience has been that a sharp edge will deteriorate within a 
week or so whether or not it is used. Accordingly it is efficient 
to employ automatic sharpening, and we now use a modified 
Fanz automatic knife sharpener (Arthur H. Thomas Com- 
pany). Grinding with white rouge and soap is followed by 
polishing with plain water. 

Very critical adjustment of the knife tilt is essential. The 
tilt is considered optimum when a section is cut with every 
cycle of the microtome. Outside the range of optimum tilt, 
which is not more than 0.5°, skipping and shredding of sec- 
tions takes place. At the suggestion of Dr. J. Hillier an optical 
lever is used to measure the tilt. A horizontal light beam is re- 
flected from the back facet of the knife bevel to a scale. The 
optimum knife position usually corresponds to an angle of 
about 2° between this facet and the block. 

The most important new development concerns the han- 
dling of the sections once they are cut. Sections were originally 
flattened with needles by hand. By using dioxane the need for 
manual dexterity has been nearly eliminated. Sections are 
transferred from the knife edge with a small brush and placed 
on the surface of a drop of 50 percent dioxane in water. This 
begins to soften the embedding material, and surface tension 
starts to flatten the section. Then pure dioxane is added to the 
original drop. As the concentration increases, both the paraffin 
and the Parlodion are softened and partly removed. Inter- 
facial tension completes the flattening. After the section has 
dried on the glass, residual paraffin is taken out with benzene, 
and finally the slide is dipped vertically into a 0.2 percent 
Parlodion in amy] acetate solution. This dissolves the original 
Parlodion and, after the slide has been drained, the sections 
are left embedded in a homogeneous film of controlled thick- 
ness. Screens are placed over the sections and the film floated 
off the slide onto a clean water surface. After picking up and 
drying they are ready for use. 

Sectioning is normally done at room temperature. The 
sectioned bacteria of the cover illustration were produced in 
the manner described. In the case of very hard material, 
sectioning may be facilitated by lowering the temperature 
of the block. This increases the hardness of the embedding 
material so that more support is offered. We have used dry ice 
for this purpose, and with its aid have sectioned hair (Fig. 1), 
which otherwise ripped from the support. It is to be empha- 
sized in this connection that although cooling always im- 
proves the cutting qualities of the block, it is never a substitute 
for a sharp knife properly tilted. 

The thickness of the sections has been measured using a 
pre-shadowed replica technique. The sections were flattened 
on glass and completely extracted before being shadow-cast. 
Dow polystyrene latex spheres of known size were scattered 
around the section to serve as an internal standard. The 
thickness measured in this way is in satisfactory agreement 
with the value expected from the microtome calibration. 


1D. C. Pease and R. F. 


Baker, Proc. Soc. Exper. Biol. and Med. 67, 
470-474 (1948). 
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